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MODIFICATION  OF  PROGRAMS  CHIEF  AND  CID  TO  COMPUTE 
SOUND  RADIATION  FROM  AN  ARBITRARY  BODY  IHH  IfTTRD 
BOUNDARY  CONDITIONS 


INTRODUCTION 

The  computer  program  CHIEF  [1] ,  written  to  compute  the  sound  radiation 
from  ?n  arbitrarily  shaped  body,  was  originally  developed  to  handle  only  the 
case  of  Neumann  boundary  conditions  on  the  surface  of  the  body;  i.e.,  the  case 
where  the  normal  velocity  is  prescribed  over  the  entire  surface. 

There  may  be  situations  where  it  is  desired  to  know  the  effect  of  putting 
some  material  on  portions  of  the  surface  of  the  radiator.  For  example,  it  may 
be  interesting  to  see  if  the  shape  or  side-lobe  structure  of  the  directivity 
can  be  favorably  changed  or  if  such  a  material  has  a  deleterious  effect  upon 
the  radiation  efficiency  of  the  source.  A  case  in  point  is  that  low- impedeuice 
material  (pressure-release  material)  is  often  applied  to  the  surface  of  an 
underwater  structure  to  effect  shielding  of  that  part  from  the  acoustic  field. 

From  a  mathematical  point  of  view,  this  class  of  boundary  value  problems 
characterized  by  mixed  boundary  conditions  [where  at  least  two  different  types 
of  boundary  conditions  (Neumann,  Dirchlet,  or  Robin),  also  called  impedance, 
are  prescribed  on  the  same  surface]  is  not  amenable  to  exact  solution.  Even 
for  a  simple  coordinate  surface  such  as  a  perfect  sphere,  the  mix  of  boundary 
conditions  means  that  the  solution  functions  of  the  Helmholtz  equation  do  not 
form  a  complete  orthogonal  basis  set  over  the  surface.  Hence,  the  classic 
solution  technique  of  expanding  the  boundary  condition  into  a  generalized 
Fourier  Series  url  then  matching  mode  for  mode  with  a  similar  series 
expression  for  the  radiation  field  is  not  applicable.  Approximate  matching 
procedures  such  as  point  collation  and  lezist-square  error  determination  have 
been  tried  and  are  discussed  in  the  literature.  [2,3] 

This  report  discusses  modifications  that  have  been  made  to  the  programs 
CHIEF  and  CID  to  accommodate  the  situation  where  part  or  parts  of  the  radiator 
are  covered  with  one  or  more  locally  reacting  arbitrary  complex  impedance 
coatings.  The  range  of  magnitude  of  these  impedances  can  be  from  zero 
(perfect  pressure  release)  to  the  largest  value  the  computer  can  handle 
without  exponent  overflow  (rigid);  any  value  that  is  purely  real,  purely 
imaginary  or  complex,  in  between  these  extremes  is  allowable.  Three  sample 
problems  are  presented  to  illustrate  many  of  the  new  improvements  to  both  CID 
and  CHIEF. 


NEW  FEATURES  AND  MODIFICATIONS  TO  CHIEF 

The  following  new  features  have  been  added  to  the  1991  version  of  CHIEF. 

1.  CHIEF  now  includes  the  capability  to  describe  multiple  locally 
reacting  complex  impedances. 

2.  The  cone  geometry  has  been  added  to  CCUNMD. 
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3.  The  routine  ZRADMl  has  been  added  to  calculate  the  radiation  matrix 
when  impedance  surfaces  are  defined. 

4.  New  storage  allocation  (Units  26  through  29)  have  been  assigned. 


NEW  miURSS  AND  MODlFICAnONS  TO  CID 

The  following  new  features  have  been  added  the  1991  version  of  CID. 

1 .  The  cone  geometry  has  been  added  to  the  SURFACE  DESCRIPTION  screen 
(see  Appendix  Bl) . 

2.  The  impedance  surface  description  screens  have  been  added  (see 
Appendix  B2) . 

3.  The  DISPLAY3D  [4]  program  used  to  display  the  CHIEF  geometries  uses 
the  CID  screen  input. 


USSR  INPUTS  AND  SUBROUTINB  GALLS 

This  section  describes  the  changes  that  have  been  made  to  the  input 
parameters  and  subroutine  calls  since  the  last  program  release. 


General  Input  Parameters  for  Impedance  Surfaces 

NUUZ  Nxuaber  of  surface  regions  in  a  symmetry  block  that  are  covered 

with  the  impedance  coating.  (Note:  The  NUHZ  impedance  surfaces 
appear  after  the  radiating  surfaces.) 

ZCOAT  Vector  that  contains  the  complex  impedance  coating  for  each 
surface  region  in  a  symmetry  block. 


Calculation  of  Radiation  Matrix 
ZRADMI(ZMTX,NUMARS) 

ZMTX  Complex  matrix  area  that  must  be  allocated  in  control  routine 
for  storage  of  the  radiation  impedance  matrix  (NUMARS,NUMARS) . 

NUMARS  Number  of  surface  areas  (including  impedance  surfaces) . 

SAMPLE  PROBLEMS 

This  section  contains  the  CHIEF  drivers  for  various  types  of  sample 
problems .  Since  many  of  the  dimension  statements  are  the  same  for  each 
control  program,  these  lines  will  be  printed  only  once  and  are  labeled 
VARIABLE  DECLARATION. 
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Conon  Variable  Declaration 


C  *****  CONTROL  88  ***** 

C  PROGRAM  CHie^88  DRIVER 

C  MXSRBQ  -  MAXIMUM  NUMBER  OF  SURFACE  REGIONS 

C  MXIPS  -  MAXIMUM  NUMBER  OF  INTERIOR  POINTS 

C  MXARS  -  MAXIMUM  NUMBER  OF  SURFACE  SUBDIVISIONS/SYM  BLK 

C  MXGAUS  -  MAXIMUM  ORDER  OF  GAUSSIAN  QUADRATURE 
C  MXQPTS  -  MAXIMUM  NUMBER  OF  QUADRATURE  POINTS/SUBDIVISION 

C  MXBLKS  -  MAXIMUM  NUMBER  OF  SYIUETRY  BLOCKS 

C  MXFFP  -  MAXIMUM  NUMBER  OF  FAR-FiaD  POINTS 

C  MXNFP  -  MAXIMUM  NUMBER  OF  NEAR-FIB.0  POINTS 

C  MXPTSC  -  MAXIMUM  NUMBER  OF  POINT  SOURCES 

C  MAXCOR  -  MAXIMUM  NUMBER  OF  FINITE  B-QIENT  NODES 

C  MXFPS  -  MAXO(MXARS*MXIPS, MXFFP, MXNFP) 

PARAMETER  (MXSREG=8»e) 

PARAMETER  (MAXCOR=iae0) 

PARAMETER  (MxiPS=2l) 

PARAMETER  (MXARS=SM) 

PARAMETER  (MXGAUS=84) 

PARAMETER  (MXQPTS^SIO) 

PARAMETER  (MXBLKS=1M) 

PARAMETER  (MXFFP=361) 

PARAMETER  (MXNFP=381) 

PARAMETER  (MXPTSC=2<) 

PARAMETER  (MXFPS=62I] 

C  PARAMETER  (MXFPS=MAXe(MXARS«MXIPS, MXFFP, MXNFP)} 

C  PARAMETER  (NfDVEC=2*MXARS) 

C 

C  Input  coaaons 
C 

COMMON  'CONST/  RH0,C 
COMMON  /PRTCOM/  NUNPRT , NUNERR 
COMON  /PRTRD/  RUNIO,DATE 
CHARACTER*4  RUNID 
.■HARACTER*8  DATE 

COMMON  /K  iA  G/  NOQPTS,NOPMXS,NDVMXS,NDDBaK,NOVB-S,NDSPS. 

1  NDPMXF,NOVMXF,NDPMXN,NDVMXN,NDPSSP,NDEXPK,NDCOMV, 

1  NOTaiP,NOZROe,NOPATB 

COMMON  /SVALS/  NSREQ,NSBQNS(MXSREQ),SUL(MXSREG),SUU(MXSReQ), 

1  SVL(MXSRB3),SVU(MXSREG},NSU(MXSREQ),NSV(MXSREG), 

1  CCS(ia,MXSReQ),TRNSS(3,MXSREQ},IZAX(MXSReQ}, 

I  IOROU(MXSREO),IOROV(MXSi^},NCCEQS 

COMMN  /CORD/  COOROS(MAXCOR,3} 

COUiON  /IPTS/  NUMIPS,IPXS(3, MXIPS) 

REAL  IPXS 

COkMON  /PTSINP/  NUMPTS , PTSRCS (4 , MXPTSC) , PTSIT (MXPTSC) , 

1  lOPTSC (MXPTSC) 

COMPLEX  PTSIT 

COMMON  /PLfINP/  AINC,THTINC,PHIINC,ISCATR 

COMMON  /BAFFLE/  INFBAF 

COMMON  /FFINP/  NUMTHP,THTPHI(2,MXFFP) 

COkMON  /NFINP/  NUkFPN,NFPXS(3,MXNFP) 

REAL  NFPXS 
C 

C  Output  coaaons 
C 

COMMON  /TAPREC/  RECR0(ia),IRECR0(3a) 

COMMON  /TAPRCl/  ARECRD(ia) 

CHARACTER#4  ARBCRD 

COMMON  /PRGVLS/  NDIMPV,NUMARS,NUMSFP,NUIFFP,NUMNFP,NfDVEC 
COMMON  /SURARS/  AREAS (MXARS) 

COMMON  /ODSVEC/  TVECT (MXARS) ,B (MXARS), IPIVTR(kCffPS) 

COMPLEX  TVKT 

COMMON  /VB^PS/  Va (MXARS), SP (MXARS) 

COMPLEX  va,sp 

COMMON  /PDISL/  POWER, DIRIND, SRCRVL 

COMMON  /FFVALS/  FFP (MXFFP), PNRMFF (MXFFP), IFFNRM,RkrNRM 

COMPLEX  FFP 

COMMON  /TSCOM/  TGTSTH (MXFFP) 

COMMON  /NFVALS/  NFP (MXNFP) ,PNRMNF (MXNFP), INFNRM,RMNNRM 
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COyPLEX  NFP 

COMMIN  /PTSCOM/  PTSSP(MXA;S} 

COyPLEX  PTSSP 

COIMON  /EXTCOM/  E)rPRS(MXFPS)  .I0CTFG 
COMPLEX  EXTPRS 

COIMON  /NBPRTC/  IRHSPT.NARSPT.NPTBLK.FRRPT 
COtHON  /NBPRTS/  SYHTPT 
CHARACTER»3  SYMTPT 


C***************  iip«dance  coating  aodif ications  «««•***«•***• 
COIMON  /COATING/  ZCQAT,  NUMZ 
COMPLEX  ZC0AT(5e«) 

C  NUMZ  -  total  nuiber  of  iapadanca  layar  surfacas  (N!>(.'  *  NSV) 

C**«**«*««««*«t««««*«*«««*««4«M*«*«**«««M*«««*«**»*««««M»** 


DIMENSION  CC(10),  TRNS(3h  IB.TS(a,  306) 

REAL  XI (1000),  Yl(1000) 

CHARACTER»3  SYMTYP 

CHARACTER*4  FLOTYP,  TAPEIO,  PRHYP 

INTEGER  XIZAX,  XNSEQNS,  XIRG,  XNSU,  XNSV,  XIORDO,  XIORDV 

INTEGER  ICOOR,  18.01 
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Saaple  Run  1:  Riulinting  Polnr  Cap  on  an  lapedance  Coated  Spherical  Baffle 

This  sample  problem  computes  the  acoustic  radiation  from  a  spherical 
source  consisting  of  a  90°  uniformly  vibrating  polar  cap  while  the  rest  of  the 
sphere  is  covered  with  an  impedance  layer.  Because  of  axisymmetry,  the 
azimuthal  coordim.te  angle  was  taken  to  be  1/30  of  the  sphere.  There  are  24 
regions;  the  first  twelve  regions  (90  )  define  the  polar  cap.  Thus,  there  are 
12  regions  with  an  impedance  coated  material . 

Three  cases  were  run  using  this  break-up.  In  the  first  case,  the 
impedance  layer  is  defined  as  a  perfect  pressure-release  material  (ZCOAT(I)^  = 
(0,0) .  In  case  2,  the  impedance  coating  is  defined  as  having  pc 
characteristics  (ZCOAT(I))  =  (/Jc,0).  And  in  case  3,  the  coating  is  defined  as 
19 

(ZCOAT(I))  =  (10  ,0),  which  simulates  a  rigid  boundary  condition.  Setting 

19 

ZCOAT  equal  to  10  (real  or  imaginary)  produces  far  field  results  identical  to 
those  generated  using  the  CHIEF  where  the  velocity  was  specified  to  be 
identically  zero  on  the  impedance  coated  regions. 

GEOMETRY  OF  THE  PROBLEM 


Fig.  1.  B0  uniformlj  Tibrsting  polar  cap  on  an  ispedanca-coated  apbarical  baffla. 

Program  Listing  of  Sample  Run  1 

C  Prograa  MRBILLT 

C  SPHERE  IITH  IMPEDANCE  COATINQ  SURFACE 

C  Variable  Declarations 

C  (COMMON  VARIABLE  DECLARATIONS  GO  HERE) 
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COtHSX  PUATX(  1752),  VMATX(  1752) 

MK  .ZE  =  1752 
RUhIO  =  ’URBl’ 

DATE  =  ’19-FEB-91’ 

CALL  INITCM 
rUl  DPNSFL 
RHO  =  leei 

C  =  ISM 

OPEN (UNIT=Nt’NPRT, FILE=RUNIO//  ’  OUT ' , aTATUS=  ’ NBI ’ 

1  F0RU=' FORMATTED') 

i«pedanc«  coating  andif ications  »»*♦•*»»»♦•*»»» 
NUMZ  =  12 

N  =  1 


C  -  load  iapodanca  for  SPHERE 


ZC0AT(N  .  •)  =  CUPLX(l.f 

,  0  0 

) 

ZCOAT^H  ♦  1)  =  C«PLX(«.f 

.  0  0 

) 

2C0AT(N  -  2)  =  CMPLX(».« 

,  0  0 

) 

ZC0AT(N  .  3)  =  CMPLX(e.e 

.  0  0 

) 

ZC0AT(N  ♦  4)  =  CMPLX(«.6 

.  0  0 

) 

ZC0AT(N  ♦  6)  =  CMPLX(«  A 

,  0.0 

) 

ZC0AT(N  -  8)  =  CMPLX(B.0 

.  0  0 

) 

ZC0AT(N  >  7)  =  CMPLX(8.e 

.  0  0 

) 

ZC0AT(N  -  8)  =  CMPU((0.0 

,  0.0 

) 

ZCaAT(N  >  9)  =  CMPLX(0.0 

.  0  0 

1 

ZC0AT(N  ♦  10)  =  CMPLX(0.0 

.  0  0 

) 

ZCCAT(N  .  11)  =  CMPU((0.i 

.  00 

) 

DO  99  JJ=1,3 

N=1 

IF  (JJ  .£«.  1)  THEN 

ZC0AT(N  .  12)  =  CMPLX(0.6 

.  0.0 

) 

ZC0AT(N  ♦  13)  =  CUPLX(0.i 

.  0  0 

) 

ZC0AT(H  .  4  =  CMPLX(0.0 

,  0.0 

) 

2C0AT(N  .  16)  =  CMPLX(0.0 

.  0  0 

) 

ZC0AT(N  .  18)  =  CUPLX(i.0 

.  0  * 

) 

ZCOAIiii  •  17)  =  CMPU(0.0 

,  0  0 

) 

ZC0AT(N  .  18)  =  CMPLX(0.0 

.  0-0 

) 

ZC0AT(N  .  19)  =  CMPLX(0.0 

.  0  0 

) 

ZC0AT(N  .  20)  =  CMPLX(0.0 

.  0  0 

) 

ZC0AT(N  -  21)  =  CMPLX(0.0 

.  0  0 

) 

ZC0AT(N  .  22)  =  CHPLX(0.0 

,  0.0 

) 

ZC0AT(N  ♦  23)  =  CMPLX(0.0 

.  0  0 

) 

END  IF 

IF  (JJ  ,ER.  2)  THEN 

ZC0AT(N  .  12)  =  CMPLX(1.5E08 

,  0  0 

) 

ZC0AT(N  .  13)  =  CMPLX(1.5E06 

,  0  0 

) 

ZC0AT(H  -  14)  =  CHPL."(1.6EB6 

,  0  0 

) 

ZC0AT(N  -  16)  =  CMPLX(1.6E06 

.  0  0 

) 

ZC0AT(N  .  18)  =  CMPLX(1.63»8 

,  0.0 

) 

ZC0AT(N  .  17)  =  CMPLX(1.5E36 

,  0.0 

) 

ZC0AT(N  .  18)  =  CMPU((1.6E06 

,  0  0 

) 

ZC0AT(N  «  19)  =  CMPU{(1.6EB6 

,  0.0 

) 

ZC0AT(N  .  20)  =  CMPLX(1.5E08 

,  0  0 

) 

ZCOATfN  .  '1)  =  CMPLX(1,5E06 

,  0  0 

) 

ZC0AT(N  .  22)  =  CMPLX(1.6E06 

,  0.0 

) 

ZC0AT(N  •.  23)  =  CMPLX(1  6E06 

,  0.0 

) 

END  IF 

IF  (JJ  EQ,  3)  THEN 

ZCOAT'N  ■*  12)  =  CUPLX(0.0 

,  1.BE19 

) 

ZC0M{h  ♦  13)  =  CMPLX(0.0 

,  1  0E19 

) 

7CjAT(N  .  14)  =  C«PLX(0,0 

,  1.0E19 

i 

ZC1AT(N  «  16)  =  CUPLX(0.0 

,  1.0E19 

> 

ZCli’TfH  -  16)  =  CMPLX(0,0 

.  1.0E19 

) 

ZCDAT(N  *  17)  =  CUPLX(0.6 

,  1.0E19 

) 

ZCOATCM  -  18)  =  CUPLX(t.e 

,  1  0E19 

) 

ZCOAT(H  .  19)  =  CMPLX(B  0 

,  l.iElR 

' 

2C0AT(N  •  20)  =  CUPLX;B  0 

,  1  0E19 

) 

ZC0AT(N  •.  21)  =:  CMPLX(C  0 

,  1  0E19 

) 
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ZC0AT(N  .  22)  =  CUPLX(g  e  ,  1  8E!9  ) 

ZC0AT(N  .  23)  =  CMPUf*  «  .  1  6E19  ) 

END  IF 

N  =  N  ♦  (24  «  1  ) 

♦*•»»  end  of  iipedance  coatirj  aodif ications  f-»*»**««»** 


C  syiietry  inputs 

PI  =  ACDS(-1.0) 

SYMTYP  =  'ROT' 
NBLKS  =  30 
IRhSYM  =  I 
CONVERT  =  l.b 

C  Surface  Region  inputs 


RQTLIM  =  PI  /  NBLKS 
XIRG  0 

00  I  I  =  1,  10 
CC(I)  =  0.0 

1  CONTINUE 

DO  2  I  =  1,  3 
TRNS(I)  =  0.0 

2  CONTINUE 

C  .  SPHERE 


XNSEONS  =  6 

CC(  1)  =  1.  »  CONVERT 

CC(  2)  =  CC(1).CC(1) 

TRNS(l)  =  (0.00E-00)  *  CONVERT 

TRNS(2)  =  (3  00E.00)  *  CONVERT 

TRNS(3)  =  (0.00E-00)  a  CONVERT 

XIZAX  =  -.3 

XSUL  =  0 

XSUU  =  3.I415926B 

XSVU  =  -0.10472 

XSVU  =  0.10472 

XNSU  =  24 

XNSV  =  1 

XIORDU  =  4 

XIORDV  =  04 

XIRQ  =  XIRG  -  1 

CAa  LDSURR(XIRG,  XNSEflNS,  CC,  TRNS,  XIZAX, 

1  XSUL,  XSUJ,  XSVL,  XSVU,  XNSU,  XNSV,  XIORDU,  XIORDV) 


NSREG  =  XIRG 

C  CALL  PLOTCHIEF(RUNID,  NBLKS,  JYUTYP,  1,  30.0,  60.0,  0.0) 


C  Define  Frequency  input. 

FRER  -  716  20 

C  Generate  surface  P  and  V  eatrices 


CALL  SUPJi(AT(FREq,  SYMTYP,  NBLKS,  PUATX,  VMATX,  MDSIZE) 
IF  (  JJ  .Eq.  1)  THEN 
CALL  PRTOUT('GEDM’,n 
CALL  PRT0UT('AREA',1) 

END  IF 
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C  Oecoaposa  Uitrices 

CALL  OECOyMCSYMTYP,  NBLKS,  IRHSYM,  PMATX,  VWIATX,  UDSIZE) 
CALL  I0SUB(K0Va-S,  H,  VO.,  *) 

C  Block  t  1 

va(  1}  =  [i.em,  t  eeu) 

VB.(  2)  =  (1  0M0,  0.0M0] 

va(  3)  =  (1.0000  ,  0.0000) 
va(  4)  =  (1.0000  ,  0.0000) 
va(  6)  =  (1.0000,  0.0000) 
va(  6)  ^  h  MU.  0  0000) 
va(  7)  =  (1.0000,  0.0000) 
va(  8)  =  (1.0000  ,  0  0000) 
va(  9)  ^  (1.0000,  0.0000) 
va(  10)  =  (1.0000,  0.0000) 
va(  11)  ^  (1  0000  ,  0.0000) 
va(  12)  =  (1.0000,  0  0000) 

Va(  13)  =  (0  M00,  0.0000) 

va(  14)  =  (0.0000,  0.0000) 
va(  IS)  ^  (0.00M,  0.0000) 

Va(  18)  =  (0.H00  ,  0.0000) 

va(  17)  ^  (0.0000  ,  0.0000) 
va(  18)  =  (0.0000  .  0.0000) 

Va(  19)  =  (0.0000  ,  0.M00) 

va(  20)  =  (0.0000  ,  0.0000) 

va(  21)  ==  (0.0000  ,  0.0000) 

Va(  22)  =  (0.0000  ,  0.0000) 

Va(  23)  =  (0.0000  ,  0.0000) 

Va(  24)  =  (0.0000  ,  0.0000) 

CALL  IQSUB(NDVaS,  1,  Va,  NfOVK) 


C  Oanerata  Surface  Pressures. 


CAa  SURPRS(FREfl,  SYUTYP,  NaxS,  IRHSYM,  PklATX,  VMATX,  MOSIZE) 

CALL  PRTOUT('Va’,l) 

CALL  PRT0UT('SP’,1) 

aOTYP  =  'FAR  ■ 

NUMTHP  =  73 
00  26  I  =  1,  NUMTHP 

THTPHI(1,  I)  =  (I  -  1)  .  6 
THTPHI(2,  I)  =  00 
26  CONTINUE 

C  Calculate  FAR  Field  Matrices 

CAa  aOUAT(FREq,  SYMTYP,  aOTYP,  NaXS,  PMATX,  VMATX,  MOSIZE) 


IFFNRM  =  0 


C  Calculate  Far-Field  Pressures. 

CAa  aOPRS(FRER,  FLOTYP,  NaXS,  IRHSYM,  PMATX,  VMATX,  MOSIZE) 

CAa  PRTOUT (FLOTYP,  1) 

99  CONTINUE 
STOP 
END 


The  output  listed  below  shows  the  input  parameters  for  case  1.  This  is 
followed  by  the  velocity  distribution  and  surface  pressures  for  all  of  the 
subdivisions.  All  of  the  symmetry  strips  are  the  same.  Following  the  surface 
pressures  is  the  normalized  far-field  patterns,  and  next  are  the  surface 
pressures  and  patterns  for  cases  2  and  3. 
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Output  Listing  of  Samule  Run  1 


PROGRAM  r  H  I  L  F  8  3  RUN  MPBI 

l/AIL  19  FEB-9 

REGION  NSE5NS 

1  6 

NSU  NSV  SUL  SUU 

24  1  3.080000F.00  a.3141B9E-01 

SVL 

-  104720E*00 

SVU 

0.104720E'>00 

IQRDRU 

4 

lORORV 

4 

REGION  -  TRANSLATION  (T;,T2J3',,  IZAX 

1  ij.000oeieE>o?  e  0.000000t>00  3 

REGION  -  CCS 

1  0.i000E>0i  0  1000E-01  0  000eE<00  0  0000E<00  0.000eE-00  0.0000E.e0  0.ee00E*0e  e.000eE«08  e.0000E«ee  0.0000E<00 

NO.  OF  SURFACE  AREAS  =24  NO.  OF  SURF.AfE  FiaO  POINTS  =  24 

1 

PROGRAM  C  H  I  E  F  8  8  RUN  URBI  DATE  19-FEB-9 


SURFACE  AREAS  fAR^AS  ARE') 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


0. 17918£-e2 
0.53447E-02 
0.88062E-02 
0.12117E-01 
0.15220E-8i 
0.18063E-01 
0.20597F.-0i 
e.2?7'3E-?: 

0.24671E-01 
0.25942E-0i 
0.26870E-ai 
0.27337E-01 
0.27337E-01 
0,26870E-01 
0.2S942£-0i 
0.24.57r  ■; 

0  227/9u 

0.20S97E-01 

0.18063E-01 

0.1S220E-01 

0.12117E-01 

0.88062E-02 

0.53447E-02 

0.17913E-02 

PROGRAM  C  H  I  E  F  8  8  RUN  MRBI  DATE  19-FEB-9 


FRE(JUENCY=  716.2  RH0=  1000.0  C=  1500  0 


SYMMETRY  TYPE=  ROT  NO  OF  BL0CKS=  30 

VaOCITIES  (BLKf  .  .AkTAI,  REAL ,  IHAG , MAG,  ANG,  REAL  ZCOAT,  IMAG  ZCOAT) 


1 

1 

0  100000EO01 

0  000000E<00 

0  i0000eE->ei 

0  00 

1 

2 

0.100000E>01 

0  e00000E'>00 

0  i0ee00E'>ei 

0  00 

1 

3 

0  100000E<01 

0  ee0000E<00 

0  i00e00E<0i 

0  00 

1 

4 

0  100000E<01 

0  00e000E>00 

0  i00000E<ei 

0  00 

1 

5 

0  I00000E*01 

0  0000e0E<00 

0  I00000E»0i 

0  00 

1 

8 

0  100000E=01 

0  000000E<00 

0  100000E*01 

0.00 

1 

7 

0  !00000E-01 

0  000000E<00 

0  I00000F.*01 

0  00 

1 

;j 

0  100000E-01 

0  000000E<00 

0  100000E-01 

0  00 

1 

9 

0  I00000E<01 

0  00e000E.00 

0  100000E*01 

0  00 

1 

10 

J.102280E.01 

0  000000E<00 

0  100000E*01 

0.00 

1 

11 

0  100000E-01 

0  000000E-0e 

0  100B00E-01 

0  00 

1 

12 

0  100000E-01 

0  00e0e0E>00 

0  i0ee00E>0i 

0  00 

1 

13 

-0  B19151E-00 

0  925700E*0a 

0  108134E.01 

119  28 

1 

14 

0.9693b0E-01 

a  223281E*00 

0  247090E.00 

86  90 

1 

IS 

0  I027SSE-00 

0  «08;40L  01 

0  137204E-00 

41  48 

1 

16 

0  839903!  01 

0  ;.3ii00E-0i 

0  850073E-01 

8  87 

0  MtfefE*ie 
0  e000fiE*n 
0  000000E«00 
0  0000ME«00 
0  000000E«M 
0.0f0000E«00 
0.000000E*N 
0  00M00E*M 
0.0000ME«00 
0  000000E«00 
0.000000E«00 
0.000000E<N 
0.N0000E«00 
0.0H000E«N 
0.000000E«00 
0.000000E>M 


0.000000E-00 
0.0000HE«00 
0.0000HE« 
0.0fll0f0E« 
0.0000HE« 
0.00MI0E« 
0.000000E* 
0.00fffME« 
0.0000ME« 
0.000000E> 
0.0N000E» 

0  0000I0E* 
0.0000NE< 

0  0000ME« 
0.000IME« 

0  0000ME* 
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1 

17 

e.626094E-81 

-0.223428E-01 

0  673414E-01 

-22.93 

0.6ee6eeE«e0 

0.000000E*00 

1 

18 

6.246284E-ei 

-6  334644E-61 

0.411894E-01 

-54.31 

e.000600E«60 

0.000000E«00 

1 

19 

6.2SS898E-02 

-6.304916E-01 

0  365987E-01 

-85.26 

e.666eB6E«e0 

0.00B000E«00 

1 

26 

-6.166377E-61 

-0.286314E-61 

0  232124E-01 

-117.28 

6.eee000E«e0 

0.000000E«60 

1 

21 

-6.168224E-61 

-0  879903E-62 

0.189847E-0I 

-162.39 

6.B000e0E«00 

0.000000EO00 

1 

22 

-B.183030E-01 

0  lB6Se6E~02 

0.I83978E-01 

174.18 

0.B00000E«00 

0.600000E<00 

1 

23 

-0.176939E-01 

e.9618S7E-e2 

0.20e61SE-01 

151.33 

0.000000E«00 

0.000000E*00 

1 

24 

-0.177B78E-01 

0.148317E-B1 

0.230093E-01 

146.51 

0.000000E«00 

0.000000E«00 

FULL  va.ocm  symjctry 

i 

PROGRAM  C  H  I  E  F  6  8  RON  MRBI  DATE  19-FEB-9 

FREQUENCY*  716.2  RH0=  1000.6  C=  1600.0 
SYIKETRY  TYPE*  ROT  NO.  OF  BLOCKS*  30 
SURFACE  PRESSURES  (BLKf,AREA|,  REAL, I MAG, MAG, ANG) 


1 

1 

0  138996E.07 

0.243849E.05 

0.13101BE<07 

1.07 

1 

2 

0.134038E>67 

0  755111E<05 

e.l34249E<07 

3.22 

1 

3 

0.13g889E<07 

0.142a90E«08 

0.140615E>e7 

5.82 

1 

4 

6.147641E-07 

0  23Sa99E<08 

0.148918E«07 

9.11 

1 

5 

0.153652E-e7 

0.354132E>06 

0.15768lE<07 

12.98 

1 

6 

6.15776aE*07 

0.491521E<06 

0.165188E«B7 

17.31 

1 

7 

6.167324E-67 

0  838299E«06 

0.  ie970SE'>07 

22.02 

1 

8 

6.15e995E'>67 

0.771658E<0a 

0.169565E«07 

27.07 

1 

9 

6.137796E*67 

0  875019E-66 

0.163230E>e7 

32.42 

1 

10 

0.117443E*07 

6.916e71E<06 

0. 149105E>07 

38.03 

1 

11 

6.699727E«6e 

0.684173E«0a 

0.124752E*07 

43.85 

1 

12 

6.537282E<06 

0.631044E«0e 

e.828787E*0e 

49.59 

1 

13 

6.666660E<06 

0.006006E<00 

0.000B00E«00 

0.00 

1 

14 

6.660666E«66 

0.660066E«00 

e.006008E«80 

0.00 

1 

15 

6.660066E«N 

0.006H0E«00 

0.000000E«00 

0.00 

1 

16 

6.600006E*06 

0.000660E«00 

0.00600BE«00 

0.00 

1 

17 

6.N6669E*66 

0.006608E<00 

0.000060E«00 

0.00 

1 

18 

6.06CW)3E<00 

6.000600£->00 

0.000000E*00 

0.00 

1 

19 

6.666666E'>06 

0  000606E«06 

0.666000E*00 

0.00 

1 

26 

6.066666E«0P 

0.606600E«60 

0  6000e6E«00 

0.00 

1 

21 

6.66006BE<00 

0.006600E«00 

0.600000E*00 

0.00 

1 

22 

6  066enE<H 

0  M0006E«00 

0.860e60E<00 

6.00 

1 

23 

6.6606e6E«66 

0.060000E«00 

0.066060E<00 

6.00 

1 

24 

6  666660E<66 

0.060006E«00 

0.000000E'*00 

6.00 

»*«  FULL  Va.0CITY  SYUCTRY  •«« 

POfER  0UT(Kf)=  78<6.260 

1 

PROGRAM  C  H  I  E  F  8  8  RUN  MRBI  DATE  19-FEB-9 


FREQUENCY*  718.2  RHO*  1000.0  C=  1600.0 
SYUIETRY  TYPE*  ROT  NO.  OF  BLOCKS*  36 


FAR-FIB.D  PRESSURES  (  PHI, THETA, FFP(RL,IM, MAG, ANG)  ,  NORMALIZED  PAHERN,  SOURCE  LEVa,  01  ) 


6.66 

0.00 

-.266623E'>67 

-.655189E*06 

0.26B949E«07 

-165.34 

6.00  * 

6.H 

B.n 

-.248029E<67 

-.66e300E*06 

e.256665E<07 

-165.18 

-0.08 

6.60 

16. N 

-.246739E«67 

-.659422E«06 

e.249607E*67 

-164.68 

-6.32 

in 

16. N 

-.229266E«67 

-.663957E*66 

6.23862BE«67 

-163.84 

-6.71 

0  60 

26.66 

- . 214270E^07 

-  6e8968E<0a 

0  224470E-67 

-162.86 

-1.24 

0.60 

26. N 

-  196984E-67 

-.673270E'>0a 

0  208163E«07 

-161.13 

-1.96 

6  60 

30. M 

-.178304E-67 

-.676639E*06 

0.196747E-67 

-159.26 

-2.66 

6  00 

35. N 

- .  1B9687E*07 

-.674413E.66 

6  173262E-87 

-157  69 

-3.49 

6  00 

46  H 

-.141487E.67 

-  666a00E«08 

0  166489E'>67 

-154.71 

-4.37 

6  00 

45  N 

-.124797E.67 

-  666963E«66 

0  141033E>07 

-152.24 

-6.28 

6.00 

56. N 

- . 109996E«07 

-.638682E*06 

0  12718eE*67 

-149  88 

-6.18 

6  00 

56  H 

-  973268E<06 

-  612816E<06 

0  116612E>07 

-147.86 

-7.66 

0  00 

66  N 

-  668210E>06 

-  B79294E-06 

0.164373E-07 

-146.29 

-7.89 

0  06 

66.06 

-.783448E«B6 

-  537995E*06 

0  950383E*06 

-146.62 

-8.71 

6  00 

76  H 

-.716364E«06 

-  489182E*0e 

0  867471E>06 

-146.67 

-9  66 
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0.80 

76.00 

-.3e3606E-96 

-.433439E-.08 

0.7S2684E-06 

-146  86 

-10.28 

0  00 

80.00 

-,821697E*06 

-.371842E-06 

0  724310E-06 

-149  13 

-11.07 

0  00 

85  00 

-.586732E.06 

-.304956E-06 

0  881296E*06 

-152  54 

-11.86 

0  00 

90  00 

-.555674E*06 

-  2.34820E-06 

0  893161E*06 

-157  09 

-12.66 

0.00 

95.00 

-  526139E-06 

-  16293 iE.0S 

0  549834E-06 

-162  78 

-13.46 

0.00 

100 . 00 

-.49306lE-0a 

-  912237£-e6 

0  501429E.06 

-169.52 

-14.26 

0.00 

105.00 

-  457477E-.06 

.21B432E*0S 

0  457998E»06 

-177.27 

-15.05 

0.00 

110.00 

-.417085E-08 

0.<29278E.06 

0.41926aC.08 

174.12 

-15.61 

0.00 

116.00 

-.371019E-06 

0.100720E.06 

3.384448E-06 

164  81 

-16.57 

0.00 

120.00 

-  319018E^06 

0  149247E-06 

8  352203E.06 

154.93 

-17.33 

0.00 

125.00 

-.261194E»06 

0 . 188452E-.06 

0.320S16E.08 

144.48 

-18,14 

0.00 

130.00 

-.198128E-06 

0  210756E»08 

0.289248E-06 

133.23 

-19.04 

0.00 

136.00 

-.  130a37E-06 

0  221192E-.06 

0  256991E-06 

120.60 

-20.07 

0.00 

140.08 

-.607a39E*06 

0 . 2i7836E'86 

0  226158E*06 

105.59 

-21.18 

0.00 

145.00 

0.101391E-0S 

0.201777E-06 

0. 202031 E.06 

87  12 

-22.18 

0.00 

150.00 

0.796466E.05 

0.175286E-06 

0.192532E.06 

65.56 

-22.57 

0.00 

165.00 

0.146143E.06 

3.141710E.06 

0  202054E.06 

44  32 

-22.12 

0  00 

160.00 

0.203864E>06 

0. 106242E-08 

0 . 229426E*06 

27  30 

-21.05 

0.00 

166.00 

0.2S3066E-06 

0  704420E<86 

0. 26268 7E<06 

15.55 

-19.88 

0.00 

170.00 

0.290267E.06 

0  417507E.05 

0  29324SE<06 

5.19 

-18.92 

0  00 

176.00 

0  313440E-06 

0.228810E.06 

0.314274E.06 

4  18 

-18.32 

0  00 

180.00 

0  32131'’E-86 

0.163323E-a£ 

0  321731E*06 

2  90 

-18,11 

0.00 

186.00 

0.313440E-86 

0.22e810E-05 

0.314274E.08 

4  18 

-18.32 

0.00 

190.00 

0  290257E.06 

0.417608E.05 

0.293245E«06 

8.19 

-16.92 

0.00 

195.00 

0 . 253066£-0e 

0  704419E-06 

0 . 2828e7E*08 

15.55 

-19.88 

0.00 

200  00 

0.203864E-.06 

0  105242E*06 

0.229427E*06 

27.30 

-21.06 

0.00 

206  00 

0.145143E.06 

0  141716E*96 

0 . 202864E<06 

44.32 

-22.12 

0.00 

210.00 

0.796486E.05 

0.176286E<06 

0.192632E.06 

65.56 

-22.57 

0.00 

216.00 

0.101392E*05 

0  201777E-06 

0  202031E-08 

87.12 

-22.16 

0  00 

220  00 

-.6078J7E*05 

0.217038E.06 

0  226157E-06 

105.59 

-21.18 

0.00 

226.00 

■.130837E-06 

0.221192E-06 

0.2S6991E-06 

120.60 

-20.07 

0  00 

230 . 00 

-.198128E-86 

0.210736E-06 

0.289248E'>06 

133.23 

-19.04 

0.00 

236  00 

-,261194E-08 

0.188452E-06 

8.320915E«06 

144.40 

-18.14 

0.00 

240.00 

-  319018E-06 

0  149247E*06 

0  3S2203E-06 

164.93 

-17.33 

0.00 

246.00 

-.371019E.06 

0.100720E>06 

0  384448E-06 

184.81 

-16.57 

0  00 

260.00 

-  417065E.06 

0  429269E-05 

0.419288E.06 

174.12 

-15.81 

0.00 

266.00 

-,457477E*06 

-.218403E-05 

0  457998E*06 

-177.27 

-15.06 

0.00 

260.00 

-,493061E*06 

-.912238E-05 

0  501429E->08 

-169.52 

-14.28 

0  00 

265.00 

625135E-.06 

- . I62931E*86 

0  649634£<06 

-162.76 

-13.46 

0.00 

270.00 

-.556574E-06 

-  234820E-06 

0.603161E-.06 

-157,09 

-12,86 

0.00 

276.00 

-.S86702E-.06 

-.3049S6E-.06 

0.661295E.06 

-152,54 

-11.86 

0.00 

280 . 00 

-.821697E-06 

-.371842E.06 

0.724310E-06 

-149.13 

-11.07 

0.00 

286.00 

-.663686E-06 

-.433439E.08 

0.792884E..08 

-146.85 

-10.28 

0.00 

290.00 

-.718384E.06 

-.489182E->06 

0.867478E-08 

-145.87 

-9.60 

0.00 

296.00 

- .  783448E->06 

- . 637994E*06 

0.95e383E'>06 

-146.62 

-8,71 

0.00 

300.00 

-.86B210E*06 

-.579294E-08 

8.104373E*07 

-146.29 

-7.89 

0.00 

306  00 

-.97325eE-06 

-  612810E*06 

0.116012E*07 

-147.80 

-7.06 

0.00 

310.00 

-.109996E.07 

-.638502E-06 

0.127188E-07 

-149.86 

-6.18 

0.00 

316  00 

-.124797E*07 

-.656962E-68 

0.141033E-07 

-152.24 

-5  28 

0.00 

320.00 

-.141487E.07 

-.668600E*06 

0.158489E-07 

-164.71 

-4  37 

0.00 

326.00 

-.1B9687E.07 

-.674413E-06 

0  173252E*07 

-167.09 

-3.49 

0.00 

330 . 00 

-.178384E»07 

-.87S639E*06 

0. 190747E-07 

-169.26 

-2  86 

0.00 

336  00 

-  196964E-07 

-.673270E-.06 

0  208153E-07 

-161.13 

-1.90 

0  00 

340.00 

-  214270E*e7 

-.688988E.06 

0.22447e£<e7 

182.68 

-1  24 

0.00 

346.00 

-.229205E-07 

-.663957E.06 

0.23e628E-07 

-163.84 

-0.71 

0.00 

350  00 

-.2<0739E»e7 

.fc59422E-06 

0  249607E*07 

-164.68 

-0.32 

0.00 

366.00 

-.248029E-07 

-  686299E-08 

0 . 256666E«07 

-166.18 

-0.18 

0.00 

360 . 00 

-.250623E»07 

-  656189E-08 

8.2S8949E-07 

-166.34 

0  00 

»  DIRECTION  FOR  PATTERN  NORMALIZATION 

1 

PROGRAM  C  H  I  E  F  8  8  RUN  MRBI  DATE  19-FEB-9 


FHEqUENCY.-  716.2  RHO-^  1808  8  1808  8 

SYMMETRY  TYPE=  ROT  NO  OF  BLOCKS--  30 

VaOCITIES  (BLKf,AREA)[,  REAL ,  IMAG.MAC.ANG,  REAL  ZCOAT,  IMAG  ZCOAT) 


11 


Thompson,  Dubus,  i  Siders 


1 

1 

8.18e8e0E«01 

8.88eeetE<00 

0.ie8000E*01 

8.88 

0.B008ME'>80 

0.808808E«B0 

1 

2 

0  088808E'8e 

0.188808E'>81 

8  08 

8.M0ee8E*68 

8.B0B868E«08 

1 

3 

wieesowE-ai 

B.80e00eE<08 

0  1880e0E..01 

0  00 

B.800B0BE*B8 

B.e0BB8eE«B0 

1 

4 

e.i00e08E*0i 

0.00e008E«08 

0 . 108088E«01 

0.00 

e.e00000E*00 

0.0008e0E«0e 

1 

5 

6.i888eeE-ei 

0  006e00E-80 

0.1000B8E>ei 

0.00 

8.B00000E«00 

0.0000B0E«00 

1 

6 

R.ieeawiE-ei 

0  eee008E<00 

0.i8e800E'>ei 

0.00 

8.800000E«00 

0.000000E«0e 

1 

7 

8.100e08E«81 

0.e8H0eE<08 

0.i0000eE«0i 

0.00 

8.800000E«00 

0.8B0800E«00 

1 

8 

s.ieeeeeE^ei 

8  l00l86E>e0 

0 . 18000eE«81 

8.00 

€.meBBE*00 

0.000000E«00 

1 

9 

8.188008E*91 

e.e8080iE^00 

0.188880E«01 

8  80 

8.80000eE«00 

0.000e00E<^00 

1 

10 

0.ie0e08E«8i 

0.00000eE«e0 

8.ie8008E<81 

0.00 

B.B8B00eE«00 

B.BBBBBBE^Be 

1 

11 

8  ie9808E*81 

0.800088E«00 

0.  leemE-Bi 

0.00 

8.e00000E«00 

0.00e000E«00 

1 

12 

e.i8e8e8E->ei 

0.888808E«88 

0.100080E«01 

0.00 

B.BB00B0E«00 

0.e00BB8E«0e 

1 

13 

-8.288403E«8e 

0.733958E-01 

0  2968S8E«00 

186.63 

B.160000E«B7 

0.000000E«0e 

1 

14 

-8.112223E*§0 

e.i63832E-ee 

e.i98682£«ee 

124.41 

0.160000E«07 

0.000e00E«00 

1 

16 

8.411214E-f2 

0.14E298E'>00 

e.l46346E<00 

88.38 

B.16e0B0E«07 

B.00B000E«00 

1 

16 

S.8683iaE-81 

e.926611E-01 

0.113692E*00 

54.49 

B.1S000BE«07 

B.00000BE«00 

1 

17 

8.e66737E-81 

e.36646eE-01 

0.93e98lE-81 

23.18 

8.160000E«07 

0.000080E«00 

1 

IS 

8.768692E-81 

-8  782738E-82 

0.762417E-01 

-6.67 

0.160000E«07 

0.000000E*00 

1 

19 

8.49e632E-01 

-8  342943E-01 

0 . 698686E-01 

-34.96 

0. 160008E«07 

0.00000eE>00 

1 

28 

8.162133E-81 

-a.448118E'01 

0.4e6663E-01 

-70.93 

8.160000E'>07 

0.000000E^00 

1 

21 

-•.181976E-81 

-f .488886E-81 

0.448883E-01 

-114.03 

0. 15e000E*07 

0.0000BBE*00 

1 

22 

-0.462319E-81 

-•0.3888S4E-81 

0.664886E-01 

-146.43 

8.160000E«07 

B.BBBBBBE-'BB 

1 

23 

-8.e08933E-81 

-8.199423E-01 

0.898364E-01 

-183.21 

0.150000E«07 

0.000000E-'00 

1 

24 

-8.7647HE-81 

-0.141181E-81 

0.777623E-01 

-169.64 

0.160000E407 

0.000000E«00 

«*«  FULL  VaOCITY  SYliCTRY  ««< 

PROGRAM  C  H 

I  E  F  8  8 

RUN  MRBI  DATE  19-FEB-9 

FREOUENCY= 

718.2  RHO=  1000 

.0  C=  1608.8 

SYU4ETRY  TYPE=  ROT  NO  OF 

BLOCKS^  38 

SURFACE  PRESSURES  (BLXf.AREAf, 

REAL, IMAG, MAG, ANG) 

1 

1 

0 . 13''’19E*07 

0.1866e9E'>86 

0.136639E*07 

4.61 

1 

2 

0.137'.33E*87 

0.163607E'>06 

0.138150E<07 

6.38 

1 

3 

8.141384E'>e7 

e.212236E''e6 

e.l4296aE«07 

8.54 

1 

4 

0.1462736*I7 

0.291468E«86 

0.149149E«07 

11.27 

1 

5 

0.16862lE*f7 

0.398219E«86 

0.15S497E«07 

14.63 

1 

6 

0.162676E-87 

8  681918E«86 

8.188620E«07 

18.21 

1 

7 

8  161836E-87 

8.8168S1E«08 

0. 163ie9E*07 

22.18 

1 

8 

8.144821E«07 

0.717e94E«06 

eieieasE^aT 

26.34 

1 

9 

8.133378E*87 

8.7S6604E'>06 

0.1S4839E.07 

30.63 

1 

10 

0.116770E*07 

0.808068E«86 

e.l41E49E«07 

34.42 

A 

11 

0.96e68lE*86 

8.724315E'0e 

e.ll9988E*07 

37.13 

1 

12 

I.7I9718E*88 

9.4e7682E«06 

8  8ei880E«0e 

34.49 

1 

13 

0.429686E«0« 

-8.11H94E«e6 

e.443487E*06 

-14.37 

1 

14 

8.168336E«8e 

-8  246748E«06 

0.297874E*86 

-66.69 

1 

16 

-8.ei6821E«04 

-8  217936E*86 

e.218023E«08 

-91.82 

1 

18 

-8.998477E«86 

-8.136827E-’88 

8  178638E<88 

-126.61 

1 

17 

~0. 12838lE'>80 

-8.649e84E'86 

8.13963SE«88 

-166.82 

1 

18 

-8.1138l4E«8e 

0.112911E'>86 

e.ll4363E«86 

174.33 

1 

19 

-8.73E947E-.05 

0.ei4414E>e6 

0.897989E«85 

146.06 

1 

28 

-8.228199E«86 

8.668168E«85 

8.698494E>85 

189.07 

1 

21 

e.2729«3E<a6 

8.ei2t99E«86 

8.678284Eo85 

66.97 

1 

22 

e.693478E«86 

0  4e8288E«0S 

8.63232eE«B5 

33.67 

1 

23 

0.9914e8E«0& 

0.299134E>06 

8.183665E«06 

16.79 

1 

24 

8.1147e6E«88 

8  211771E*06 

8.116643E«06 

18  46 

Fua  vaociTY  syimetry  »». 


POWER  0UT(KW)=  7883.428 


1 


PROGRAM  C  H  1  E  F  8  8  RUN  MRBI  DATE  19-FEB-9 


FREQUENCY^  716.2  RH0=  18800  C=  1508.0 
SYVMETRY  TYPE=  ROT  NO  OF  BLOCKS^  38 

FAR-FiaO  PRESSURES  (  PHI, THETA. FFP(RL,IM, MAG, ANG)  ,  NORMALIZED  PATTERN,  SOURCE  LEVB.,  DI  ) 


12 


NRL  UEyORANDUM  REPORT  6813 


0  00 

0.00 

-.225230E-07 

-.339591E*06 

0  234135E*87 

-184.15 

0.00  * 

0.00 

5.00 

-.223294E*07 

-.339S14E-06 

0.232271E.07 

-164.02 

-0  07 

0.00 

10.00 

-.217627E-e7 

639204 E*06 

0 . 226820E.07 

-163.63 

-0.28 

0  00 

15.00 

-  208635E’07 

-  e38425E.06 

0  218184E.07 

-162.99 

-0.61 

0.00 

20.00 

-.196939E-07 

-.636783E.06 

0  206978E.07 

-162.08 

-1.07 

0  00 

25.00 

183302E*07 

-  633738E*06 

0.193948E<07 

-160.93 

-1.64 

0.00 

30.00 

-.168B44E*07 

-.828609E.06 

0  179885E>07 

-159.55 

-2  29 

0.00 

35.00 

- .  153480E-07 

-.620606E-06 

0.165534E>07 

-157.98 

-3.01 

0.00 

40.00 

- .  138751E-.07 

-.608855E'>06 

0.151522E-07 

-156.31 

-3.78 

0.00 

45.00 

- .  124974E-0r 

- .  692467E*0e 

0  130306E.07 

-154.64 

-4.57 

0.00 

50.00 

-.112521E-07 

-.570537E«06 

0  128159E-.07 

-153.11 

-5.37 

0.00 

55.00 

101514E.87 

-.542312E-06 

0  115180E.07 

-151  91 

-6.16 

0.00 

60.00 

-.923236E-06 

-  507165E-.06 

6.105337E.07 

-151.22 

-6.94 

0.00 

65.00 

-.845913E»06 

-.484707E-.06 

0.96S154E<08 

-151.22 

-7.70 

0.00 

70.00 

- . 782670E»06 

-.414851E*06 

0 . 885818E«06 

152.07 

-8.44 

0.00 

75.00 

-.731383E»06 

- .  357862E-.06 

0 . 814239E>06 

-153.93 

-9.17 

0.00 

80.00 

-.689582E*06 

- .  294408E<06 

0 . 749800E'>06 

-156.86 

-9.69 

0.00 

85.00 

-.654661E-06 

- .  22S586E*06 

0.692438E>06 

-160.99 

-10.58 

0.00 

90.00 

-.623999E-06 

-.162931E*06 

0  642467E>06 

-166.23 

-11.23 

0.00 

95.00 

-.695031E^06 

- . 783976E*06 

0.600173E'>06 

-172.49 

-11.82 

0.00 

100  00 

-.586256E»08 

-.430427E*04 

0  56S273E.06 

-179  56 

-12.34 

0.00 

105.00 

-.532241E*06 

0.667534E>05 

0.536411E«06 

172.85 

-12.80 

0  00 

110.00 

-.493619E-06 

0 . 132036E'>06 

0  510973EO06 

165.02 

-13.22 

0.00 

115.00 

-.447141E-06 

0.18e833E'>06 

0.485379E<06 

157.11 

-13.67 

0.00 

120.00 

-.390785E-06 

0.234669E'>06 

0  455e32E->06 

149.01 

-14.21 

0.00 

125.00 

-.322930E*06 

0.267532E>06 

0.419360E«06 

140.36 

-14.94 

0.00 

130.00 

-.242650E«06 

0.286109E*06 

0.3761E0E<06 

130.30 

-15.91 

0.00 

135.00 

-.149922E<06 

0.269987E«06 

0.326449E<06 

117.34 

-17.11 

0.00 

140.00 

-.459879E.05 

0  279807E.06 

0 . 283581E*06 

99  33 

-18.34 

0.00 

145.00 

0.664792E<05 

0  257322£<06 

0.266776E«06 

75.51 

-18.90 

0.00 

150.00 

0.183306E>06 

0.22E345E'>06 

0.290486E<06 

50.87 

-16.13 

0.00 

155.00 

0.298966E-06 

0.137667E*06 

0.352933E«06 

32.10 

-16.44 

0.00 

160.00 

0.406947E-06 

0.14&245E>06 

0.433108E«06 

20.02 

-14.66 

0.00 

165.00 

0.500352E«06 

0.111798E.06 

0.51269«E>06 

12.60 

-13.19 

0.00 

170.00 

0.572635E<06 

0.8233S5E«05 

0.578624E-.06 

8.18 

-12.14 

0.00 

175.00 

0.616392E>06 

0.631948E«05 

0.621613E«06 

5  83 

-11.62 

0.00 

180.00 

0.634060E<06 

0.565625E'05 

0  636578E<66 

5.10 

-11.31 

0  00 

185.00 

0.  il  392E-00 

0.631947E*eB 

0  821613E-06 

5.63 

-11.62 

0.00 

190.00 

0.572635E*06 

0.823334E«05 

0.578523E«06 

8.18 

-12.14 

0.00 

195.00 

0.500352E-06 

0.111797E*06 

0.S12690E'>06 

12.60 

-13.19 

0.00 

200.00 

0.406948E-06 

0.14e245E'>06 

0.433108E<06 

20.02 

-14  66 

0.00 

205.00 

0.298966E'>06 

0.ie7566E<06 

0.352933E«06 

32.10 

-16.44 

0.00 

210.00 

0.163307E-06 

0.225345E*06 

0.290485E«0e 

50.87 

-18.13 

0.00 

215.00 

0  664793E*06 

0.2S7321E'>06 

0 . 265770E«06 

75.51 

-16.90 

0.00 

220  00 

-.459875E-.05 

0.279607E<06 

0.283561E<06 

99.33 

-18.34 

0  00 

225  00 

-.149921E*08 

0  289987E>06 

0 . 326449E<0e 

117.34 

-17.11 

0.00 

230.00 

- . 242660E*0e 

0.286109E<06 

0.375150E*06 

130.30 

-15.91 

0.00 

235.00 

-.322938E>06 

0.267532E*06 

0.419360E<06 

140.36 

-14.94 

0  00 

240.00 

-.390785E<06 

0.234e69E'>06 

0.455831E«0e 

149.01 

-14.21 

0.00 

245.00 

-.447141E-06 

0.188833E.06 

0 . 485379E*06 

157.11 

-13.67 

0.00 

250.00 

-.4936igE«06 

0.132036E«06 

0  510973E«06 

165.02 

-13.22 

0.00 

255.00 

-.532241E»06 

0.667533E«05 

0. 53641 1E<06 

172.86 

-12  80 

0.00 

260.00 

- . 565256E«06 

-.430442E'>04 

0.565273E<86 

-179.66 

-12  34 

0.00 

265.00 

-.595031E<06 

-.783976E.05 

0.660173E«06 

-172.49 

-11  82 

0.00 

270.00 

-.823999E.06 

-.162931E-96 

0.842466E-08 

-166.23 

-11  23 

0.00 

275.00 

-.664661E-06 

- . 225686E<06 

0.692438E«06 

-160.99 

-10.58 

0.00 

280.00 

-.689582E«06 

-.2g4408E«06 

0.749806E<0e 

-156.88 

-9.89 

0.00 

285.00 

-.731302E^66 

-.367062E.06 

0.814239E«06 

-153.93 

-9.17 

0.00 

290.00 

-.782670E*06 

-.414850E<06 

0.e858ieE«06 

-162.07 

-8.44 

0.00 

295.00 

-.845913E«06 

- . 464707E.06 

0.9651S3E<06 

-151.22 

-7.70 

0.00 

300.00 

-.923238E-06 

-.607164E-06 

0 . 105337E*07 

-151.22 

-6.94 

0.00 

305.00 

-.101614E«87 

-.642312E-06 

0.115180E<07 

-151.91 

-8.16 

0  00 

310.00 

-  .  U2521E-07 

.570B36E-08 

0.126169E*07 

-163.11 

-5.37 

0.00 

315.00 

-.124974E.07 

-.692467E*00 

0.13e306E«07 

-164.84 

-4.57 

0.00 

320.00 

-.138761E*07 

-.608855E406 

0  151S22E«07 

-158.31 

-3.78 

0.00 

325.00 

-  ■'.53460E*07 

-.620605E«00 

0  165534E>07 

-157.98 

-3.01 

0.00 

330.00 

- . 168544E*87 

-.628609E*06 

0  179886E.07 

-159.55 

-2.29 

0.00 

335.00 

-.183302E-07 

-  633736E^06 

0.193948E*07 

-160.93 

-1.64 

0.00 

340.00 

- . 196939E«07 

- . 636703E*06 

0.296878E>07 

-162.08 

-1.07 

0.00 

345.00 

-.288634E*07 

63842  <£•>  06 

0  218184E«07 

-162.99 

-0  61 

0  00 

350.00 

-,217627E.07 

-  639204E-06 

0.226820E*07 

-163.63 

-0.28 

247.4 


7.7B 


13 


Thompson,  Dubus,  &  Siders 


i.M  35S.M  -  223294E-07  -.839514E*88  «.232271E*«7  -184. #2  -0.07 

0.M  361.  M  -.228230E*07  -.e39B91E<06  0.234135E<e7  -164.1$  0.00 


*  DIRECTION  FOR  PAHERN  NQRKALIZATION 

1 

PROGRAM  C  H  I  E  F  8  8  RUN  MRBI  DATE  19-FEB-9 


FREQUENCY:  716.2  RH0=  1000.0  C=  1500.0 

SYMNETRY  ’YPE=  ROT  NO.  OF  H.0CKS=  30 


vaOCITIES  (BLK|,AREAf,  REAL,  I  MAG.  MAG,  ANG,  REAL  ZCOAT,  IMAG  ZCOAT) 


1 

1 

0.100000E«01 

0.e00000E«M 

1 

2 

0.1B000eE-01 

0  000000E«00 

1 

3 

0.100000E<01 

0.M0e00E*00 

1 

4 

0  100000E«01 

0.000000E*H 

1 

5 

0.10e000E«01 

0.000000E*00 

1 

6 

0.100000E«01 

0  000000E*00 

1 

7 

0.100000E>01 

0.000000E«00 

1 

8 

0.100000E«01 

0  000000E*00 

1 

9 

0.100000E«01 

0.000000E«00 

1 

10 

0.100000E«01 

0.000000E->00 

1 

11 

0.100000E*01 

0.000000E«00 

1 

12 

0.100000E«01 

0.000000E«H 

1 

13 

0.3296B0E-14 

0.566012E-13 

1 

14 

0.243163E-13 

0  363465E-13 

1 

IE 

0.312216E-13 

0.16e352E-13 

1 

16 

0.30068EE-13 

0.17S188E-14 

1 

17 

0.238301E-13 

-0.8SB669E-14 

1 

18 

0.146888E-13 

-0  141628E-13 

1 

19 

0.439377E-14 

-0.1ES414E-13 

1 

20 

-0.669067E-14 

-0.136699E-13 

1 

21 

-0.145778E-13 

-0.979933E-14 

1 

22 

-0.216270E-13 

-0.S34356E-14 

1 

23 

-0.2640  8£- 13 

-0.ISB91SE-14 

1 

24 

-0.269v<)4E-13 

0.694412E-16 

»»» Fua  vaociTY  syiiietry  ♦** 

PROGRAM  C  H 

I  E  F  8  6 

FREUUENCY=  716.2  RH0=  1000 

.0  C::  1600.0 

SYMMETRY  TYPfe 

ROT  NO .  OF 

0LOCKS=  30 

SURFACE  PRESSURES  (BLKf.AREAf, 

REAL,  IMAG, MAG 

1 

1 

0.13S092E*07 

0.137472E<06 

1 

2 

0.136950E«07 

0.101612E«06 

1 

3 

0.140664E«07 

0.234871E'>06 

1 

4 

0.14E111E*07 

0.3063e3E->0e 

1 

5 

0.148962E*07 

0.396436E->06 

1 

6 

0.160868E«07 

0.49e214E-0e 

1 

7 

0.149497E*07 

0.SB8979E«0e 

1 

8 

0.143898E«07 

0.e9012SE«06 

1 

9 

0.133693E<07 

0.762193E«M 

1 

10 

0.118709E«07 

0.76340eE*06 

1 

11 

0.999980E«06 

0.694102E*06 

1 

12 

0.767446E«0e 

0.483704E*0e 

1 

13 

0  E66012E*06 

-0  329060E<06 

1 

14 

0.363465E*06 

-0.243163E«06 

1 

16 

0.1O6362E«06 

-0.312216E<06 

1 

16 

0.17E188E«0E 

-0.30068EE<06 

1 

17 

-0.858669E>0E 

-0.238301E*06 

1 

18 

-0.141O2eE*06 

-0  146880E«0e 

1 

19 

-0.1SE414E<06 

-0  439377E->06 

1 

20 

-0.136699E>0e 

0.569067E«06 

1 

21 

-0.B79933E«06 

0  14S770E<06 

1 

22 

-0  634365E>06 

0.216270E<06 

1 

23 

-0.16S915E<0$ 

0.2e4818E>06 

0.1000e8E<01 

0.00 

0.000080E«00 

0.108080E«01 

0.00 

0.000000E«00 

8.100008E«01 

0.00 

0.000000E«00 

0.100B00E*01 

0.00 

0.060000E«00 

0.100000E-01 

0.00 

0.000000E*00 

0.i00e00E-0i 

0.00 

B.000000E*00 

0.100800E«01 

0.H 

0.000000E«00 

0.100000E*01 

0.00 

0.000000E«00 

0.1M000E«ei 

0.00 

B.090000E«00 

0.100000E«01 

0.00 

0.000000E*00 

0.1B0000E->01 

0.N 

0.000000E^00 

0.1000e0E^01 

0.00 

0.000000E'»00 

0.666968E-13 

86.67 

0.000800E«00 

0.429020E-13 

65.47 

0.000000E«0B 

0.353769E-13 

28.05 

0.000000E«00 

0.30119SE-13 

3.33 

0.000000E*00 

0.2B3197E-13 

-19. 7B 

0.000800E'>00 

0.204f44E-13 

-43.96 

0.000000E«00 

0.161B0SE-13 

-74.21 

0.000000E«00 

0. 147978E-13 

-112.62 

0.0M000E«00 

6.17B8S3E-13 

-146.09 

0.000B06E*00 

0.222774E-13 

-166.12 

0  000000E«00 

0.20S277E-13 

-176.63 

0.000000E*00 

0.29004SE-13 

178.63 

0.00B000E-00 

RUN  MRBI  DATE  19-FEB-9 


.ANG) 


0.135790E«f7 

6.81 

0.138149E>07 

7.55 

0.142611E*07 

9.48 

0.148310E«07 

11.92 

0.1B4141E«07 

14.86 

0.168819E«07 

18.21 

0.161050E*07 

21.83 

0.169S91E«07 

26.62 

0. 163313E«07 

29.38 

0. 141137E*07 

32.74 

0. 121727E*07 

34.77 

0.924144E*06 

31.66 

0.6e69e8E«0e 

-3.33 

0  429028E<06 

-34.63 

0.3B3769E«06 

-61.96 

0.30119BE*0e 

-86.67 

0  263197E«0e 

-109.76 

0.204044E'>06 

-133.96 

0.161S06E«06 

-164.21 

0.147978E«06 

167.38 

0.176e53E«06 

123.91 

0  222774E-06 

103.88 

0  265277E«06 

93.37 

0.B000e6E<00 
0.000000E«00 
0.000000E«00 
0.e00000E«00 
0.000000E«00 
0.000000E«00 
0.00000eE^00 
0.000000E*00 
0.0e000eE*00 
0.090000E«00 
0.00000eE'>00 
0.000000E«00 
0. 100000E'»20 
B.100000E'»20 
0.ie0000E«20 
0.100008E«20 
0.100000E«20 
0.100000E«20 
0.10000BE'>20 
0.100000E*20 
0.100000E«20 
0.100000E*20 
0.1B0B09E^2B 
0.100000E«20 
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1  24  0.S94412E-B4  3 .  JS9S8«F.- 06  0.2902J5L-06  06.83 

M»  FULL  VELOCITY  SYMMETRY  .»* 

POWER  OUT(KW)-  7956  770 


1 


program  C  H  I  t  r  8  3  RU'.V  MPHi  OATE  19-TEfl-9 


FREHUENCY=  716.2  RHO=:  1000.0  C=  1.500  0 
SYkWETRY  TYPE=  ROT  NO,  OF  BLOCKS-  .38 


FAR-FiaO  PRESSURES  (  PM t  THeTA , Frp(Rt  ,  IM,  MAC.ANG)  .  NORMALIZED  PATTERN,  SOURCE  LEVB.,  01  ) 


0.00 

0.00 

-  218313E*g7 

-  300630E.06 

0.22e426F:.07 

-134  62 

0.00 

0.00 

5.00 

-.216S23E*07 

-  601108E.06 

a  224712E.8/ 

-154  48 

-0  07 

0.00 

10.00 

-  211282E*67 

-.602237E.06 

0.219S98E.07 

-164  09 

-0.28 

0.00 

15  00 

-.202958E-07 

-.603609E.08 

0  211743E.07 

-183  44 

-0.58 

0  00 

20  00 

-.192112E*e7 

-.604498E.08 

0  20139SE.07 

-182  53 

-1  02 

0.00 

25  00 

-.179443E-67 

-  803972E-00 

0  189335E.07 

-181  40 

-1.65 

0  00 

30.00 

-  :95704;'e7 

8.17e255£.a7 

-180.07 

-2.18 

0  00 

35.00 

-.  151637E.07 

-  S94323E.08 

0.1628e8E.07 

-158.80 

-2.86 

0.00 

40.00 

-.137904E-07 

-  38297.3F.08 

0.149720E.07 

-157.08 

-3.59 

0.00 

45.00 

-.125048E-07 

-.36S931E.08 

0.137258E97 

15S.65 

-4.36 

0.00 

60.00 

-.113462£*e7 

-.542414E.08 

0.1257eiE.07 

-154.45 

-6.11 

0.00 

55.00 

-.103388E-07 

-.511916E.08 

0.115367E.07 

-153  66 

-5.86 

0.00 

60.00 

-.949194E-08 

-.474217E.06 

e.l08106E.07 

-153.46 

-8.58 

0.00 

65.00 

-.880240E-06 

-.429479E.08 

0  97942eE.08 

-153.99 

-7  28 

0.00 

70.00 

-,a26874'r>06 

-,3?8199E.06 

8.908186F..05 

-155.39 

-7.94 

0.00 

75.00 

-.78337SE-06 

-.321181E*0a 

0.846661E.06 

-157.71 

-8.64 

0.00 

80.00 

-,750675£»86 

-  259583E-08 

0 . 794290E.06 

-180.92 

-9,10 

0  00 

85.00 

-,724546£.06 

-,191771E.86 

0.756268E.06 

-164  95 

-9.59 

0.00 

90.00 

-.701752E-06 

-  128293E.06 

0  713383E.06 

-169.64 

-10,03 

0.00 

95.00 

-.670962E-08 

-.617906E.06 

0.68i767E.06 

-174  80 

-10.43 

0.00 

100.00 

-.852840£*06 

0.309076E.04 

0  652847E.88 

179.73 

-10. 80 

0  00 

105.00 

-.62ei50£‘06 

0  847910E-06 

0  62352SE.08 

174.04 

-11.20 

0  00 

110.00 

-.B''7883E*06 

0, 121919E.08 

0  590585e.86 

166  09 

-11.87 

0.00 

116.00 

-.S'’  .09E-08 

0. 173327E.03 

0.S5126eE.06 

161.67 

-12.27 

0  00 

120  00 

-.454364E*06 

e.21816eE.06 

0 . 504027E.06 

154.36 

-13.06 

0  00 

125.00 

-.36988 IE* 06 

0.256928E.e8 

0 , 449626E.06 

145.31 

-14.04 

0.00 

130  00 

-.260939E*06 

0.28846eE.06 

0.392922E.08 

133.19 

-15  21 

0.00 

136.00 

-  153082E-06 

0.3099eiE.06 

0.346702E.08 

116.28 

-16.32 

0  00 

140  00 

-  246i33t-06 

0.326947E.08 

0 . 327866E.e6 

94.29 

-16.78 

0.00 

146.00 

0.112824E-06 

0.330207E.06 

0.366529E.06 

71.56 

-16.06 

0.00 

150.00 

0.263612E-.08 

0.344729E.06 

0  427909E.08 

53.87 

-14.47 

0  00 

156.00 

0.390911E-06 

e,.347622E.06 

0  6231i8£.06 

41.65 

-12,73 

0.00 

160.00 

0.617617E-88 

0.34e02eE.06 

0  623748E.86 

33.92 

-11.20 

0.00 

165.00 

e.626080E*06 

0.34 7032E.ee 

0.715026E.06 

29.00 

-10.00 

0.00 

170  00 

0.709337E.06 

0,345578E.06 

0 . 789039E.e6 

25.97 

-9.16 

0.00 

176  00 

0.761749E-06 

0.344396E.08 

8.e35986E.06 

24.33 

-8  65 

0.00 

180.00 

0.779648E*06 

0.34396eE.06 

0.8S2144E.06 

23.81 

-8  49 

0.00 

185.00 

0.761749E*06 

0.344397E.06 

0 . 835986E.06 

24.33 

-8  86 

0.00 

190.00 

0.709337E*e8 

0.345578E.e6 

0.789e39E.06 

25,97 

-9  16 

0.00 

196.00 

0  626080E*06 

0  347032E.08 

0.716826E.06 

29.00 

-10  00 

0.00 

200.00 

0,617817E.06 

0.34802eE.08 

0 . 623740E.86 

33.92 

-11  20 

0.00 

205.00 

0  390912E->06 

8.347822E.06 

0.523118E.B8 

41.65 

-12.73 

0.00 

210.00 

0  253612E-06 

e.344729E.e8 

e.427909E.06 

53.67 

-14.47 

0.00 

215.00 

0  11282BE*08 

0.338207E.e6 

0.358S29E  06 

71.65 

-18,06 

0  00 

220.00 

-  246130E»05 

0.326947E.06 

b.32788EE.06 

94.29 

-16.78 

0.00 

225.00 

- ,  163082E-B6 

0.3e9961E.06 

6.345702E.66 

118  28 

-16.32 

0  00 

230  00 

- . 268939E*88 

0  288460E.e8 

0.392922E.06 

133  19 

-15.21 

0.00 

235.00 

-  369681E-06 

0  266928E.e6 

0.449626E.86 

145  31 

-14.04 

0.00 

240,00 

.4543841.06 

0.21B186E.e8 

0.504027E.e6 

154.36 

-13.05 

0.00 

246.00 

- . B23309E.06 

0.173327E.08 

0.65126eE.06 

181.87 

-12.27 

0.00 

250.00 

-.577863E.B6 

0,121919F.08 

e.590585E.06 

188.09 

-11.67 

0  00 

265.00 

-  6261B0E'>06 

0.647910E.05 

0  623526E.e6 

174.04 

-11.20 

0  00 

260  00 

-.6b2840E*06 

0  309e63F..04 

0  652847E.06 

179  73 

-10.80 

0  00 

266  00 

-  678981E.06 

-  617906E.0B 

0,681767E.e6 

-174.80 

-10.43 

0.00 

270.00 

-.7B1762E-06 

-  128292E.e6 

0.713383E.e6 

-169.64 

-10.03 

0  00 

275  00 

-.724646E-06 

-,194770E.06 

0.750268E.08 

-164.95 

-9.59 

0.00 

289.00 

-.75a874r.e5 

-,289S83E.65 

0  794289E.05 

-160.92 

-9.10 
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0.H 

286.  H 

-.783376E*a6 

-.321181E*86 

a.84668aE'>a6 

-167.71 

-8.64 

•  •• 

29f.aa 

-.826674E«86 

-.37819aE«a6 

B.9a8186E->e6 

-166.39 

-7.94 

i  N 

295. aa 

-.88a248E«a6 

-.429478E«ae 

0.97942EE-e6 

-163.99 

-7.28 

a  M 

aaa.H 

-.949194E«a6 

-.474216E*a6 

0  ia8ia8E'>e7 

-163  45 

-6.68 

•  ft 

aas.aa 

- .  ia33S8E*a7 

-.611918E'>aa 

0.116387E*a7 

-153.68 

-6.66 

1  M 

3ia.H 

-.113482E«87 

-.642414E'>ae 

0.126781E^87 

-164.46 

-6.11 

a.M 

316. aa 

-.126848E«a7 

-.66693aE«a6 

a  13726aE«07 

-166.66 

-4.36 

a.H 

32a. aa 

-.137984E*a7 

-.682973E*a6 

a.i4972aE«a7 

-157.08 

-3.59 

a.aa 

325. aa 

-.161837E*a7 

-.694322E>a8 

a.l82888E«e7 

-158.60 

-2.86 

a.aa 

338. aa 

-.1867B4E*87 

-.6aa96BE«a6 

0 . 176286E«a7 

-160.07 

-2.18 

a.aa 

336. aa 

- .  179443E*a7 

-.8a3972E«a8 

0.189334E-07 

-161.40 

-1.66 

a.aa 

34a. aa 

-,192U2E*a7 

-.6a4498E'>a6 

0  281390E«07 

-162.63 

-1.92 

a.aa 

346. aa 

- . ZB2958E«a7 

-.fia36a8E*a6 

0.211743E-07 

-183.44 

-0.68 

a.aa 

35a. aa 

-.2I1282E*a7 

-.8a2237E*a6 

a.219698E'>07 

-164.69 

-0.26 

a.aa 

356. aa 

-,2I8623E-.a7 

-aaiiaaE^es 

0.224712E.07 

-164.48 

-0.07 

a.aa 

3ea.aa 

-.218313E-a7 

-.saaeaaE^ee 

0.220428E*e7 

-164.82 

a.aa 

*  DIRECTION  FOR  PAHERN  NORMALIZATION 


The  agreement  between  the  far-field  directivity  pattern  computed  by  the 
technique  described  by  Butler  [2]  and  the  modified  CHIEF  program  is  so  close 
that  no  difference  would  be  discerned  from  a  plot  of  two  sets  of  data.  Hence, 
the  data  is  listed  in  tabular  form  in  Table  1. 
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Table  1  -  Directivity  Patterns  For  Three  Different  Impedance  Values. 


■ 

Bear i ng 
Angle 

Surf«c« 

Z  =  0.0 

Surface 

II 

o 

Surface 

1  =  RIGID 

Ref.  (2) 

Mod  1 f i «d 

CHIEF 

1 

Ref.  (2) 

Mod i f i ed 

CHIEF 

Ref.  (2) 

Mod i f 1 ed 

CHIEF 

e° 

0.0  dB 

0.0  dB 

0.0  dB 

0.0  dB 

0.0  dB 

0.0  dB 

10° 

-0.4  dB 

-0.3  dB 

-0.3  dS 

-0.3  dB 

-0.3  dB 

-0.3  dB 

20° 

-1.2  dB 

-1.2  dB 

-1.1  dB 

-1.1  dB 

-1.1  dB 

-1.0  dB 

30° 

-2.6  da 

-2.7  dB 

-2.2  dB 

-2.3  dB 

-2.2  dB 

-2.2  dB 

40° 

-4.4  dB 

-4.4  dB 

-3.6  dB 

-3.8  dB 

-3.8  dB 

-3.6  dB 

50° 

-6.1  dB 

-6.2  dB 

-5.2  dB 

-5.4  dB 

-6.1  dB 

-5.1  dB 

60° 

-7.7  dB 

-7.9  dB 

-6.9  dB 

-6.9  dB 

-6.7  dB 

-6.6  dB 

1  70° 

-9.1  dB 

-9.6  dB 

-8.3  dB 

-8.4  dB 

-8.1  dB 

-7.9  dB 

80° 

-11.0  dB 

-11.1  dB 

-9.6  dB 

-9.9  dB 

-8.9  dB 

-9.1  dB 

90° 

-12.8  dB 

-12.7  dB 

-11.0  dB 

-11.2  dB 

-10.0  dB 

-10.0  dB 

100° 

-14.4  dB 

-14.3  dB 

-12.0  dB 

-12.3  dB 

-10.6  dB 

-10.8  dB 

110° 

-16.9  dB 

-16.8  dB 

-13.0  dB 

-13.2  dB 

-11.2  dB 

-11.7  dB 

120° 

-17.8  dB 

-17.3  dB 

-14 .0  dB 

-14.2  dB 

-13.0  dB 

-13.1  dB 

'  130° 

-19.5  da 

-19.0  dB 

-16.6  dB 

-15.9  dB 

-16.7  dB 

-16.2  dB 

140° 

-22.4  dB 

-21.2  dB 

-18.7  dB 

-18.3  dB 

-16.2  dB 

-16.8  dB 

150° 

-22.4  dB 

-22.6  dB 

-17.6  dB 

-18.1  dB 

-14.6  dB 

-14.6  dB 

160° 

-20.8  dB 

-21.1  dB 

-14.7  dB 

-14.7  dB 

-12.1  dB 

-11.2  dB 

170° 

-19. B  dB 

-18.9  dB 

-12.3  dB 

-12.1  dB 

-9.6  dB 

-9.2  dB 

1  180° 

-18.9  dB 

-18.1  dB 

-11.6  dB 

-11.3  dB 

-8.6  dB 

-8.6  dB 

Thompson,  Dubus,  &  Siders 


O 

Fig.  3.  Rapraaentatire  directiritf  patterna  for  a  D0  uniformlx  Tibrating  polar  cap  on 
an  lapedanca-coatad  aphera  for  the  Taluea  of  the  iapadanca  gi^an  in  Table  1. 
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Sample  Ruu  2:  Impedsoice  Coated  Wed£e  With  Spherical  Source 

This  sample  computes  the  acoustic  radiation  from  a  spherical  source  near  a 
wedged-shape  surface  covered  with  an  impedance  material.  The  interesting  part 
of  this  example  is  that  the  geometry  is  generated  using  equation  sets  6 
(spherical  region) ,  1  (rectangular  planar  region) ,  and  10  (finite  element 
inputs) . 


Fig.  8.  Spherickl  kource  near  a  wedged- shaped  baffle. 


Program  Listing  of  Sample  Run  2 


C  Progra*  WEDGETES 
C  WEDGE  WITH  VIBRATING  SPHERE 

C  (COMMON  VARIABLE  DECLARATIONS  GO  HERq 

COMPLEX  PMATX(14884),  VMATX(14884) 

MDSIZE  =  14884 
RUNID  =  'WEDO' 

DATE  =  '23-0Cr-9«’ 

CALL  INITCM 
CALL  OPNSFL 
RHO  -  1000 
C  =  1500 

lipedsnce  coating  aodif Icalpons 
NUMZ  =  9 

N  =  1 
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C  — '  claar  iipedance  for  sphere 

DO  II  I  =  N,  N  -  (8  «  4  )  -  ? 

ZCOAT(I)  =  CMPI.X(0.e,  10) 

10  CONTINUE 

N  =  N  .  (8  «  4  ) 

C  -  clear  iapedance  for  SIDEl 

DO  ?0  I  =  N,  N  .  (3  «  8  )  -  1 

ZCOAT(I)  =  "MPLX(0.0.  00) 

20  CONTINUE 

N  =  N  .  (3  .  8  ) 

C  -  claar  iapedance  for  SI0E2 

DO  30  I  =  N,  N  .  (6  *  3  )  -  1 

ZCQAT(I)  =  CMPLX(0  0,  00) 

30  '■TNTINUE 

N  =  N  «  (8  «  3  ) 

C  -  clear  irpedance  for  BOTTQU 

00  40  I  =  N,  N  .  (3  *  3  )  -  I 

;C0AT(I)  =  CUPLX(0  0,  00) 

40  CONTINUE 

N  =  N  .  (3  *  3  ) 

C  -  claar  iapedance  for  EE  SIDEl 

00  60  I  =  N,  N  *  (1  •  18  )  1 

ZCOAT(I)  =  cun  X(0  0,  0  0) 

58  CONTINUE 

N  ::  N  •  (1  «  18  ) 

C  -  clear  iapedance  for  EE  SIDE2 

DO  60  I  =  N,  N  ♦  (1  •  18  )  -  1 

IvOAT(I)  =  CNPU(0.l,  0.1) 

80  CONTINUE 

N  =  N  .  (1  *  18  ) 

C  -  load  iapedance  for  SENSOR 

00  70  I  =  N,  N  ♦  (3  *  3  )  -  1 

ZCOAT(I)  =  CMPLX(.3o6  ,  00  ) 

70  CONTINUE 

N  =  N  .  (3  .  3  ) 


Ctaaaataaaa*  end  of  iapedance  coating  aodif ications  ••«*(«•«•••• 

OPEN (UNI r=NUNPHT, EiLE=aUNID// OUT ' .STATUS^ ’NE»’ , 

1  E0K»t='P0RMAnED') 


C  s/aaetr/  inputs 

PI  =  ACQS(-1  0) 

SYkfTYP  =  'RF' 
NBLKS  =  2  «•  2 

IRHSYM  =  1 
CONVERT  =10 

C  Surface  Region  inputs 

ROTLIM  =  PI  /  NBLKS 
XIRG  =  0 

00  1  I  =  1,  10 
CC(I)  =  0  0 

1  CONTINUE 

00  2  I  =  1,  3 
TRNS(I)  =  0.0 

2  CONTINUE 

C  -  - sphere 
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XSStQNi 

.  .  q 

cr(  !) 

-  009S2S 

«  CONVERT 

;c(  2) 

■  CCCnytC'i 

> 

’  ■  ■  •  ^ 

»  lONVffiT 

rRNS(Z) 

^  ,;0  00e.00) 

*  CONVERT 

I. ‘■'■IS  (3) 

^  (0.00E.00) 

»  CONVERT 

XIZAX 

^  .3 

XSIX. 

=  0 

XGiAJ 

=  3  1A1.S926 

XSVL 

^  0  000000 

X3'*U 

=  1  a. '0795 

-  9 

■•• 

xinROu 

-  4 

XICSOV 

-  4 

XiRC  ; 

XIR-L  .  1 

CAl_  LD 

SjRR(  IRQ,  XNSEI^NS,  CC,  TRNS 

1 

XS‘X  XSUU,  XSVL.  XSVU 

- SI 

DEI 

cc;  ') 

0508 

•  CONVERT 

1KN3(1) 

-  ;0  00E.00) 

♦  CONVERT 

|RNS;2) 

=  (0  00E*00) 

*  CONVERT 

TRNS(3) 

=  {0286.00) 

•  CONVERT 

XTZAX 

=  .1 

xsm. 

^  0 

*  CONVERT 

XSUJ 

^  0508 

♦  CONVERT 

A,  ■■ 

-  ?’ 

•  CONVERT 

xs  zu 

-  i?zib2& 

•  CONVERT 

XRSU 

=  3 

XNSy 

=  6 

XIORDU 

-  4 

XIOROV 

^  i 

xiRG  -  ; 

XIRG  •  1 

CALL  LDSURR(XIRG,  XNSEHNS,  fC,  TRNS, 

1 

XSiA_,  XSUU,  XSVI.,  XSVU, 

XASEQNS 

=  1 

etc  1) 

-  0508 

«  CONVERT 

TRRSCI) 

=  (0  00E.00) 

•  CONVERT 

TRNS (2) 

=  (0  00E.00) 

*  CONVERT 

TRNS (3) 

=  (.0280.00) 

*  CONVERT 

XIZAX 

=  -2 

XT'S. 

=  0  < 

‘  CONVERT 

XovAj 

=  123825  « 

1  CONVERT 

XS'Zl. 

=  0 

*  CONVERT 

aSVU 

-  0508 

*  C-'.-IVdlT 

XNSU 

=  8 

XNS'y 

-  3 

XIORDU 

-  4 

XIORDV 

-  4 

X'HG  =  XIRC  .  1 

CA:  .  DO 

■jrr;x,rg,  xnsehns,  cc.  trns 

] 

XSIX,  XSUU,  XSVL.  XSVU 

-  aairoM 

XN'L.rN3 

--  1 

CC'  I) 

=-  25 '025 

»  CONVERT 

IRNSfl) 

^  (0  00F.80) 

•  CONVERT 

IRN0;2, 

^  (8  00f,.00) 

*  CONVERT 

r '  . 

=  f  .00, 

<  CONVERT 

X ;  /AX 

-  *3 

X'L* 

-  0  t 

CONVERT 

XT,;. 

026A  « 

CONVERT 

X'  l; 

.  0 

*  CONVERT 

X'iVj 

0508 

.  CONVERT 

XM,' , 

-  3 

XN:.V 

-  3 

yi/:  L -A- 

4 

,  XIZAX, 

XNSU,  XNSV,  XIORDU,  XIORDV) 


XIZAX, 
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c 


XIOROV  =  4 
XIRG  =  XIRG  »  1 

CALL  LOSURRtXIRG,  XNSEHNS,  CC,  TRNS,  XIZAX, 

1  XSUL,  XSUU,  XSVL,  XSVU,  XNSO,  XNSV, 

-  FE  SIDEl  - 


XrORDU,  XIOROV) 


XNSERNS  -10 
CC(  1)  =  0 

CC(  2)  =  0 
TRNS(l)  =  (0.00E«00) 
TRMS(2)  =  (0.00E*00) 
TRNS (3)  =  (0.00E*00) 
XIZAX  =  *3 
XSUL  =  -1 
XSUU  =  1 
XSVL  =  -1 
XSVU  =  1 
XNSU  =  1 

XNSV  =  1 

XIOROU  =  4 
XIOROV  -  4 


«  CONVERT 
»  CONVERT 
»  CONVERT 
«  CONVERT 
»  CONVERT 


C  R<a(f  in  coordinita  and  aleaant  data. 


ICOOR  =>  0 

0PEN(UNIT=8,FILE=’  lOMAQC.WED’ ,STATUS=’OLD’ .F0R«=’F0RyAnED’) 

1401  READ(8.1402,ERR=1403)  IDUM,  (COORDS(IDUI, J),J=1.3) 

1402  FORyAT(BNn0,3(BNE10.0)  ) 

IF  (lOUU  GT.  ICOOR)  ICOOR  =  IDlAil 
IF  (IDUM  .GT.  MAXCOR)  THEN 

WRITE  (8,*)  'Error:  Coordinata  indax  in  coordinata  ‘  // 

.  'data  ia  greater  than  MAXCOR. ’ 

STOP 
ENO  IF 

GOTO  1401 

1403  aOSE(B) 

C  Diaplay  coordinata  data. 

IRITE(NUNPRT,1407) 

1407  F0RMAT(1H1,'  COORDINATE  DATA',/) 

00  1409  1=1,  ICOOR 

fRITE(NUNPRT,1408)  I,  (COORDS(I, J), J=l,3) 

1408  F0RyAT(I6,6X,3E16.4) 

1409  CONTINUE 

IQ.QI  c  0 

OP9t(UNIT=0  FILE?’  lOHJB  IOr,STATUS=’OLO’,FORJil=’FORIlATTED’) 

1404  READ(8,1406,LRR=1406)  IDUH.  (ia.TS(J,ia.EM*l),J=l,8) 

1406  FORMAT (T7, 9 {BNI6)  ) 

IB3I  =  laai  <  1 
GOTO  1404 
1408  aOSE(8) 

C  Display  aleaant  data. 

fRITE(NUNPRT,1410) 

1410  FORMAT///,’  aaiEHT  DATA’,/) 

DO  1412 1  =  1,  laa 

tRITE(NUNPRT,1411)  I,  (iaTS(J,I),J=l,8) 

1411  FORMAT(I6,10X,8I5) 

1412  CONTINUE 

00  1414  I  =  1,  lELBt 
DO  1413  J  =  1,  8 
CC(J)  =  iaTS(J,  I) 

1413  CONTINUE 

XIRG  =  XIRG  ♦  1 

CAa  LOSURR(XIRG,  10,  CC,  TRNS,  XIZAX, 
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1  XSUL,  XSUU,  XSVL,  XSVU,  XNSU,  XNSV,  XIORDU,  XIORDV) 

U14  CONTINUE 

C - F£  <IIDE2  - 


XNSERNS  =10 

CC(  1)  =  0  ♦  CONVERT 
CC(  2)  =  0  •  CONVERT 
TRNS(l)  =  (0,00E.00)  ♦  CONVERT 
TRNS(2)  =  (0.00E.00)  »  CONVERT 
TRNS(3i  =  (0.00E-00)  ♦  CONVERT 


XIZAX 

=  -3 

XS'JL 

=  -1 

XSUU 

=  1 

XSVL 

=  -1 

XSVU 

=  1 

XNSU 

=  1 

XNSV 

=  1 

XIORDU 

=  4 

XIORDV 

=  4 

lELEM  =  0 

0 PEN  (UN I  T=  9 ,  F I1.E=  ’  I OB-lilS .  102  ’ ,  STATUS= '  OLD ’ ,  FOR«=  ’  FORUAHED  ’ ) 
1504  READ(8,1605,ERR=1506)  IDUM,  (ia.TS(J,IB.EM.l),J=l,8) 

1508  FORMAT  (;T7,9(BNI5)  ) 

ia.ai  =  lEiai  ♦  i 

GOTO  1504 
1506  CL0SE(8) 

C  Display  elesent  data. 

fRITE(NUNPRT,lS10) 

1510  FORMAIC//,  ‘  aEMENT  DATA’,/) 

00  1512  I  =  1,  laa 

IRITE(NUNPRT,1S11)  I,  (IB-TS(J,I),J=1,8) 

1511  FORMAT(I6,10X,8IS) 

1612  CONTINUE 

DO  1514  I  =  1,  laEM 


00  1513  J  =  1,  8 

CC(J)  =  IELTS(J,  I) 

1513  CONTINUE 

XIRG  =  XIRG  ♦  1 

CALL  LOSURR(XIRG,  10,  CC,  TRNS,  XIZAX, 

1  XSUL,  XSUU,  XSVL,  XSVU,  XNSU,  XNSV,  XIORDU,  XIORDV) 

1514  CONTINUE 

C  -  SENSOR  - 


XNSEHNS  =  1 


CC(  1) 

= 

.0286  *  CONVERT 

TRNS(l) 

- 

(0  00E*00)  *  CONVERT 

TRNS (2) 

r 

(0  00F.00)  •  CONVERT 

TRNS(3) 

= 

(0  00E-00)  «  CONVERT 

XIZAX 

-3 

XSUL 

= 

0  •  CONVERT 

XSUU 

0508  •  CONVERT 

XSVL 

= 

0  *  CONVERT 

XSVU 

0508  »  CONVERT 

XNSU 

3 

XNSV 

- 

7 

XlOPljU 

_ 

4 

XIQPDV 

= 

4 

XIRO  =  XIRG  •  1 

CALI.  LDtURR(XlRG,  XNSEONS,  CC,  trwc_  XIZAX, 

1  XSUL,  XSUU,  XSVL,  XSVU,  XNSU,  XNSV,  XIORDU,  XIORDV) 

NSRPC,  =  XIRG 
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C  CAU  PLOTCHieF(RUNID,  NBLKS,  SYIOTYP,  1,  30.0,  60.0,  0,0) 
CALL  PRT0Lrr('GE0M',l) 

C  Define  Frequency  input. 

DO  99  FREI)  ^  2009  ,  3000  ,  1000 

C  Generate  surface  P  and  V  eatrices. 


CALL  SURMAT(FREq,  SYUTYP,  K0LKS,  PUATX,  VMATX,  MDSIZE) 
CALL  PRT0UT(‘AREA',1) 


C  Decoapose  Matrices 


CALL  DKQkWCSYMTYP,  NBLKS,  IRHSYU,  PMATX,  VMATX,  MDSIZE) 
CALL  ICSUB(NDVa^,  10,  VH.,  0) 

C  Block  I  1 

Va.(  1)  ^  (1.0000,  0.0000) 

VB.(  2)  =  (1.0000,  0.0000) 
va(  3)  =  (1.0000  ,  0.0000) 
va.(  4)  =  U  MM,  0.0000) 

VB.(  E)  =  (1.0000,  0.0000) 

va(  6)  =  (1.00H,  0.0000) 

VB.(  7)  =  (1.0000,  0.0000) 

va(  8)  =  (1.0000  ,  0.0000) 
va(  9)  =  (1.0000  ,  0.0000) 
va(  10)  =  (1.0000  ,  0.0000) 

VH.(  11)  =  (1.0000,  0.0000) 

va(  12)  =  (1.0000  ,  0.0000) 

Va(  13)  =  (1  8000  ,  0.0000) 

va(  14)  =  (1.0000,  0.0000) 

Va(  16)  ::  (1.0000  ,  0.0000) 

VB.(  16)  =  (1.0000,  0.0000) 

VB.(  17)  =  (1  0000,  0.0000) 

va(  1.6)  1.0000  ,  0.0000) 

va.(  19)  (1.00N,  0.0000) 

VEL(  20)  >  (1.0000,  e.M00) 

va.(  21)  =  (1.0000,  0.0000) 
va(  22)  ==  (1.0000  ,  0.0000) 
va(  23)  =  (1.00M,  0.00M) 

VB.(  24)  =  (1.0000,  0  0000) 

va.(  26)  =  (1.B0M,  0.0000) 

VB.(  26)  =  (1.0009,  0.0000) 

VE.(  27)  =  (1.09N,  0.0000) 

VB.(  28)  =  (1.00N,  0.0H0) 

Va(  29)  =  (1  00N,  0.0000) 
va.(  30)  =  (1.00M,  0.00M) 
va(  31)  =  (1.0000  ,  0  0000) 
va(  32]  =  (1.0000  ,  0.0000) 
va.(  33)  ^  (0  00M,  0.0000y 
va(  34)  =  (9.9000  ,  0.0000] 
va(  36]  =  (0  0000,  0  00N) 
va.(  36)  =  (0.00M,  0.0000] 
va(  37]  =  (0.N00  ,  0.0000) 
va(  38)  =  (0.0000  ,  0.0000] 
va(  39)  =  (0.0000,  0.0000] 

va(  40}  =  (0  0000  ,  0.0000) 

va(  41)  =  (0  N00  ,  0.0000) 
va(  42)  =  (0.00M,  0.0000) 
va(  43)  =  (0.0000  ,  0.0000) 
va(  44)  =  (0.0000,  0.0000) 
va(  46)  =  (0.M00  ,  0.0000) 
va(  46)  =  (0.0000  ,  0.0000) 
va(  47)  ^  (0.00N,  0.0000) 
va(  48}  =  (0.0000  ,  0.0000) 
va(  49)  3  (0.00N,  0.0000) 
va(  60)  =  (0.0000  ,  0.0000) 
ya(  61)  >  (0  0000  ,  0.0000) 
va(  62)  =  (0  0000,  0  0000) 
va(  63)  r  (0  0000,  0.0000) 
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va(  s.) 

d  ;H?3) 

VEL(  Sb) 

- 

-J;'"  J  .  .  SOt'iJ) 

VEL(,  S6) 

- 

-j 

vac  so 

■) 

va{  5m'! 

•2V.  1  )?v3-3) 

vac  69) 

O.:?60a) 

vac  83) 

= 

V£L(  ol) 

= 

vac  82) 

';3.0:\'3,  K.0€30) 

vac  83-1 

- 

>'  •»  4,^,3  1-*  nn.tn\ 

^  .  1  V  ,  V  .  } 

vac  80 

3.0 ’=8,  3.3S33) 

vac  85) 

Va'  58'i 

- 

1  7  ^ 

V5LC 

•ii 

va; 

Ci  -  1,  ■i.k‘3CiJ) 

va;  ?=;) 

- 

5  -  '  ^ 

va;  7i) 

:  i  :  '0.  V03e0) 

vac 

,  i  0  .  .5''00) 

vac 

,  ,  >5  3303) 

vac  ’J) 

- 

3000) 

vac  'O 

- 

,3  ,  >).:,'0g3) 

va  ■  , 

= 

i, 

vac  '5) 

:r 

IV'  '’.t/-,100) 

vac  ■  i'l 

.^.5033) 

vac  n) 

- 

0, 

vac 

- 

0.j0eg) 

Vo.;  80) 

- 

C0.„3;nl,  4,3030) 

VB-C  81) 

= 

;0  '  •rr'i),  0.3«30) 

vac  82) 

- 

..0.  ';i30,  -.'i.kees) 

vac  '-v  n 

i 

•■'.?'333) 

va.  6.1, 

3,3o'30) 

vac  85) 

= 

C0.0»J0,  3.0000) 

vac  88) 

(0,.i0,!3,  '<.0000) 

vac  80 

0,0000) 

vac  88, 

= 

Ce.-’sjg,  3.0000) 

vac  89: 

= 

,;0.?0J0,  0.?330) 

vac  90) 

= 

(0.00i0‘  0.0603) 

vac  91; 

= 

!'3.0003,  a  0000) 

Vfu,  6"; 

■.;••,  •  t‘m0«i) 

VG-C  93} 

= 

C*}. iliji'i*,  3.3003) 

vac 

- 

,t’.’W00,  3,0000) 

vac  98) 

= 

I3.ii3r0,  3,0000) 

vac  96) 

= 

!0.!Vf0,  0.0030) 

vac  921 

= 

(3.3«O0,  0.0000) 

vac  98) 

- 

(3.3303,  3.0000) 

VELC  99) 

= 

; 0.3030,  0,0380) 

vac  130) 

= 

(ii.'.’prr,  i'.ri000} 

va  c  101  j 

- 

ie.ib.'fii,  lJ,bW00) 

vaci02) 

- 

(8.0808,  0.0080) 

va.ci03) 

(0,8038,  6,0080) 

vaci04) 

n 

(0.80613  3  0000) 

varies) 

= 

C0.80V-P,  6.0060) 

vac  lob) 

-- 

(0.Ci'36,  0,0006) 

vacie?) 

= 

(0.8003,  e,e000j 

vaci08) 

- 

(0.8r0.',  0,0030) 

ve.citii; 

= 

(''  .(■■  ,  8  .  r’Srr  ) 

vaciie) 

= 

(3.0('8£,  0.0000) 

vaciii) 

= 

(3.0i'80,  0.6000) 

varn?) 

- 

C0,8'-3-,  6, 0030) 

vaciisi 

- 

(3  f  : ' ,  f .8000) 

vfL(;;<) 

- 

(i  .8  '  ,  ,  3  6'' 33) 

vacns) 

= 

(,  .»''i?,  8  3;:30) 

VElCIIB'; 

(3,8  ,  8, 0303) 

VGOO^' 

C.'  .r'---,  !  8>>53; 

VB.CIIS, 

--- 

(0  it-.';,  6.60i'/0) 

VB.CllSj 

(i  f 6,0300) 

va.riie) 

(t-  6,f'30e) 

va.cij!  j 

(k  f83  ^  6  0300) 

vacao 

fk  ,  6  003w) 

C*li  iCSUfcl(f':  V!:  ,1,  '  ,  NWDVLC) 
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C  Generit*  Surface  Pressures. 

CALL  SL«PRS(FREH,  SYlfTYP,  NBLKS,  IRHSYM,  PMATX,  VMATX,  MOSIZE) 
CAU  PRT0UT(’VEL’,1) 

CAL'  PRT0i"'('SP’.n 


FLOTYP  =  ’FAR  ’ 

NUfTHP  =  19 

00  26  I  =  1,  NUMTHP 

IHTPHICl,  I)  =  (I  -  1)  »  10 
THTPHI(2,  I)  =  0.6 
25  CONTINUE 


C  Calculate  FAR  Field  Matrices 

CALL  FLDMAT(FREQ,  SYMTYP,  aOTYP,  NBLKS,  PMATX,  VMATX,  UDSIZE) 

IFFNRM  =  8 


C 


Calculate  Fir-Field  Pressures. 

CALL  FLOPRS(FRa|.  FLOTYP,  NBLKS,  IRHSYM,  PMATX.  VMATX,  MDSIZE) 


CALL  PRTOUT (FLOTYP,  8) 
99  CONTINUE 


STOP 


Output  Listing  of  Sample  Run  2 


1  COORDINATE  DATA 


1 

8.3636'^  48 

2 

8.482Vi:  08 

3 

8. 68 286-88 

4 

8  88126-88 

6 

8.78856-88 

a 

8.85976-88 

7 

8  95986-88 

a 

8.18936-81 

9 

8.15526-81 

18 

8.14116-81 

11 

8.12786-81 

12 

8.11296-81 

13 

8.98786-82 

14 

0.84676-82 

16 

8.33876-81 

16 

8.31846-81 

17 

8.28226-81 

18 

8.25486-81 

19 

8.22586-01 

28 

8.19766-81 

21 

8.16936-81 

22 

8.58606-81 

23 

8.46576-81 

24 

8.42336-81 

25 

8.36186-81 

28 

8  33876-01 

27 

8  29636-01 

28 

8.25466-81 

29 

8  58686-81 

30 

8.46576-81 

31 

8.42336-01 

32 

8.38186-81 

33 

0  33876-81 

34 

8.29636-81 

35 

8  25486-81 

38 

8  58806-01 

37 

8.46576-61 

8.50886-01 

0.1524E«B8 

8.58886-81 

0.1694E<8e 

8.58686-01 

8.1863E'00 

8.50886-01 

0.2032E*M 

0.50606-81 

6  22026*00 

0.50686-01 

8  23716*88 

0.50886-01 

8  25486*08 

8.50886-01 

0  1524E*00 

8.58886-01 

8  1694E*H 

0.50606-01 

8  18836*08 

0.50686-01 

8  20326*00 

0.50886-01 

0  22626*88 

8.58886-01 

8.237lE*80 

8.58686-01 

0.2548E*M 

8 . 58666-01 

8.15246*86 

8.50806-81 

8  16946*80 

0.58886-81 

0. 18636*80 

8 . 58806-81 

8.20326*00 

0.60886-01 

0.22026*00 

0  50886-01 

8  23716*80 

0  50886-01 

8.25406*00 

0.50886-81 

8.15246*80 

8  50886-81 

8  16946*88 

0  50806-01 

8.16836*08 

0.58886-01 

0  28326*08 

0.50886-61 

0.22026*80 

0.50806-01 

0.23716*00 

0  50806-81 

0.25406*80 

0.33876-01 

0  15246*08 

0.33876-01 

0.16946*00 

0  33876-01 

0.18636*00 

8.33876-81 

8  28326*80 

8.33876-01 

0  22026*80 

0  33876-81 

8  23716-86 

8  33876-01 

0.25486*00 

6  16936-01 

0.15246*80 

0  16936-01 

0.16946*08 

28 
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38 

0, 

42^3E^0i. 

e 

-fc 

0. 

1b'' 38-06 

39 

0 

3Oi0E-0i 

e 

'•f 

•0' 

0 

7032E-@0 

40 

0. 

3387E-0 

e. 

.'5C 

0' 

0 

41 

0 

29638-0; 

q 

1C' 

0 

?3’1E.0P 

42 

0. 

2fi8E-?' 

& 

1  /  ‘ 

' 

•  »■- 

0 

25i0E>0e 

43 

0 

5ee0E-p- 

0 

li 

A-_ 

0 

15;'4E>00 

44 

0. 

46S7E-0: 

0 

1^  - 

0. 

0 

1594E.e0 

45 

0 

4233E-e’ 

0 

i  ^ 

a  . 

3 

136 38 *08 

46 

0 

38i0E-6; 

0 

2C 

4  C 

-06 

0. 

20:-2£.0e 

47 

B. 

3387E-S: 

0 

1  ■ 

-.L 

0f. 

0 

22t'2£-00 

48 

V 

4>. 

ir- 

~rJo 

0 

23'1E»0I' 

49 

0 

.2540E--0; 

4  7 

•0' 

0 

25-i0E.0E 

ELEMENT  DATA 


1 

1 

V- 

0 

U 

0 

2 

;! 

H 

10 

> 

3 

C 

0 

4 

4 

i) 

1  !•; 

i’ 

5 

6 

? 

7 

n 

9 

A 

0 

t; 

3 

8 

9 

i  ? 

0 

C 

0 

9 

la 

1 1 

1 3 

0 

0 

10 

L 

‘,2 

1  V  •  ; 

, 

0 

fc 

0 

11 

13 

5e 

0 

0 

0 

12 

[3 

^  •> 

o  .V 

;• 

0 

0 

13 

' ; 

l-l 

0 

3 

1* 

>;■ 

15 

3 

s 

»? 

0 

18 

18 

19 

9v‘ 

3 

0 

0 

17 

1  c> 

2  3 

/7  ■'*  - 

0 

*4 

0 

18 

?0 

21 

?B 

0 

0 

19 

29 

23 

30 

■:) 

0 

0 

20 

23 

24 

31  Vi 

0 

,.i 

0 

21 

2 » 

2S 

32  r 

9 

0 

0 

22 

9" 

■'  ' 

A 

) 

t* 

23 

^  ■ 

.1- 

'J 

A 

0 

24 

V 

2B 

3S 

0 

3 

g 

25 

29 

30 

3/  :i-> 

3 

0 

H 

0 

26 

30 

31 

33  i 

■3 

0 

0 

27 

3: 

32 

39  hi 

0 

3 

0 

28 

32 

33 

40  V' 

0 

3 

0 

29 

33 

34 

i  '  ;  1 

1 

« 

0 

3 

30 

34 

3'> 

<? 

i 

31 

.3  • 

;-:i 

32 

sr 

33 

4  3 

3 

'4 

d 

0 

33 

33 

39 

45 

*; 

gi 

0 

0 

34 

39 

40 

4  7 

0 

35 

40 

41 

<>' 

0 

0 

36 

41 

yi 

t  \i  1  • 

.4 

0 

uROGi-.Aiil 

f'  I  '  • 

t>  o 

R'JN  WEDG 

DATE  23-0CT-9 


EGION 

NSEQNS 

nSU 

NSV 

lU 

suu 

SVL 

svu 

lORDRU 

IQRDRV 

1 

6 

3 

4 

0 

'•860!?  -00 

0 

;i;4i59E.0i 

0  0000008.00 

0.167080E.01 

4 

4 

2 

1 

3 

C 

0  ■ 

’  ‘(!".80  -00 

0 

5.'38007-01 

0.000000E.00 

0  12382SE.00 

4 

4 

3 

1 

0 

.'i 

0  ■ 

■0ti0?  .00 

3, 

I  ■3825E.00 

0.000000E.00 

0  508000E-01 

4 

4 

4 

1 

3 

3 

0  i 

1  ’0008  .00 

0 

?■  4000»--01 

0  eeeeeeE.ee 

0  SBeeaeE-Bi 

4 

4 

5 

10 

1 

L 

.'0000  .01 

0 

100800E.01 

-.ie0000E.0i 

0.ie000eE.0i 

4 

4 

c 

10 

- 

'■0303'  .01 

0 

1,  0003F.01 

..100000E.B1 

e.i00000E.0i 

4 

4 

1 

10 

-01 

1  '0000  -01 

- . i000e0E.0i 

0.100000E.01 

4 

4 

6 

18 

i 

] 

•.‘0«02:,.0l 

0 

130000E.01 

-.100000E.01 

0.i00ee0E.0i 

4 

4 

9 

10 

1 

1 

'  771607:  .01 

0 

1?0006-:.01 

-  ie00e0E.ei 

0.100000E.01 

4 

4 

10 

10 

i 

] 

■'3000..  01 

S 

lP0000h.01 

..i00000E.ei 

0.100000E.01 

4 

4 

11 

10 

1 

- 

■-0003  ,.01 

0 

1000001 .01 

-.100000E.B1 

0.100000E.01 

4 

4 

12 

10 

1 

- 

•.•0000  -01 

0 

100000F.01 

-.100000E.01 

0.100000E.01 

4 

4 

13 

10 

i 

? 

■  '0ft0»  -01 

0 

130000E.01 

.  100000E.01 

0.100000E.01 

4 

4 

;  A 

10 

-3: 

0 

;->0000S'.0! 

.100000E.01 

e.i00000E.0i 

4 

4 

1  5 

10 

- 

-■■  .0. 

7 

■.70030  -01 

-  i0e000E.0i 

0.100000E.01 

4 

4 

/ 


Thompson,  Dubus,  &  Siders 


16 

10 

1 

1 

- .  100000E'>01 

0.i0000eE->0i 

-.100000E«01 

0.10M60E'>01 

4 

4 

17 

10 

1 

1 

- .  100000E'>01 

0  100000E«81 

-.100000E*01 

0.1000ME«01 

4 

4 

18 

10 

1 

1 

- . 100000E«01 

0  i00000E*ei 

-.100000E«01 

0.1000e0E«01 

4 

4 

19 

10 

1 

1 

-  100000E«01 

0  100000E>0I 

- . 100000E*01 

0.100000E«01 

4 

4 

20 

10 

1 

1 

-.100000E<01 

0.100000E*01 

-.100000E«01 

0.100000E'>01 

4 

4 

21 

10 

1 

1 

- .  100000E'>01 

0.100008E«01 

-.100000E«01 

0.100000E«01 

4 

4 

22 

10 

1 

1 

-  100000E«01 

0.100000E<01 

-.100000E«01 

0.100000E-01 

4 

4 

23 

10 

1 

1 

- . 100000E«01 

0.100000E<01 

-.100000E«01 

0. 100000E«01 

4 

4 

24 

10 

1 

1 

-.100000E«01 

0. 100800E'>01 

-.100000E-01 

0.100000E*01 

4 

4 

26 

10 

1 

1 

- . 100000E«01 

0.10e000E*01 

-.100000E«01 

0.100000E*01 

4 

4 

26 

10 

1 

1 

- . 100000E>01 

0.100080E«01 

-.100000E«01 

0.100000E«01 

4 

4 

27 

10 

1 

1 

- .  100000E'>01 

0.100800E->01 

-.100000E<01 

0. 100000E«01 

4 

4 

28 

10 

1 

1 

-.100000E-01 

0.100000E«01 

-.100000E*01 

0.100000E«01 

4 

4 

29 

10 

1 

1 

-.100000E<01 

0  100000E'>ei 

-.100000E'>01 

0.100000E«01 

4 

4 

30 

10 

1 

1 

-.100000E'>01 

0.100000E«01 

-.100000E»81 

0.100000E«01 

4 

4 

31 

10 

1 

1 

-.100000E401 

0  100000E«01 

-.100000E«01 

0.100000E«01 

4 

4 

32 

10 

1 

1 

-.100000E-01 

0  i0000eE'>0i 

-.100000E«01 

0.100000E'»01 

4 

4 

33 

10 

1 

1 

- .  100e00E'>01 

0.100000E«01 

-.100000E'>01 

0.100000E«01 

4 

4 

34 

10 

1 

1 

- .  100000E'>01 

0.i000e0E«0i 

-.100000E'>01 

0.100000E«01 

4 

4 

36 

10 

1 

1 

-.100000E'>01 

0  100000E«01 

-.100000E«01 

0.100000E«01 

4 

4 

36 

10 

1 

1 

- .  100000E'>01 

0.i00000E«ei 

-.100000E>01 

0.100000E«01 

4 

4 

37 

10 

1 

1 

-.100000E-01 

0  i00000E«ei 

-.100000E«01 

0.100000E«01 

4 

4 

38 

10 

1 

1 

- .  100000E«01 

0  100000E«01 

-.100000E«01 

0.100000E*01 

4 

4 

39 

10 

1 

1 

-.100000E'>01 

0.i000e0E«0i 

-.100000E«01 

0.100000E*01 

4 

4 

40 

10 

1 

1 

- .  1000e0E'>01 

0.100000E«01 

-.100000E'>01 

0.100000E*01 

4 

4 

41 

1 

3 

3 

0.000000E->00 

0.608080E-01 

0.000000E'>00 

0.508000E-01 

4 

4 

REGION  -  TRANSLATION  (T1,T2,T3),  IZAX 


1 

0  N0000E<00 

0.000000E*00 

0 . 000000E«00 

3 

2 

0.000000E>00 

0  000000E«00 

0.2e6000E-01 

1 

3 

0.000000E*00 

0.000000E«00 

0.286000E-01 

2 

4 

0.000000E>00 

0.000000E>00 

0.000000E<00 

3 

6 

0.000000E->00 

0.000000E«00 

0.000000E<00 

3 

6 

l.0fl0000E«0e 

0.000#00E«00 

0.000000E«00 

3 

7 

0.0M000E'>00 

0  000000E«00 

0.000000E<00 

3 

8 

0  0f000fE>00 

0.000000E«M 

0  000000E->00 

3 

9 

0.f0ff00E>0<°: 

0  000000E<00 

0.000000EO00 

3 

10 

0.000f00E<00 

0.0f0000E«00 

0.000000E*00 

3 

11 

0.000f00E«00 

0.000000E<00 

0.000000E*00 

3 

12 

0.0f00ME<00 

0.000000E«00 

0.000000E«00 

3 

13 

0.00W00E»00 

0.H0M0E^00 

0.000000E«00 

3 

14 

0.00M00E«00 

0.000M0E-00 

0.000000E«00 

3 

16 

0.000000C«H 

0  N0000E>00 

0.000000E«00 

3 

16 

0.006000E^H 

0  000000E«00 

0.000000E<00 

3 

17 

0.0f0M0E>00 

0  000000E<00 

0  000000E«00 

3 

18 

0.000000E«W 

0  000000E<00 

0.000000E*00 

3 

19 

0.000fME<M 

0  000000E*00 

0.000000E«00 

3 

20 

0.0M00eE«00 

0.000000E«00 

0.000000E«00 

3 

21 

0.0f0fME«00 

8  000000E«00 

0.080000E-00 

3 

22 

0.000#HE«00 

0  000000E<00 

0.000008E«00 

3 

23 

0.00M00E«00 

0  000000E«00 

0.000000E«00 

3 

24 

0.000000E«00 

8.000000E«00 

0  000000E'>00 

3 

26 

0.HNHE«00 

0.000000E*00 

0.000000E'>00 

3 

26 

0.H0000E<06 

0.000000E-00 

0.000000E'>00 

3 

27 

0.900000E«00 

0  000000E<00 

0  0000e0E«00 

3 

28 

0.000000E>00 

0  000000E«00 

0.000000E<00 

3 

29 

0.6f0000E<00 

9.000000E«00 

0.000000E<00 

3 

30 

0.000000C«00 

0.000000E«00 

t.Btems.-ee 

3 

31 

0.000000E«00 

0.000000E«00 

0  000000E«00 

3 

32 

0.000000E«H 

0  000000E<00 

0  000000E'>00 

3 

33 

0.000000E«00 

0.000000E<00 

0  000000E<00 

3 

34 

0.00M00E'00 

0  000000E>00 

0  000000E«00 

3 

36 

0.00f000E<00 

0.000000E->00 

0.000000E<00 

3 

36 

0.0M000E<00 

0  000000E>00 

0.000000E«00 

3 

37 

0.000N0E«00 

0  000000E«00 

0.000000E<00 

3 

38 

0.000080E«00 

0  000000E«00 

0.000000E'>00 

3 

39 

0.000000E«00 

0  000000E<00 

0  000000E<00 

3 

40 

0  0e0000E«00 

0  000000E>00 

0.00000eE<00 

3 

41 

0  000000E<00 

0  0000e0E<00 

0  000000E>00 

-3 

28 
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REGION  -  CCS 

1  0.952&E-02  ,0.9w'  L,  ej  0.e030E‘0«  ?.;)92«E.00  0.0000E*00  0.0000E-00  I.SeME^H  •.•000E<'M  e.0000E.00 

2  0  50B0E-0;  0  3s  0  0.0000Es00  0  0a00E'B0  0.000eE«0e  0.B00eE«BB  e.BeB0E«M  •.••bbe^m  •.0000E'>00 

3  0.B030E-01  3  9£  A‘  3  0^  'v:  :  :  0000Es0«;  0  0000E«00  f.M00E*M  •.00B0E>00 


4  0.2540Es03  ‘  91'  -  s  g  0  j.,'9;30Ls00  i'.93»)0E>00  0.000eE>H  0.032SE*M  I.MBBE^M  ■.0000E«00 

5  0.1000ES01  0.2>)00>'..  3:  '1.3'9's  :  'i.B000Es01  0.£000E<0?  C.000BE>IW  e.lBME^BB  •.MME«N  f.HeBE^fB  9.MeBE^M 

6  0.2000E*01  0.2  '.;  <  .r,..,.  s  3.e000E<00  0  0000 E* 00  0.0000E«00  0.g000E«00  B.BfBBE-'H  0.0000E>00 

7  0.3000E>01  3  4  4  .'  VCs  .:  9  :000Es02  0  000eE>e0  0.0000E«00  0.00e0E«0e  0.0000E«00  0.0e00E«H  0.0000E>00 

8  0.4e00EsOl  3  b023E'  '  C.:'-O0'  s0?  Z.:?v0E<02  0.0000E<00  0.0000E*00  B.BBBBE^BB  0.0000E*00  0.t000E«H  B.BBBBE^BB 

9  0.S000Es0'  C.C.?-'9r-r.  E.  J:  «  123?Es02  0.00e0E<00  e.0090E>00  0.0000E'>00  0.0000E«00  0.0000E-00 

10  0  ‘"’'.■Cc«01  3  -A'.  3.:300E*02  0.0000ES00  0.0000E>00  0.000#E«BB  0.M00E«00  0.0000E«00  0.0000E>00 

11  0.B000E<01  3  ] 9008-02  0.0000E-00  B.eeME^tt  0.0000E'>00  B.BBBBE^Bg  0.0000E'>H  0.0000E<00 

12  0.90008*21  3  IC  :■  ■  '  -^2  ''•?!"*02  g  « '0E.30  0.0006E.00  0,0000E*00  0.0000E«00  0.0000E-00  0.0000E-00 

13  0  1300t*01'  2  :  -  :...!^L*02  0.00e8E<0e  0  0000E«00  0.0000E*g0  0.0000E«e0  a.»B9»E*U  0.0000E-00 

14  0.n90E.B2  ri.l.  .'2  ::.i8:38E.02  0.0000E-00  0.0000E-00  0.0000E'.00  0.0600E-00  0.0000E-00  0.0000E-00 

15  0  1230E'02  e  s  21  *1  *  '  .-.1930E<02  0.0000E*00  0.0000E'>00  0.0000E«0e  e.BeBeE-^BB  0.0000E«00  0.0000E-00 

16  0.1300E*02  0  ]'•.<.  .  **.  2  2';20**12  1!.2000E*02  0.0000E>00  0. 00008*00  0.00008*00  0.0000E*00  0.0000E*00  0.0000E*00 

17  0.1500E*02  0  b<<  -  '  :!.2*»'  6.2200E*82  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

18  0.16008*02  2  ivl'sbi  •*  /  3.23008*02  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

19  0.17008*02  3  'b  *\  1.2  *'12  3.2<008*92  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

20  0.18008*02  ‘1  •"  V..;*i  -  -’-i  j.2;*',)0**02  0. 00008*00  0  00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

21  3  19008*02  ?f‘  *>  2  •.'8*02  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

22  0.20008-62  3.2itr  :-.  ■/  5.2, '098*02  6.00038*09  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

23  0.22008*92  0.230'  :-*:2  3.3399'/*t2  .5  2 2908-02  9.03008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

24  0.23008.-02  9  lia.v:* 5',:  *.3199'M'2  3  ;-:f998-32  0.0000L-00  0.00008-00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

25  0.24008*02  9 . .BEaO**!'  2.31068-02  0  00068-00  0  00008-00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

28  0.26008*02  Jl.CSiti— u'.3;;0f  8-,  ’  •'  ,32308*92  6  00008*00  0.00008*00  0.00008*00  0.00ME*00  0.00008*00  0.00008*00 

27  0.26008-02  0.27«9E->2  » 3.3308*02  0.00008-00  0.00008-00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

28  0.27008*02  2  2'102’'-8:'  1  3!30~-r2  8  34068-02  0.00008-00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.06008*00 

29  0  29008-?2  ■'  '  ■  9.;  ■  '2  r-;i‘3f*0?  . *'950E-9,1  (1.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

30  0.30098-62  3 . .'I- .  1.37008-02  e.00008'e0  0.00008-00  0. 00008-00  0.00008-00  0.00008-00  1.00008-00 

31  0.31008-02  .!.32-'5,;'  1. 50308-02  3.80008-00  0.00008-00  0.00008-00  0.00008-00  0.00008*00  0.00008*00 

32  0.32008*02  9.  •'.•i- .?7  "..32'j?8*92  3.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008-00 

33  3.33308*02  6, 3491  >.'2  9. •! : 'C'.-f: 2  0  43<5e£*02  9.90008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

34  0.34008*02  9.553'3''.k1:  9.4'2:M..-02  •'.41908*02  0.00008*00  0.00008*00  0.00008*00  1.00008*00  0.00008*00  0.00008*00 

35  0  33008*02  •9.3r-208'92  9.4«.K.1*32  3  43968-' 02  0.00008*30  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

36  0  3/008*02  0  31098-02  9.450e:.-<.2  <'.44068*02  6.60908*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 

37  0,38008*02  9  9,4e00i--0?  9  4  5(?0E-02  «  0000E-00  9  00908*00  0.00008*00  0.00008*00  0.00008*00  0.00008-00 

36  0  39008*02  J  •  :  9  :•  -  7  . «4''0r:. 02  9  90008*00  0 . 00008*00  0. 00008*00  0.00008*00  0.00008*00  0.00008*00 

39  0.40008*02  0.4.h^f->:-02  3  4eo0r:.,’2  •’.47008*02  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008-00 

40  0,41008-02  0  429»! ->)•>  -I  >>  49008.-02  0.90008-00  0.00008-00  0.00008*00  0.00008*00  0,00008*00  0.00008*00 

41  0.28808-01  9  ■1.499C8-0'-  «  48008*02  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00  0.00008*00 


NO,  OF  SURFACE  A5E>E  .22  N  : ,  OF  SL-Rl' ACE  Fia.0  POINTS  =  122 

1 

Pi-;;.  ■.  i  -,188  RUN  »EDG  DATE  23-0CT-9 


FREQUENCY;,  2090  0  i(‘00  »  >>  690  6 

SYMMETRY  rYPE=  RFF  Nl  !9!:  ;i;  Y2,  4 

VELOCITIES  (BLKf.AFt-*,  i  .  ,M  ,  IF,-,  ‘NO', 


1 

0  100009F-01 

0  00 

1 

2 

i-  . 

0  iP3C00t-01 

0.00 

1 

3 

y , 

0  1000008-01 

0  00 

1 

4 

V  . 

6  1096008-02 

0  00 

1 

b 

0  1096008-01 

0.00 

1 

tj 

{ 

:  ’  A'  , 

0  1000098-01 

0.00 

1 

7 

f'. 

I,/'  *3, 

0.1060008-01 

0.60 

i 

8 

i 

6  1096098-01 

9.60 

1 

G 

j  .,1’ 

6  1090098-01 

9  00 

1 

1  *5 

r 

0.',P30''>J;'.01 

0.90 

1 

i  i 

e 

;  7y.  - 1^:7'  •■d': 

t z*  ?  •? 

0.1006008-01 

0  00 

1 

1/ 

R 

0  1096098.01 

0.00 

i 

n 

e 

9  !69600':.0I 

0.00 

1 

14 

e 

■^5- 

9  1036008.01 

0  00 

1 

I  b 

R 

] ' 'y. 

3  i03e00':.0i 

6.90 

1 

16 

i . 

.1; 

0  '006008-01 

6.00 

9  '  ’•  ■.■>pr,.7j  j  '096091-01  0  00 

>  ;  •"  '■  '•  '  'If -1  i  0,00 


1 


Thompson,  Dubus,  ft  Siders 


1 

19 

6.ie6M6E«61 

0.666000E*00 

0.100006E'>ai 

6.06 

1 

2« 

6.166666E«61 

0.000e06E*e6 

0.10060eE-01 

6.00 

1 

21 

6  166666E«61 

6.066006E«66 

0  100000E«01 

0  H 

1 

22 

6.ie6666E>61 

6.609e06E*06 

0.100066E>01 

0.00 

1 

23 

6.166666E*61 

0.606006E«06 

0.10e666E<01 

0.M 

1 

24 

e.l66666E«01 

0.6e6006E^60 

0.100660E-01 

0.00 

1 

2S 

i.l33e96E«6l 

e.606e06E>e0 

0  166660E'>01 

0  00 

1 

26 

6.ie66ME'>61 

6.666636E>60 

0.100066E>01 

0.00 

1 

27 

6.1666ME«61 

6.606606E*66 

0  160066E«01 

0.N 

1 

28 

6.186e66E*61 

0.666606E«66 

0  106666E«01 

0.00 

1 

29 

6.166606E«61 

0.006606E*06 

0.18066eE<01 

0.H 

1 

38 

6.ie6e66E'>ei 

0.066000E«06 

0.100060E'01 

0.00 

1 

31 

6.166006E'>01 

0.006006E«66 

0.160000E.01 

0.00 

1 

32 

6.16666eE«01 

0.600006E«00 

0.100060E*01 

0.00 

1 

33 

6.066000E«00 

6.000006E«6e 

0.060006E>00 

0.00 

1 

34 

0.e6e66eE<06 

6  H6006E-00 

0.060000E«00 

0.00 

1 

36 

6.060e00E<66 

0.066606E<06 

0  068666E<00 

0.00 

1 

36 

6.066M0E<96 

0  060606E«06 

0.800060E«00 

0.00 

1 

37 

6.660666E«00 

6.666666E<00 

0.060660E<00 

0.00 

1 

38 

6  800600E«00 

0.006006E*00 

0.000000E>00 

0.00 

1 

39 

6.660606E«06 

6.600660E->00 

0.606606E'.00 

0.00 

1 

41 

6.666006E<06 

0  0060e6E<06 

0.660000E«00 

0.00 

I 

41 

6.600666E«66 

6  666666E«06 

6.666606E-00 

0  N 

1 

42 

6.666006E<M 

6.666006E«06 

0.660066E'>8e 

6.N 

1 

43 

6.666666E«60 

0  66M66E'>66 

0.000666E«e0 

0  M 

1 

44 

6.6666#6E>M 

6  6066e0E<06 

0.606666E>00 

0  00 

1 

46 

6.666666E<66 

6.666606E'>60 

V .  vtnwv0c-*v0 

0.00 

1 

46 

6.6H666E>66 

6.666606E«66 

0.0066e6E<00 

0  00 

1 

47 

6.666666E'>N 

0.666606E«60 

0.ee6666E-00 

0  00 

1 

48 

6.666666E«66 

6.666606E«66 

0.0666eeE*00 

0.H 

1 

49 

6.666666E«66 

6.666666E«00 

0.660666E-80 

0.N 

1 

66 

6.6666NE«N 

6.666606E«66 

0.666066E-60 

0.M 

1 

61 

6.666666E«66 

6.006606E^66 

6.606606E«e6 

0.00 

1 

62 

6.666f66E«H 

6.666e06E«N 

0.666660E«00 

0  M 

1 

53 

6.6666f6E«66 

6.606606E->N 

0.N 

1 

64 

6.6666HE«66 

6.666006E>66 

0.666666E*60 

0.M 

1 

66 

6.06666«f  60 

6.660606E^66 

6.660660E>00 

0.00 

1 

56 

6.66660Cc<66 

0.666006E'>66 

0.660066E*60 

0  00 

1 

67 

6.666666E«06 

6  e000e6E<00 

0.660666E«e0 

0.00 

1 

66 

0.6e60e6E«66 

6  666006E<06 

0.000666E<e0 

0.M 

1 

69 

0.0f0660E>66 

0.e066e0E<06 

0.006008E'>0e 

0.00 

1 

66 

6.e6H66E«66 

0.606600E*66 

0.066666E<00 

0.00 

1 

61 

0.606666E«e6 

6.060006E«66 

0.6e6068E'>60 

0.H 

1 

62 

6  666666E«06 

0.606e00E->H 

0.660006E«60 

0  00 

1 

63 

6.666666E<06 

6  0He06E>60 

0  6e6666E*e0 

0.H 

1 

64 

6.006006E<e6 

0.6000e0E«H 

0.008808E>00 

0  60 

1 

66 

6.666666E>06 

6.666e06E«66 

0.000e00E«00 

0.00 

1 

66 

6.606066E^6e 

6.060606E'>N 

0.000e06E<00 

0  00 

1 

87 

6.066666E«N 

6.666e06E*N 

0.000606E*0e 

0.00 

1 

68 

6.600606E«M 

6.666606E«N 

0.6666e0E'>66 

0.00 

1 

89 

6.66f660E>66 

6  ee6606E>66 

0.000680E*06 

0.N 

1 

76 

6.666666E*N 

6.6H600E*66 

0.660066E*0e 

0.00 

1 

71 

6.666666E>e6 

6.666e06E>00 

0.000060E«00 

0.M 

1 

72 

6.6060e6E«06 

6.666600E«06 

0.0006e0E«08 

0.00 

1 

73 

0.6f66NE*06 

6.60e0e0E«66 

0.0e006eE«00 

0.00 

1 

74 

6.666M6E«e6 

6  660006E«66 

0.660068E«e0 

0.H 

1 

76 

6.6f6060E«66 

6.e00666E*60 

0.680668E«0e 

0.M 

1 

76 

6.666M6E«H 

6.60600eE«66 

0.000686E<80 

0.H 

1 

77 

6.666M0E«66 

6.606e06E«N 

0.0600nE*0e 

0.00 

1 

78 

6.e66666E«66 

6.666606E*60 

0.006606E<06 

0  H 

I 

79 

6.6060e0E«66 

6.006600E«66 

0  000006E<80 

0  N 

1 

86 

6.066666E*06 

6.666606E<06 

0  066066EO00 

0.M 

1 

81 

6.66606eE>06 

0  e0e006E«66 

0.00H08E<00 

0.H 

1 

82 

6.066666E>06 

6.H6606E<H 

0.606866E«06 

0.00 

1 

83 

8  666666E«M 

6.666606E'>00 

0.060660E‘H 

0.00 

1 

84 

6.60666eE«H 

6.660606E«66 

0.066666E*n 

0.H 

1 

86 

6.66M60E*66 

6.606606E<66 

0.066006E<66 

0.H 

1 

86 

6.666608E«66 

e.666660E«06 

0.606668E'>06 

0.M 

1 

87 

6.6666e6E'>06 

6.0f66e6E*60 

0.606606E«n 

0.00 

1 

88 

6.666660E«00 

0  609006E«0e 

0  606000E«66 

0.H 

1 

89 

0  066e69E>00 

0.006e00E*H 

0  000600E<00 

0  00 

1 

96 

0  666660E*00 

0  00e000E'>00 

0  000666E<06 

0.H 

30 
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1 

91 

0.000000F'00 

0  000000E*00 

0  0e6000E<80 

0  00 

1 

92 

0.000000E*00 

0.000000E*00 

0  000000E<00 

0  00 

1 

93 

0.003000E>30 

e.000000E<00 

0  000000E>60 

0  00 

1 

94 

0.000000E>00 

0  00e000E'>00 

0  000000E<00 

0.00 

1 

95 

0.000000E>80 

0.000000E-00 

0  000000E«00 

0  00 

1 

96 

0.000000E<00 

9.030000E'>00 

0  000000E-00 

0.00 

1 

97 

0.000000E-00 

0.000000E'>80 

0  000000E>00 

0  00 

1 

98 

8.000000E>00 

0.000000E<e0 

0 . 000000E<00 

0  00 

1 

99 

8.000000E<00 

0  000000E«00 

0  090080E«80 

0  00 

1 

100 

0.000000E<00 

0 . 000000E'>00 

0 . 000000E.00 

0.00 

i 

101 

0.000000E*00 

0.000000E«00 

0.000000E<00 

0  00 

1 

102 

0.000000E<00 

0.000000E«00 

8 . 090000E.00 

0.00 

1 

103 

0  000000E<00 

0  000000E-00 

0  000000E*e0 

0.00 

1 

104 

0 . 000000E'>00 

0 . 000000E«00 

0.00 

1 

105 

0.000000E'>00 

0 . 000000E-00 

0  0e0000E'>00 

0.00 

1 

106 

0.0e000eE'>00 

0 . 000000E-00 

0.000000E«00 

0.00 

1 

107 

0.0000e0E<00 

0.000000E«00 

0.00080eE«00 

0.00 

1 

108 

0.000000E<00 

0.000060E-00 

0.00000eE«00 

0.00 

1 

109 

0.000000E*00 

0.000000E>00 

0 . 000000E«00 

0.00 

1 

110 

0 . 000000E>00 

0 . 000000E<00 

0 . 000000E*00 

0.00 

1 

111 

0.000000E«00 

0.000000E<00 

0 . 000000E<00 

0.00 

1 

112 

8.000000E<08 

9 . 000000E.00 

3  000000E>00 

0  00 

1 

113 

0 . 000000E>00 

0.000000E*00 

0.000000E*00 

0.00 

1 

114 

-0.919748E-01 

-0.681387E-01 

0  114495E-00 

-143.45 

1 

115 

-0.710678E-01 

-0.344366E-01 

0.789716E-01 

-154.15 

1 

118 

-0.495600E-01 

'0.135340E-01 

0  513755E-01 

-184.73 

1 

117 

-0.710071E-01 

-0.344389E-01 

0.709711E-61 

-164.15 

1 

118 

-0.583367E-01 

-0.191811E-01 

0.614036E-01 

-161.82 

1 

119 

-0.43460IE-01 

-0.a47442E-e2 

f.442788E-81 

-168.97 

1 

120 

-0.495S83E-01 

-0.135344E-01 

0  513732E-01 

-164.73 

1 

121 

-0.4346aS£-01 

-3,847467E-02 

0  442771E-01 

-168.97 

1 

122 

-0.353640E-01 

-0.477029E-02 

0  3S6843E-01 

-172.32 

*.»  FULL  VaOCITY  SYWETRY 

PROGRAM  C 

H  I  £  F  8  8 

RUN  fEDG  DATE  23-OCT- 

FREHUENCY=  2000.0  RHO=  1000 

SYMMETRY  TYPE=  RF  NO.  OF 

SURFACE  PRESSURES  (BLKf,AREAf, 

.0  C=  1600.0 

BLOCKS^  4 

REAL, IMAG, MAG, ANQ) 

1 

1 

0.190689E*05 

0.104885E«06 

0.106604E>0e 

79.70 

1 

2 

0.190689E-06 

0.104885E<0e 

0. 106e04E«06 

79.70 

1 

3 

0 . 190689E«06 

0.104885E«06 

0.106604E>06 

79.70 

1 

4 

0 . 190689E-05 

0.104885E-06 

0  106e04E'>66 

79.70 

1 

5 

0.ie0495E«0S 

0.ie5601E<06 

0 . 107338E-06 

60.32 

1 

6 

0.180497E*06 

0.106801E'>06 

0  107329E-0e 

80.32 

1 

7 

0.180497E-05 

0.105801E'>06 

0.107330E>06 

80.32 

1 

8 

0.180496E-05 

0.105801E-06 

0  107330E*06 

80  32 

1 

9 

0.165834E-e5 

0.108442E*08 

6.107723E*08 

81.16 

1 

10 

0.16S638E'>05 

0.106442E*06 

0  107723E*08 

61.15 

1 

11 

0.166638E.05 

0.106442E-06 

0  107723E-06 

81.15 

1 

12 

0.165634E-66 

0.10e442E«06 

0  107723E*06 

81.16 

1 

13 

0.153097E-05 

0,107242E-06 

0  108287E-06 

62.03 

1 

14 

0.150101E«05 

0.107241E-06 

0.108207E*08 

82.03 

1 

16 

0.150101E*05 

0.107241E*06 

0. 108266E«06 

82.03 

1 

16 

0.150096E*05 

0.107241E-06 

0. 168287E*66 

82.03 

1 

17 

0.136903E'>05 

0.108067E-06 

0.I06931E«06 

82.78 

1 

18 

0.138905E*05 

0.1080eeE*06 

0.108930E<06 

82.78 

1 

19 

0.13690bE-05 

0.108066E<66 

0.108930E«06 

82.78 

1 

20 

0.138902E-05 

0.100066E-06 

0 . 108930E-06 

82.78 

1 

21 

0.127172E»05 

0.100774E*06 

e.ie9515E-06 

83  33 

1 

22 

0.127172E.05 

0.108773F*06 

0  109514E-06 

83.33 

1 

23 

0.127173E*06 

0.108774E*06 

0.10951SE.06 

83.33 

1 

24 

0.127173E-05 

0.108774E«06 

8.109515E*08 

83.33 

1 

25 

0.1209eeE-05 

0.109246E«0e 

0.1099I4E*06 

83.68 

1 

26 

0.120969E<06 

0.ie9247E.06 

0.109914E*06 

83.68 

1 

27 

0.120969E-0S 

0.109247E-06 

0  109914E*0e 

83.68 

1 

26 

0.120969E-06 

e.l09247E.06 

0  109914E-06 

83.68 

1 

29 

0.nai85E-06 

0.10892eE.fe 

0  109689E-06 

83  79 

31 


Thompson,  Dubus,  k  Siders 


1 

3< 

6.118464E«65 

e.l66926E>66 

0.10956gE<06 

83.79 

1 

31 

6.118485E«66 

6  16B926E*M 

0.169569E<06 

83.79 

1 

32 

6.118485E<6S 

6.166926E'>M 

0.169569E-06 

83.79 

1 

33 

6.732883E>64 

6.258366E>64 

0.777672E-04 

19.42 

1 

34 

6.585173E*64 

6.168433E«e4 

0.695135E>04 

10.60 

1 

35 

6.486771E-64 

6.886709E-B2 

0.486852E>84 

1.04 

1 

36 

6  4M883E«64 

-6.821921E<e3 

0  413388E>04 

-8  66 

1 

37 

6.346278E*64 

-6.114628E*B4 

0.369M7E«04 

-18.62 

1 

38 

6.276665E«64 

-6.153748E<64 

0.316378E>04 

-29  06 

1 

39 

6  7M366E'>64 

6.229893E->64 

0  742824E»04 

18.03 

1 

46 

6.572969E«64 

e.962306E*03 

8.S80994E>04 

9.63 

1 

41 

6.486182E«64 

6  277307E<02 

e.486190E«04 

0.33 

1 

42 

6.464664E>64 

-6.656642E«63 

6.48Sa82E<84 

-9.20 

1 

43 

6.33759eE«64 

-6.116E42E*64 

6.367140E«04 

-19.06 

1 

44 

6.274827E*64 

-6.154748E>64 

6.31E406E<04 

-29.38 

1 

46 

6.664396E«64 

6.192619E«64 

8.691838E-B4 

16.19 

1 

46 

6.553895E«64 

6.781282E«03 

e.55937BE«e4 

6.03 

1 

47 

6.469735E*64 

-0  7B2261E«e2 

0.489787E*04 

-0.66 

1 

48 

6.396668E«64 

-6.7124S6E«03 

0.4e4393E'>04 

-10.16 

1 

49 

6.333319E'>64 

-6.119883E'>64 

0  354222E*04 

-19  78 

1 

56 

6.272385E«64 

-6.156E16E-64 

0  314148E<04 

-29  88 

1 

51 

6.733214E«64 

0.25840eE«04 

0.777417E*04 

19.41 

1 

52 

6.766e36E«64 

6.229974E«64 

0  743111E^04 

18.03 

1 

53 

6.664566E<64 

6.192976E->64 

0.692666E<04 

16.19 

1 

54 

6.586742E«64 

6.168616E«64 

0.595728E<04 

10.61 

1 

66 

6.673438E'>64 

6.963e04E«63 

0.581481E<04 

9.64 

1 

56 

6.654166E«64 

0.782138E->63 

0.559659E'>04 

8.03 

1 

67 

6.487716E.64 

6.918S27E«62 

0  4B78e3E^04 

1.08 

1 

68 

6.48e959E«64 

6.3B3ei6E«02 

0.4669e8E*04 

0.36 

1 

59 

6.476178E'>64 

-6.6B7496E«62 

0.470228E«64 

-0.84 

1 

61 

6.416283E-64 

-6.616467E«63 

0.41488eE>04 

-8.64 

1 

81 

f.46S918E'>64 

-6.6E1134E«63 

0.411107E<84 

-9.11 

1 

62 

6.398816E«64 

-6.769B96E«63 

0.465664E«04 

-16. M 

1 

63 

6.343146E«64 

-6.113622E>64 

e.38i4e8E«e4 

-18.32 

1 

64 

6.339943E-64 

-6.11S717E-64 

0.359098E«e4 

-18.80 

1 

66 

6.3346S2E'>64 

-6.119413E'>64 

6.365319E*04 

-19.64 

1 

66 

6.282i7  E»64 

-6.161727E-64 

0  326379E«04 

-28.27 

1 

67 

6.279527E-64 

-6. 163062E«64 

0.3ie696E<84 

-28.70 

1 

68 

6.276683E«64 

-6.156636E*64 

0  316M6E*04 

-29.48 

1 

69 

-6.224617E-.63 

-6.2N137E<64 

0.201393E«04 

-96.40 

1 

76 

-6.1999e4E«63 

-6.2N113E«64 

0.261116E*e4 

-96.71 

1 

71 

-6.136684E<63 

-6.2N526E*64 

0.26M52E-04 

-93.73 

1 

72 

•6.199928E«63 

-6.2N538E«64 

0  261E33E«04 

-96.69 

1 

73 

-6.1774E6E-63 

-6.2N476E«64 

0.201269E«04 

-95.  W 

1 

74 

-6  iHsasE^es 

-6.2N746E«64 

6.2B1667E>04 

-93.27 

1 

76 

-6.146792E<63 

-6.262964E«64 

0.263487E>64 

-94.11 

1 

76 

-6.127149E<63 

-0.262835E*64 

0.203233E«04 

-93.69 

1 

77 

-6.7e9867E*62 

-6.262799E«64 

e.2e2945E«04 

-92.17 

1 

78 

6.23f684E*64 

-6.178769E*64 

0.2NS39E'>04 

-37.47 

1 

79 

6.186ei6E<64 

-6.192686E>64 

0.268366E«04 

-46.88 

1 

86 

6.145415E«64 

-6.263076E<64 

0.249776E*04 

-64.40 

1 

81 

6.1N272E«64 

-6.26e6NE«64 

6.234116E«04 

-63. N 

1 

82 

6.692N9E«63 

-6.2W777E*64 

0  220S9gE'>84 

-71.74 

1 

63 

6.336647E<63 

-6.267686E«04 

0.209795E«84 

-80  78 

1 

84 

6.22e388E'>64 

-6,177866E*64 

6.289318E<64 

-37.87 

1 

85 

6.ie5141E«64 

-6.193147E*64 

0  267650E«64 

-46.21 

1 

66 

6.144146E'>64 

-6.263328E«64 

0.249240E*04 

-64.67 

1 

87 

6.166353E«64 

-0  268721E*64 

0.233663E«64 

-83.22 

1 

88 

6.68e285E«63 

-6.269626E'>64 

8.2267S8E<64 

-71.89 

1 

89 

f.333962E<f3 

-6.2#7628E-64 

8. 289764E'f4 

-80.84 

1 

96 

6.223927E«64 

-6.179368E«64 

6.286B76E«64 

-38.69 

1 

91 

6.181596E'>64 

-6.194219E«04 

0.26669eE<04 

-48.92 

1 

92 

6.141376E-64 

-6  264026Eo64 

0.24822IE*64 

-66.28 

1 

93 

6.163369E«64 

-6.269267E«64 

6.23337aE>B4 

-63.73 

1 

94 

6.672972E.63 

-6.216494E'>64 

6  2269ME«04 

-72.27 

1 

96 

6.32e846E«63 

-6.26821SE«04 

0.210796E«04 

-81  03 

1 

96 

6.214264E«64 

-6.183116E'>64 

0.26le62E<64 

-40.63 

1 

97 

6.171256E«64 

-6  197426E«64 

0.261356E«04 

-49  M 

1 

98 

6.131466E<64 

-6.266216E<84 

0.244658E<04 

-67.48 

1 

99 

6  941758E«63 

-6.216342E<64 

0  230462E*64 

-66.86 

1 

IN 

6.69N66E>63 

-6.216391E«64 

0.218526E>04 

-74.32 

1 

161 

6.2666E7E<63 

-6  266788E«64 

0  208374E«04 

-82.92 
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1 

102 

0  214010E.04 

-0  183033E-.04 

0.281611E*04 

-40,54 

1 

103 

0.169623E«04 

-0,197963E-04 

0  260693E-e4 

-49.41 

1 

104 

0.129266E<04 

-0  206865E*04 

0 . 243932E<S4 

-58.00 

1 

105 

0  917170E.03 

-0  210938E'>04 

0  230015E>04 

-66.50 

1 

106 

0  564417E-03 

-0  210880E«04 

0  218303E«04 

-75  02 

1 

107 

0.226739E.03 

-0  207149E*04 

0 . 2083e6E«04 

-83  75 

1 

108 

0.214650E-04 

-0.102729E-04 

0.281894E<04 

-40.41 

1 

109 

0.169558E<04 

-0  19e024E«04 

0 . 260698E>04 

-49.43 

1 

110 

0.128817E-04 

-0.207080E.04 

0  243877E.04 

-58.12 

1 

111 

0.910571E.03 

-0  2U201E-04 

0  229994E-04 

-68  68 

1 

112 

0.556222E.03 

-0  211142E-04 

0  218348E-04 

-75  24 

1 

113 

0.21631!E*03 

-0  207386E-04 

0  208511E*04 

-84  05 

1 

114 

0  275924E-05 

0  204566E<05 

0  343405E.05 

36.55 

1 

115 

0.213203E*05 

0.103310E<05 

0  236915E»05 

25.85 

1 

116 

0  140682E*05 

0  406019E«04 

0.16412eE-05 

15.27 

1 

117 

0.213201E*05 

0  103311E-05 

0  236913E<0S 

25.85 

1 

118 

0  175010E-05 

0  574834E-04 

0.184209E<85 

18.18 

1 

119 

0.13e380E-05 

0  254233E*04 

0 . 132e36E->05 

11  03 

1 

120 

0.148875E.05 

0  40S031E>04 

0  154120E»05 

15  27 

1 

121 

0.130376E*06 

0.254240E<04 

0  132631E'>05 

11.03 

1 

122 

0.106092E^05 

0.143109E>04 

0  107eS3E*05 

7  68 

♦  .»  FULL  VaOCITY  SYVWETRY 
POWER  OUT{KW)=  0.003 
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Sa«pl«  Rao  3;  Conical  Source  With  Circular  Planar  Baffle  Radiation  Matrix 

The  ge^/inetry  for  this  problem  was  generated  using  Di‘5PLAY3D^,  which  is 
available  to  all  CHIEF  users.  Run  3  cciputes  a  radiation  matrix  using  both 
the  subroutines  ZRADUX  and  ZRADUl  since  the  circular  planar  baffle  has  an 
impedance  coating.  Also  the  cone  is  a  new  surface  in  CCUNMD. 


2  ,,3 - L.  X 

Fig.  4.  A  rkdiatiog  coce  with  a  circular  planar  baffle. 


Program  Listing  of  Sample  Run  3 

C  Prograi  CONE 

C  Exaipla  Problaa  of  a  rone  topped  by  an  iapedance  surface 

C  (COMUON  VARIABLE  DECLARATIONS  GO  PERET 

COMPLEX  PMATX(  152),  VaATX(  152) 

M0Si7E  =  l.,2 

RUN ID  =  ’C0NE‘ 

DATE  =  ■25-FEB-9r 

CAa  INITCJ 
CALL  OPNC'^ 

RHO  =  1»0« 

C  =  1500 


Ctaaaattaaaaaaa*  iapedance  coating  aodit ications  «•«***««««*«••* 

NUMZ  =  4 

N  =  1 
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C - clear  mpedapce  fv.r  The  Cons 

DO  10  1  -=  N,  N  •  .;?  -  2  )  -  1 

JCOAIflJ  =  CWP'^X(0  0  0  0) 

10  cosriNUi-' 

>'  :  ■  "  .  2  ) 

IC'.  iBp  c  ir  fp'  CprcuUr  I»p,  Cap 
DO  20  I  r  S,  ti  ■'  *  2  )  -  1 

ZCGAT(I)  p  CMPLX, .  ,  1.0E.4  ) 

20  CONTISCE 

h  =  N  -  12  *  2  ) 

C?***?*?*  I  -  T  'j  o'  ctpecjnce  cc  ■  ng  modifications 

OPENvcM  ipNCh:'  J/  IlE-RUNZO/,  -JT’,srATUS=’Nt»’, 
1  rCiTW-  rORWAffEt)’) 


C  sycmetry  inputs 

r'l  =  ,-tCSM  0) 

SYMITP  ^  -’ROT’ 

Na^XS  ^  12 
IRHS-J  =  1 
.jNV  -  •  0 

C  Surf aca  i  Puts 

’or:  :>J  .  -i  /  Na.xs 

XiRti  r  0 

DO  i  1  r-  !  10 

CCCj  r  0.0 

1  CONTINUE 

DO  2  I  =  1,  3 

TRNSi'n  ■  0  0 

2  ;.TN-;NUi 

C  -  (he  Cone  - 

“TSEHNS  =13 

u.(  1)  1  .  CONVERT 

CC(  2)  =  2  .  CONVERT 

TRNS(l)  =  (0  00E-00)  *  lONVERT 

'rne;:;  ,,0  00.'. 00)  •  convert 

rRNS(3)  =  (0  00f-00)  .  convert 

XICAX  =  *3 

XSU.  =  -0  28100 

XSUU  =  0  26180 

XSVL  =  0 

XSVc  =  2  *  convert 

XNSU  --  2 

XNSV  -  2 

X!0K„'j  4 

XIGRD.  =  4 

XIRG  =  XIRG  r  1 

CALL  i.DSURRfXIRb,  XNSEQNS,  CC,  TRNS,  XIZAX, 

1  XS’X,  XSUU,  XSVL,  XSVU,  XNSU,  XNSV,  XIORDU,  XIOROV) 

C  -  -  -  —  c  I  rc" i ar  Imp  Cap  - 


XnSLSNS 

1  2 

‘I(  ' 

0 

.  convert 

1RNS(U 

^  ;0 

00E-08;  •  CONVERT 

;kni(:') 

-  (0 

00E.00)  •  CONVERT 

:rn:(3) 

(0 

00E.00)  f  CONVERT 

X 1  TAX 

.  .3 

XSU 

0 

.  CONVERT 

XS-fJ 

1 

.  CONVERT 

XuV. 

-  0 

20  10 

.'L 

i) 

/'■ns 
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XNSU  =  2 
XNSV  =  2 
XIOROO  =  4 
XIORDV  =  4 
XIRG  =  XIRG  .  1 

CAa  LDSURR(XIRG.  XNSBJHS,  CC,  TRNS,  XIZAX, 

1  XSUL,  XSUU,  XSVL,  XSVU.  XNSU,  XNSV,  XIORDU,  XIORDV) 

NSREG  =  XIRG 

C  CALL  PLOTCHIff(RUNID,  NBLKS,  SYMTYP,  1,  30.0,  60.0,  0.0) 

C  Define  Frequency  input. 

FREq  =  100 

C  Genente  surface  P  and  V  aatrices. 

CALL  SURkAT(FREH,  SYVfTYP,  NBLKS,  PMATX,  VMATX,  KDSIZE) 

C  Decoepose  Matrices 

CALL  DKOMM(SyMTYP,  NBLKS,  IRHSYM,  PMATX,  VMATX,  MDSIZE) 

CALL  I0SU8(NDVELS,  10,  va,  0) 

C  Block  I  1 

va(  1)  =  (1.0000,  0  0000) 

va(  2)  =  (1.0000  ,  0.0000) 

ve.(  3)  =  (l.NH,  0.0000) 

va(  4)  =  (1.0000  ,  0.00H) 

va(  s)  =  (0.0000  ,  0.0000} 

va(  6)  =  (0.0H0,  0.0000) 

va(  7)  ==  (0.0000  ,  0.0000) 

va.(  8)  =  (0.0000,  0.0000) 

CALL  iosue(NDvas,  1,  va,  nwdvec) 

C  Generate  Surface  Pressures. 

CAU  SURPRS(FREH,  SYMTYP,  NaKS,  IRHSYM,  PMATX,  VMATX,  MDSIZE) 

aOTYP  =  'FAR  ’ 

NUMTHP  X  19 
DO  25  I  X  1,  NUMTHP 

THTPHI(1,  I)  =  (I  -  1)  «  6 
THTPHI(2,  I)  =  0.0 
25  CONTINUE 

C  Calculate  FAR  Field  Matrices 

CALL  aDMAT(FREH,  SYMTYP,  FLOTYP,  NBLKS,  PMATX,  VMATX,  MDSIZE) 

IFFNRM  =  19 

C  Calculate  Far-Field  Pressures. 

CALL  FLDPRS(FRER,  aOHP,  NaKS,  IRHSYM,  PMATX,  VMATX,  MDSIZE) 

CALL  PRTOUT(FLDTYP,  0) 

STOP 

END 


Output 


of  Sample  Run  3 


FREOUENCYx  100.0  RH0=  1000.0  Cx  1500.0 
SYV0IEIRY  TYPEx  ROT  NO.  OF  aOCKSx  12 
FAR-FiaD  PRESSURB  (  PHI,THETA,FFP(a,IM,MAG,ANG) 


NORMALIZED  PATTERN,  SOURCE  LEVa,  DI  ) 


0.00  0.M  -.186832E*03  0.310326E<04  0  31088eE-04  9345  002 
0.00  5.M  -  .188244E<03  0.31B326E'>B4  0  310885E*04  93.43  0.02 
0  00  10.00  -  184483E'>03  0  310328E*04  0  310B75E>04  93.40  0.02 
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0.00 

15.00 

-.181564E.03 

0.310330E-04 

0  318860E-04 

03 . 35 

0.02 

0.00 

20.00 

- . 177508t-03 

0.310332E-04 

0  310839E.04 

93  27 

0.01 

0.00 

25  00 

- .  172346E-03 

0.310335F.04 

0  310813E-04 

93  18 

0.01 

0.00 

30  00 

-  166118E.03 

0  310338E.04 

0  310782E.04 

93  06 

0.01 

0.00 

35.00 

-  168870E-03 

0  31034lE*0‘i 

e.310/48E-a4 

“2  ZZ 

0  01 

0.00 

40.00 

-  150659E-03 

0  310344E*04 

0.3I0710E<04 

92  78 

0.01 

0.00 

45.00 

-.141646E-03 

0.31034eE'>04 

0.310669E>04 

92.61 

0.01 

0.00 

50.00 

-.131601E-03 

0  310348E-04 

6  310627E.04 

92  43 

0.01 

0.00 

55.00 

-.120899E-03 

0.310348E*04 

e.3105e3E>04 

92.23 

0.01 

0.00 

60.00 

-.109623E*03 

0  310347E-.04 

0.310640E<04 

92  02 

0.01 

0.00 

65.00 

-.975579E.02 

0  310344E*04 

0.310490E<04 

91  80 

0.01 

0.00 

70.00 

-.850956E^02 

0.310340E-04 

0.3104S3E*04 

91  57 

0.00 

0.03 

75.00 

'  .  •  4.i.n^U  1 

S.310333E-04 

0  310417E.F4 

91  33 

0.00 

0.00 

80.00 

-  590613E.02 

0.310325E-04 

0.310381E>04 

91  09 

0.00 

0.00 

85.00 

-.466876E-02 

0.310314E*04 

0 . 310348E>04 

90  84 

0.00 

0.00 

90.00 

-  322114E-02 

0.310301E«04 

0  310318E*04 

90.59 

0.00  • 

»  DIRECTION  FOR  PATTERN  NORMALIZATION 


LIMITATIONS 

The  folloKfing  existing  features  of  CHIEF  have  not  been  addressed  in  the 
context  of  a  partially  impedance  coated  radiator: 


•  Point  source  check.  There  seems  to  be  no  reason  that  this  option  cannot 
be  used  in  its  unmodified  form;  the  results  from  the  modified  CHIEF 
program  supply  the  actual  surface  velocities  all  over  the  radiator  which 
will  aid  one  in  assigning  relative  strengths  to  the  hypothetical  point 
sources . 


•  Radiator  in  an  infinite  rigid  baffle.  This  feature  has  not  been 
modified  for  the  case  of  an  impedance  coated  baffle;  furthermore,  the 
Green’s  fuu-jion  used  in  the  existing  formulation  of  this  problem  is 
applicable  only  to  the  case  of  an  infinite-plane,  rigid  baffle. 

•  Scattering  problems.  This  feature  has  not  been  modified  for  the  case  of 
a  partially  or  totally  impedance-coated  scatterer;  the  present 
development  in  CHIEF  is  applicable  only  to  rigid  scatterers. 

The  following  feature  of  CID  exists  for  a  partially  impedance  coated 
radiator ; 

•  Impedance  symmetry  must  be  the  same  as  body,  even  if  there  is  no 
velocity  symmetry. 
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APPENDIX  A  -  THEORETICAL  CONSIDERATION 


Al:  Basic  Equations  for  Impedance  Siirfaces 

The  CHIEF  program  develops  a  matrix  equation  for  the  unknown  surface 
pressure  P  of  the  form 


[A]P  =  [B]V.  (Al.l) 

Assuming  the  frequency  is  not  at  or  near  one  of  the  so-called  forbidden 
frequencies  where  the  coefficient  matrix  A  becomes  ill-conditioned,  or  in  fact 
its  determinant  vanishes,  then  one  can  solve  Eq.  (Al.l)  for  the  surface 
pressure  P.  Now  knowing  both  P  and  V  on  the  surface  of  the  radiator,  CHIEF 
proceeds  to  compute  the  radiated  pressure  field  at  any  point  external  to  the 
surface;  in  particular,  one  usually  wishes  to  know  the  far-field  pressure 
distribution. 


Now  consider  a  radiator  part  whose  surface  has  a  prescribed  normal 
velocity  while  the  remainder.  Region  I,  is  coated  with  the  Eirbitrary 
complex  locally  reacting  impedance  layer  Z.  Neither  of  these  regions  need  be 
a  single  contiguous  region;  i.e. ,  the  active  and  the  impedance-coated  areas 
could  be  a  number  of  noncontiguous  or,  in  fact,  even  disjoint  surfaces.  The 
subscripts  R  and  I  merely  indicate  a  convenient  sabdivision  of  the  problem 
into  the  two  different  boundary  condition  regions.  Accordingly,  matrix 
Eq.  (Al.l)  can  be  partitioned  as  follows: 


\  ^RR  ^RI 

) 

Pr' 

II 

®rr  \l 

[  Ajjj  Ajj  . 

[Pl  J 

®IR  ®II  • 

(A1.2) 


Note  that  the  impedance-coated  regions  are  relegated  to  the  bottom  positions 
in  the  matrix  equation,  a  choice  that  is  followed  throughout  this  discussion 
and  in  the  modifications  of  the  computer  code.  There  would  appear  to  Le  three 
sets  of  unknown  quantities  in  Eq.  (A1.2),  the  pressures  Pj^  and  Pj  as  well  as 

the  velocity  Vj-;  i.e.,  only  velocity  Vj^  has  been  specified.  However,  while 

neither  P^  nor  Vj  is  known,  their  ratio  defines,  or  is  defined  by,  the 

impedance  coating  on  Region  I  as 

Pj/Vj  =  -  Z.  (A1.3) 


The  minus  sign  in  this  equation  deserves  some  discussion.  Impedance  (specific 
acoustic  impedance)  is  defined  as  the  ratio  of  the  pressure  acting  on  the 
material  to  the  particle  velocity  flowing  into  the  material.  However,  the 
development  inherent  in  program  CHIEF  assumes  that  the  positive  direction  of 
the  velocity  is  outward  from  the  radiator  into  the  surrounding  medium.  Hence, 
to  use  the  velocity  designation  implicit  in  CHIEF  in  the  impedance  definition, 
the  minus  sign  is  required.  [The  authors  wish  to  acknowledge  particularly 
valuable  personal  communication  on  this  subject  with  Professor  Alan  Pierce, 
Oenartment  of  Mechanical  Engineering,  Pennsylvania  State  University. ] 
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Equation  (A1.3)  can  be  used  to  eliminate  either  of  the  unknowns  Pj  or  Vj 
from  Eq.  (A1.2)  as  a  function  of  the  other.  It  seems  more  appropriate  to 
eliminate  Py  since  the  opposite  choice  would  entail  a  division  by  Z  and  it  is 
recognized  that  Z  =  0  is  a  value  to  be  used  in  the  subsequent  calculations. 
Accordingly,  one  may  reformulate  the  matrix  equation  for  unknown  surface 
quantities  as 


^RR  ^~^RI^'®RI^  ^R  ^  ®RR  °  ^R 

^IR  1  i  Vj  J  I  0  J  [  0  . 


As  this  is  formally  the  same  as  Eq.  (Al.2),  one  may  proceed  to  have  the 
existing  routines  of  CHIEF  solve  this  modified  set  of  equations.  Once  P^  and 

Vj  are  known,  Pj  is  obtained  from  Eq.  (A1.3).  The  now  known  surface  pressure 

and  velocities  are  resorted  back  into  their  respective  array  positions,  and 
CHIEF  proceeds  from  this  point  with  the  calculation  of  the  pressure  at  points 
exterior  to  the  radiator  surface  in  the  usual  manner. 


In  this  development,  the  non  active  portions  of  the  surface  of  the 
vibrator  are  assumed  to  be  coated  with  a  true  locally  reacting  impedance 
material  which  creates  a  boundary  condition  at  those  surfaces  as  described  by 
Eq.  (Al.3).  This  impedance  can  have  a  value  (real,  imaginary,  or  complex) 
from  zero  (perfect  pressure  release  material)  to  the  largest  value  the 
computer  in  question  can  accept  (which  approximates  a  rigid  boundary) .  It 
should  be  noted,  in  particular,  that  specifying  this  impedance  value  to  equal 
the  characteristic  impedance  of  the  medium  does  not  mean  that  those  portions 
of  the  boundary  behave  as  if  they  were  acoustically  transparent;  i.e.,  because 
this  boundary  condition  does  not  guarantee  continuity  of  both  pressure  and 
normal  particle  velocity  at  the  interface  between  the  medium  and  the  impedance 
coating,  there  will  be  reflections  from  the  boundary  even  though  the  impedance 
value  matches  the  characteristic  impedance  of  the  medium,  unless  the  boundary 
is  effectively  planar  in  shape. 


A2:  Radiation  Impedance  Matrix 

The  radiation  impedance  matrix  (Z^jy^)  relates  the  forces  and  velocities  on 
the  various  subdivisions  as  follows 

?  ^  ^^RAD^I  (A2.1) 


However  when  part  of  a  radiator,  Region  R,  has  a  normal  velocity  v^  while  the 

remaining  K  surfaces,  Region  I,  are  coated  with  the  arbitrary  complex  locally 
reacting  impedance  layer,  Eq.  (A2.1)  can  be  written  as 


■  ?R  ]  t^RRi  tZRi]  I  f  IR  • 
Ff  [Zjr]  [Zjj]  Vj 
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The  subscripts  R  and  I  merely  indicate  a  convenient  subdivision  of  the  problem 
into  the  two  different  boundary  condition  regions. 

The  radiating  force  is  derived  from  Eq.  (A2.2)  as 

/r  "  *  (A2.3) 

In  order  to  determine  the  radiation  impedance  matrix  for  a  partially 
coated  radiator,  one  must  write  Eq.  (A2.3)  in  the  form  of  Eq.  (A2.1).  One 
approach  is  to  write  v^  as  a  function  of  Vj^.  By  substituting  the  equivalent 

form  of  Vj  in  Eq.  (A2.3),  a  product  of  Vj^  and  the  desired  radiation  matrix 

results . 

Therefore  the  first  step  is  to  reformulate  the  velocities  of  the  impedance 
surfaces  Vj.  Because  the  individual  terms  of  the  matrices  are  needed, 

summation  notation  will  be  used. 

For  the  impedance  surface,  the  force  takes  the  form 


F.  =  ^Z.  V  +  ^Z.  V  , 
1  ir  r  ip  p' 

r  P  ^ 


(A2.4) 


where  i  and  p  vary  from  (N  -  K  +  1)  to  N  and  r  varies  from  1  to  N  -  K.  Note 
that  v^  is  the  velocity  of  the  radiating  surface  and  v.  is  the  velocity  of  the 
impedance  layer .  ^ 

With  the  new  modified  CHIEF  program  a  complex  impedance  quantity  (Z.)  is 
specified  such  that  ^ 


F.  =  -Z.v. . 
1  11 


(A2.5) 


Setting  the  right  sides  of  Eqs.  (A2.4)  and  (A2.5)  equal  to  one  another,  we  get 


-Z.v.  -  V  +  21z.  V 

1  1  “  ir  r  “  ip  p 


(A2.6) 


Upon  bringing  the  term  on  the  left  side  of  Eq.  (A2.6)  under  the  summation  sign 
and  rearranging,  the  following  equation  results 


•Z  )V  =  Xz. 

»■!  r\*  r\  *'■ '  ^  1 


(A2.7) 


where  5^^  =  1  if  p  =  i  and  5^^  =  0  if  p  ^  i.  Equation  (A2.7)  relates  the 

velocities  on  the  impedance  surfaces  to  the  velocities  of  the  radiating 
surfaces  and  can  be  rewritten  in  matrix  notation  in  the  following  way 
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{-[Zii]  - 

0*4 

where  Z  is  the  prescribed  complex  impedance;  i.e.,  Z.  of  Eq.  (A2.5) 
Equation  (A2.8)  can  be  simplified  in  the  form 


^N-K+1  .  ® 


0 


•  Z. 


-I  “ 


or 


-I  “  ^t^IR^-R* 

Substituting  Vj  [Eq.  (A2.10)]  into  Eq.  (A2.3),  we  obtain 


or 

Equation  (A2.12)  can  br  abbreviated  as 

?R  "  ■^^RAD^YR’ 


where  Z^ 


is  the  desired  radiation  impedance  matrix. 


(A2.8) 


(A2.9) 


(A2.10) 


(A2.11) 


(A2.12) 


(A2.13) 
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APPENDIX  B  -  CID  SCREENS 


B1 :  Cone 

The  screen  shown  in  Fig.  B1 . 1  gives  the  geometry  description  for  a  conical 
surface.  Most  of  the  screen  is  self  explanatory  except  for  a  few  items.  The 
cone  is  rotated  about  the  X3  axis  with  the  point  of  the  cone  in  the  +X3 
direction.  In  order  to  flip  the  cone  about  the  Xl,  X2  plane,  the  maximum  cone 
height  must  be  negative.  The  minimum  cone  height  is  always  set  to  zero.  The 
cone  is  located  at  the  (0,0,0)  position  and  cein  be  relocated  by  changing  the 
local  origin. 


B2:  Impedance  Surfaces 

Figures  B2.1,  B2.2,  and  B2.3  relate  to  the  CHIEF  impedance  surface 
modification.  Figure  B2.1  looks  like  the  CID  screen  of  the  previous  revision 
with  the  exception  of  the  "IMPEDANCE  DBFS*  entry.  When  the  prompt  is  placed 
on  this  entry,  the  screen  shown  in  Fig  B2.2  appears.  The  user  can  pick  any 
surface  to  be  an  impedance  surface.  The  letter  "I"  will  appear  to  the  left  of 
each  surface  description  that  has  been  selected.  Furthermore,  if  there  are 
any  surfaces  after  the  one  selected,  they  will  automatically  be  impedance 
surfaces  and  the  letter  *1"  will  also  appear  to  the  left  of  these  surfaces. 
When  the  "EXIT*  entry  is  chosen  in  the  screen  shown  in  Fig  B2.2,  the  user’s 
screen  is  updated  to  the  screen  shown  in  Fig  B2.3. 

The  screen  shown  in  Fig.  B2.3  appears  for  every  surface  that  was  selected 
as  an  impedance  surface  from  the  previous  screen  (Fig.  B2.2).  The  user  has 
the  option  of  selecting  one  impedance  for  an  entire  region  or  selecting 
several  different  impedances.  If  different  impedances  are  selected  for  a 
region,  the  user  has  to  indicate  which  subdivision  is  the  first  to  have  an 
impedance  coating.  The  velocity  corresponding  to  the  impedance  subdivisions 
is  automatically  set  to  zero  in  the  driver  program  by  the  CID  program.  The 
CHIEF  program  calculates  the  velocities  of  the  impedance  surfaces  and  replaces 
the  zeros  in  the  velocity  vector  with  the  calculated  values. 
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Fig.  B2.1.  CID  screen  for  selecting  impedance  surfaces  option. 
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APPEND.TX  G  -  PROGRAM  LISTINGS 


Cl :  CURMAT 


The  subroutine  SURMAT  generates  the  matrices  A  and  B  of  Eq.  (A1.4).  In 
addition  to  the  new  variable  declarations 

C ( 1 MMU N / C DATING / ZCQ AT  NUMZ 
COMPLEX  ZC0AT(5OO), 

the  only  othe*  modification  is  to  call  the  new  subroutine  MTXSWZS. 


C  CHIPr  Vsrsior  3  a-P!J5-i991 


SUBRCUT:nE  GtRyUr'  RtR,S''MTYP,NBLKS,PMATX,V«AIX.UTXDIM) 

C  in'  ''  AND  V  SURFACE  MATRICES 

C  17  Sap  1990  -  i  u; 'iilaiiLO  coapifg  modifications  made  by  0U8 

PiARiMwCIlK  (MXSREG- S60'i 
PARAMETDl  'WXIPS-- ^0) 

PAPAWETER  (MXA;  o=.‘-''0) 

PARAMETER  (MXGAUS=b.; 

PaRamETEP  (WX?PTS:^512) 
c*SijrT[:p  jvX3(.vC:’n«) 
ov:  (VXf'Ap;,,  ^ 

CA"AM-rER  ;WXNFP--361) 

C  1  » 

PARAMETER  (MXrpS-  520) 

C  PA'.AME-ER  (Mt'PS  MAt'(MXARS*MXIPS,MXFFP,MXNFP)) 


Cl'A  ■■■':  -MA  x;d'  .;  IMi,  VMATX(M'rXDiM) 

CHARACTER»3  SYMT'rP 

r...«iMtN/C0N31/RH"  .C 
CaiiWON/PRTCnM/NUNPRr.NUNERR 

C  3  Jan  1991  -  addsJ  temporary  logical  unit  numbers  NTEMPVEt,  NTEMPl  -  NTE1IP3 

:  G Mi! ,  1  / 1 . AS  Cy  A  ^  J P :  3 ,  kD,  ■MX3 ,  NDVMXS ,  NDDECM ,  NOVELS ,  NDSPS , 

.  NDPMXF  N0VMXF,N0P«XN,N0VMXN,NDPSSP,NDB<PR,NDC01ilV, 

*  NDTEi(P,NDZROB,NoeATB,NTMPVEL,NTEllPl  NTEMP2.NTaiP3 
COM,!" N /SVALS/ SSREG ,  NSEQNS (MXSREO) , SUL (MXSREG) , SUU(UXSREG) , 

*  oVL (MXSREG) , SVU (MXSREG) , NSU(MXSREG) ,NSV (MXSREG) , 

.  t -XSREG) ,  TRNSS (3 , MXSREG) ,  IZAX(MXSREG) , 

»  IQROU(UXSREG) , lORDV(MXSREG) .NCCEHS 

CnVMON/lP 'S/NUMIPS, IPXS(3,MXIPS) 

Re*  , 

CGMMoN/JP''3/'t'St7,MXHPrS) 

CojJi(jN.'GAU3Ko/BKGAUS(UXaAU'S),R0uTU(MXGAUS),R00rV(MXGA','S), 

.  *ilJV(MXGAUS*l,2) 

C  JMWuN/rL.  -u  J /:  P9K !  (A  ,  MXFPS)  ,  FPC  (4 ,  MXpPS) 

CjMMuN/BGRK/VECKMXARS)  ,VEC2(UXBLXS),VEC3fMXaXS),VEC4(MX8L\S). 

*  VitS'MX:  "S)  ,VFCa(MXFPS) 

CCmi vEC  /  ,  ? u  3 ,  '/EC4 ,  VlCS,  VECfl 
O^'JENSIGN  FACTRSfUXARSI 
F^JWAlnNCE  (i'AciHS(l)  ,VEC1(1)) 
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CailttN/PRGVLS/NDIMPV , NUUARS , NUMSFP , NUUFFP , NUUNFP , NWDVEC 
COUMON/SURARS/ AREAS (MXARS) 

COWON/TAPREC/RECROClf] ,  IRECR0(3«] 

COIMON/TAPRCl/ARECRDCll} 

C***4*******  iipedance  coating  aodif ications  -  17  Sep  1990  ••«*««****»»« 
COMUON  /  COATINO  /  ZCOAT,  NUMZ 
COMPLEX  ZCOAT (MXARS) 

C«««««««««*««*«««««*«***««***«««*«««***»t««««««**«**4*«««**«««***t***»«« 


CHARACTER*4  ARECRD 

COMPLEX  JRHOWP 
C 

C  SURMAT  ERROR  CHECKS 

C 

IERR=0 

IF  (NSREG.GT.MXSREQ.OR.NSREG.LE.e)  THEN 
IERR=2 

WRITE(NUNERR,1M1)NSREQ 

1001  F0RMAT(/,’  BtROR  IN  SURMAT  -  NSREQ=',I6] 

END  IF 

IF  (NBLKS.QT.MXBLKS.OR.N8LKS.LE.0]  THEN 
lERRi:! 

IRITE(NUNERR,  1H2]NELKS 

1002  FORMAT (/,’  ERROR  IN  SURMAT  -  NBLKS=’.I6) 

END  IF 

IF  (SYMTYP.NE.’Rff’.AND.SYMTYP.NE.’ROT’)  THEN 
IERR=1 

IRITE(NUNERR, 1003)SYMTYP 

1003  F0RMAT(/,’  BWWR  IN  SURMAT  -  SYMTYP=  '.AS) 

END  IF 

IF  (NUMIPS.GT.MXIPS.OR.NUMIPS.LT.0)  THEN 
IERR=1 

IRITE(NUNERR,1004)NUMIPS 

1004  FORMAT (/,’  BIROR  IN  SURMAT  -  NUMIPS=’.I6) 

END  IF 

IF  (IERR.BR.2)  STOP 

NAR=0 

DO  1110  I=1,NSREG 
NAR=NAR«NSU(I)«NSV(I} 

IF  (NSERHS(I).GT.MCCEaS.OR.NSEQNS(I).LE.0)  THEN 
IERR=1 

0RITE(NUNERR,  1006)  I ,  NSEC|NS  (I) 

1008  FORMAT (/,’  ERROR  IN  SURMAT  -  REG|.NSERNS=’,2I6) 

END  IF 

IZABS=1ABS(IZAX(I)) 

IF  (IZABS.GT.3.OR,IZABS.EQ.0)  THEN 
IEkR=l 

fRITE(MUNERR, 1007) I, IZAX(I) 

1007  FORMAT (/,’  ERROR  IN  SURMAT  -  REG#, IZAX=' ,216) 

END  IF 

IF  (NSU(I)»NSV(I).LE,0)  THEN 
IEaR=l 

WRITE(NUNERR, 1008)1, NSU(I),NSV(I) 

1008  FORMAT(/,’  ERROR  IN  SURMAT  -  REG#,NSU,NSV-’,3I6) 

END  IF 

IORDRU=:IOROU(I) 

IORDRV=:IORDV(I) 

IF  (lORORU . GT . MXGAUS . OR . lORDRV . GT . MXGAUS  OR. lORDRUalORDRV 
•  .GT.MXHPTS.OR.IORDRU*IORDRV.LE.0)  THEN 

IERR=1 

IRITE(NUNERR, 1009) I , lORDRU, lORORV 

1009  FORMAT(/,’  ERROR  IN  SURMAT  -  REG#, lORDRU, IORORV=’ ,316) 
END  IF 

1110  CONTINUE 

IF  (NAR. GT. MXARS  OR. NAR.LE.0)  THEN 


SO 
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1EHR-=1 

IRHf '■SVNroR^  inilNAR 

nil  ^0RWAT:/,■  ERROR  IN  SURMAT  -  }0F  AREAS=',n0) 

(iNO  li' 

NFPirNAR^NKMIPS 

IF  (NFP  GT.UXFPS  OR.NFP  LE.01  THEN 
IERR-1 

*RITE(NUNbRR,  U12}NFP 

1112  FORVATC/,’  ERROR  IN  SURMAT  -  fOF  FIELD  POINTS^’  1101 
END  IF  '  ■' 

MTXT=NARtNFP 

irjNTxr.C'  N’XDIU  CR  WTXT.LE  0}  THEN 

*?.i:£fNUNLRR,ni3)MTXT 

SU.RMAT  -  ’,110,*  WORDS  NEEDED  FOR  P,V’) 

iF  (lERo.NE  a;  STOP 


FQURP’.-.i.e.pr 

RQMEGA=fl«OPI*FRE!} 

JRHCWP  C«PLX(a  0,ROWEGA<RHO)/FOURPI 

SMALLK=ROUEGA/C 

NOIMPV^UTXDIW 

'°“nuE‘*cJ^FiSi?^N5Sp^f  SUBDIVISION, AREA) 

'^’’“^^-''l^ARWr^FPRK?  PO^NT 


CALL  ICSIjB;NDClPTS,10,QPS  01 
ISUB=3 

DC  10  l-’,NSPEO 

V  ,  Roo  Tu,  iTuv  (1 , 1) ,  wkqaus) 
»TUV(MXSAJS.l,l)nOR[XJ(I) 

CALL  CIAoSS(!OROV(I},ROOTV,ITUV(1,2),WKGAUS) 

I TUV  (NXO  A'.iS  •  1 , 2 1  n  C  RDV  ( I )  ’ 

0£LS!J=(3Lt'n)-SU.(I))/NSUa) 

Da..SV=(SVU{I}-SVL(I))/NSV(I) 

HDasu=a  6»Dasu 
HOEi,‘;;’/=0  s.cFi  ?v 
FAcn--HDas:j»Hoasv 
UAVE-HOaSU*SUL(I) 

DO  8  J=1,NSU(I) 

VAVF=HDFLSV.SVL(I1 
DO  6  K-1,NSV(I) 

ISUB^ISuBn 

FACTRS(ISU9)=FACTi 

CAL;,  Ap'or.N;uAVE,VAV£,HDaSU,HOasV,NSEHNS(I)  CCSfl  1) 

•  U,  izAX(i),  iordu(I),iorovF'),ro6tu'rootv,qps) 

call  ;C'iUH(>(DQPIS,  l,fTUV,2*(k(XGAUS.l)) 

cal'  n3L'S;:Nfn-'’rE, i,qps,7*iordu(i)»iorov(i)) 

CAl'.,  n-;.‘iVD:.lAVt..'AVF,NSEUNS(I),CCS(l,I),TRNSS(l  1)  RAX(') 

•  '  PCK 1(1,1  SUB) , FPBKl (2, ISUB) , FPBKl (3, ISUB) . Dl, 02,03; DDJ, 0) 

VAV  '  Df 
LAVlJ^UAv;  .Df.-LSJ 
CONI  INUP- 

NJMAPS-ISDB 

NUMEp-i-iSUB 

NKOVEr.?.^! 


'  Nfc  P,  iHlN 
DO  /0  1  ;,N'.,A(('’S 
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isue=isuB*i 
DO  16  1=1,3 

16  FPeKl(I,ISU8)=IPXS(I,J) 

29  CONTINUE 

END  IF 

NUrPS=NUUFPS-NU«PS 

NUIlSFP=NUIfPS 

irXTOT=NUMARS»NUMFPS 

C  GENERATION  OF  SURFACE  IIATRICES 

C  P  MATRICES  TO  FILE/UHIT  NOPMXS 

C  V  MATRICES  TO  FILE/UNIT  NDVMXS 


INTTYP=1 

CALL  IOSUe(NOPMXS,ie,PMATX,l) 
CAa  I0SUB(NDVMXS,ie,VMATX.6) 


C 


c 


DO  60  IBLK=1,NBLKS 

««««*•««**«*«*•  SYMMETRY  CONSIDERATIONS  ***************** 


IF(SY¥TYP.Bl. ’Rff ’)  THEN 
CALL  FPRffC(IBLK,NUrPS. 
END  IF 


INTTYP.FPBKl.FPC) 


IF  (SYMTYP.EH. ’ROT’)  THEN 
<iAMMA=T¥0PI/N8LKS 

CALL  FPROTC(IBLK, GAMMA, NJf PS, INTTYP,FPBKl,FPCj 
END  IF 


«»*#»*»*»«»«»««««»«*««»»*»*««»»«««<****»»»***************** 


CALL  MATCENfNUMARS, NUMFPS,FACTRS,  INTTYP,WTUV,MXGAUS, 
NOQPTS  PC ,  SMALLX ,  PMATX ,  VMATX) 


IFdBLK  EQ.IT  CAU  SLFINT(PMATX,VMATX,NUMARS,NUMFPS.SMALLK,INTTYP, 
.  IORDU.IOROV.ROOTU, ROOTY, WTUV,*KGAUS,MXGAUS, 

,  NSREG,NSEQNS,CCS,TRNSS,IZAX,SU.,SUU,SVL,SVU, 

,  NSU,NSV,FPC, AREAS) 


DO  30  I=1,MTXT0T 
WIATX(I)=VMATX(I)»JRHOfP 
30  PMATX(I)=-PMATX(I)/FOURPI 

IF(IBLK.EH.l)  THEN 
DO  40  I=l,NUMARS 
K=(I-l)«NUiFPS*I 
40  PMATX(K)=PMATN(K)*0.6 
END  IF 


C«««**«**««*  iipedince  coiling  lodif ications  -  17  Sep  1990  «♦»»«»♦»♦♦♦»• 
CALL  MTXSf2S(PMATX,  VMATX,  NUMARS,  NUMZ,  NU4IPS,  ZCOAT) 


CALL  lOSUe (NOPMXS, l,PMATX,2iMTXT0T) 
CALL  lOSUB (NDVMXS,!, VMATX, 2»MTXT0T) 

60  CONTINUE 

RKRO(l)=FREll 

RECRD(2)=RH0 

RECR0(3)*C 

IRECRO(l)=Na.KS 

IRKR0(2)»MTXDIM 

IRECRO(3)*NUMARS 

IRBCRD(4)=NUMSFP 

ARECRD(1)=SYMTYP 

RETURN 

END 
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G2:  SURPRS 

The  subroutine  SURPRS  solves  matrix  equation  (A1.4)  for  unknown 
quantities.  In  addition  to  the  new  variable  declarations 

GQMMON/COATING/ZCOAT ,NUMZ 
COMPLEX  ZC0AT(5OO), 

the  subroutine  SURPRS  computes  the  remaining  surface. pressures  from  Eq.  (A1.3) 
and  then  places  the  values  of  the  velocities  on  the  impedance-coated  regions 
back  into  the  array  VEL  (into  the  last  NUMZ  positions) .  The  [B]  matrix  of  Eq. 
(Al.2)  is  not  actually  modified  to  create^ the  indicated  right  hand  side 
coefficent  matri’'  of  Eci.  (A1.4);  rather  the  indicated  elements  Bgj  and  Bjj  are 

left  intact.  This  can  be  corrected,  however,  by  inputting  zero  values  into 
corresponding  positions  in  the  velocity  vector  on  the  right  hand  side  of  Eq. 
(A1.4).  It  is  critically  necessary,  therefore,  that  one  initially  specify 
zero  values  for  what  would  appear  to  be  the  velocities  on  the  impedance  coated 

regions . 


c .  . . . . - 

c 

C  CHIEF  Vetsion  3,1«  8-FEB-1991 
C 

C . .  . - . 


SURROtT I NE  SI IRPRS (FREQ , SYKTYP, NBLKS , IRHSYM, CUTX. CMTEMP, MTXDIM) 

C  BUILDS  RIGHT  HAND  SIDE, 

C  SGt.VFJ  EQUATIONS,  AND  COMBINES  SOLUTIONS  USING  SYMMETRY  lEIOHTS 

C  TO  OBTAIN  SURFACE  PRESSURES 

C  17  Sep  1990  -  lepedance  coating  eodif ications  aade  by  DMB 

C  7  Jan  1990  -  btg  fix  SURPRS  lade  by  PJK 

PARAMETER  (MXARS=500) 

PARAMETER  (MX8LKS=100) 

PARAMETER  (MXIPS=20) 

PARAMETER  (MXFFP=361) 

PARAMETER  (MXNFP=361) 

C  •*»« 

PARAMETER  (MXFPS=52e) 

C  PARAMETER  (MXFPS=MAX0(MXARS*MXIPS,MXFFP,MXNFP)) 

C 

COMPLEX  CMTX(MTXDIM),CMTCMP(MTXDIM) 

CHARACTER*3  SYMTYP 


tOMMCN/CCNST/RHO.C 

COMMON/PRICOM/NUNPRT.NUNERR 


C  3  Jan  1991  -  added  teeporary  logical  unit  nuebers  TEMPVB.,  TBJPl  -  TEMP3. 


COkWO  N/  NDASG/  NDQPTS ,  NDPMXS ,  NDVMXS ,  NOOECM ,  NDVaS ,  NDSPS , 

.  NDPMXF , NDVMXF , NDPMXN , NOVMXN , NOPSSP, NOEXPR, NDCOMV, 

«  NDTEMP , NOZRDB , NDPAT8 , NTMPVEL , NTEMP 1 , NTEMP2 , NT  EMP3 

COMMCN/PRGVLS/NDIMPV, NUMARS , NUMSFP, NUIFFP, NUMNFP, NWDVEC 
COMMON/PLiINP/AINC,THTINC,PHIINC,ISCATR 
COMMON/BAFFLE/ INFRIG 
CCMMON/MODAL/IMQOAL,  lUSFO 

COMMON/WORK/VEC1(MXARS),VEC2(MXBLKS),VEC3(MX0LKS),VEC4(UXBLKS), 
*  VECS(MXFPS),VEC6(MXFPS) 

COMPLEX  VEC2,VFC3,VK4,VEC$,VEC6 

COMPLEX  RC0TUN(MXGLK3),  SYMfTS(MX8LKS) ,  PRSf TS (MXBLKS) 

COMPLEX  X(MXFPS),  CRHS(MXFPS),  CSOLN(MXFPS) 

EQUIVALENCE  (ROOTUN(l) , VEC2(1)) 
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ERUIVALBKE  (SYKITSCl) ,VEC3(1)) 

EHUIVALDICE  (PRS»TS(1),VEC4(1)) 

EIIUIVAL9ICE  (xa),VBCS(l)) 

BJUIVALBICE  (CS0LN(l),ve:8(l)) 

SJUIVALENCE  (CRHS(l) .VEC6(1)) 

COHIQM/aDSVEC/TVKTCliIXARS)  ,0{MXARS) ,  IPIVTR(MXFPS) 

COMPLEX  TVKT 

comoM/vasps/va.(uxARS)  ,sp(mxars) 

COMPLEX  VEL.SP 

COMMOM/POISL/POIER . DIRINO, SRCLVL 
COIMN/SURARS/ AREAS  (MXARS) 

COIilOM/EXTCOM/EXTPRS(MXFPS) ,  lEXTFO 
COMPLEX  EXTPRS 

COMMON/NBPRTC/IRHSPT, NARSPT , NPTBLK , FRQPT 

COMMON/NBPRTS/SYITPT 

CHARACTER*3  SYMTPT 

C************  lipedanca  coating  lodif ications  -  17  Sep  1993  <«««••*««««* 
COMMON  /COATINC/  ZCOAT,  NUM2 
COMPLEX  ZCOAT(MXARS] 


IRHSPT=IRHSYM 

FROPT=FREH 

SYirPT=SYMTYP 


C 

C  SURPRS  ERROR  CHECKS 
C 

IERR=f 

IF  (NaXS.GT.MXBUS.OR.NBLKS.LE.B)  THEN 
IERR=«1 

IRITE(NUNERR,1U2)NBLKS 

1002  F0RMAT(/.’  ERROR  IN  SURPRS  -  NBLKSs'.IB) 

END  IF 

IF  (SYMTYP.NE. 'Rff’.AND.SYMTYP.NE. ’ROT’)  THEN 
lERRrl 

fRITE(NUNERR, ief3)SYMTYP 

1003  F0RMAT(/,’  ERROR  IN  SURPRS  -  SYMTYP=  ’,A3) 
END  IF 

IF  (IERR.NE.0)  STOP 

mDMTX=2«NUMSFP*NUMARS 

NWDVEC=2«NUMARS 

CALL  IOSUB(NDOECM,lB,CMTX,f) 

CAa  IOSUB(NOTaP,10,CSOLN,0} 

CALL  IOSUB(NOCOMV,1B,CMTX,S) 

IF(IRHSYM.BI.I)  THEN 
NBS0LN=N8LKS 
DIVSR=N8LKS 
B-SE 
NBS0LN=1 
DIVSR=1.0 
END  IF 

DO  m  IS0LN=1,N8S0LN 


C  •••«••«•***«**«  SYMMETRY  CONSIDERATIONS 


IF(SYMTYP.EH.’Rff’)  CAa  RffSWT(ISOLN,NBLKS,SYMITS) 

IF(SYMTYP.EH,’ROT’)  THEN 

IF(ISOLN.Eq.l)  THEN 
PI=ACOS(-1.0) 
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Ti|OPI=2.0*PI 
ARG=0 . 0 

DELARG^TWOPI/NBLKS 
DO  10  I=1,NBLKS 

ROOTUN(I)=ClilPLX(COS(ARG)  SIN(ARG)) 

10  ARG=ARG-.DB-ARG 
END  IF 

CALL  ROTSf T (ISOLN , NBLKS , ROOTUN , SY«*TS) 

END  IF 


IF(ISCATR  ER.0)THEN 

CALL  IOSUB(;NOCaMV,2,CMTX,N»DMTX) 

CALL  COMBMX (X , VcL , NUMARS , 1 , NOVELS , NWDVEC , NBSOLN . SYUWTS) 

CALL  MX«LTV(CUTX,X,CRHS,NUIilSFP. NUMARS) 

ELSE 

DO  16  I=1,NUMSFP 
16  CRHS(I)=(0.0,f  0) 

ENDIF 

IF(1EXTFG  NE  0)  THEN 

CALL  COUBMX (X , EXTPRS , NUMSFP, 1 , NDEXPR, 2»NUUSFP, NBSOLN, SYMWTS) 
DO  18  1=1, NUMSFP 
18  CRHS(I)=CRHS(I).X(I) 

END  IF 

DO  20  1=1, NUMSFP 
20  CRHS(I)=CRHS(I)/DIVSR 

IF(INFRIG  E5,0)  TH  N 

CAU  I0SUB(NDDECM,2,CMTX,NWDMTX) 

CALL  IOSUB(NDDECM,2,TVECT,3*MXARS«MXFPS) 

CALL  QDS (2 . CMTX , NUMSFP, NUMARS , TVECT , 8, IPI VTR, CRHS . CSOLN) 
asE 

DO  25  1=1, NUMARS 
25  cr'”‘';:)-:.0»CRHS(i) 

END  IF 

CAa  IOSUB(NDTaiP,l, CSOLN, NIDVEC) 

100  CONTINUE 


CALL  IOSUB(NDSPS,10,SP,0) 

CALL  IOSUB(NDVaS,10,Va,0) 

C  this  rewind  idded  as  part  of  7  Jan  19'H  bug  fix  «*««««*«=«**** 

CAa  iosuB(NTMPva,i0,va,0) 


POWER=0.0 

DO  200  IS0LN=1, NBSOLN 


C  SYMMETRY  CONSIDERATIONS  ♦♦«»**»*♦•♦♦*♦*♦**♦ 


IF(SYMTYP.Ea. 'REF')  THEN 
CALL  REFSWT (ISOLN, NBLKS, SYM0TS) 
DO  6  1=1,  NBLKS 
5  FRS»TS(I)=SYM»TS(I) 

ENDIF 

IF(SYMTYP  EU  'ROT')  THEN 
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CALL  ROTSIT(ISOLN,N0LKS,ROOTUN,SYUITS) 
PRSWTS(l)=SYyiTS(l) 

IF(NBLKS.NE.l)  THEN 
00  12  I=2,NBLKS 

12  PRS»TS(1)=SYMITS(NBLKS-I.2) 

END  IF 
END  IF 


C 


>«« 


CALL  COIAIX  (SP ,  CSOLN .  NUMARS ,  1 ,  NOTBilP .  NIOVEC ,  NBSQLM ,  PRSNTS) 
IF  (ISCATR  .EH.  OjCALL  I0SU9(NDVa.S,2.Va..NWDVEC) 


C************  iip«(ltnc«  coating  aodif ications  -  17  Sop  1998  ************ 
DO  411  I  =  NUMARS  -  NUUZ  «  1,  NUMARS 
va.(I)  =  SP(I) 

SP(I)  =  -SP(I)  *  ZCOAT(I) 

481  CONTINUE 


C*********************************************************************** 
C«****«**«««Bug  fix  of  problaa  with  Velocities  starts  here  -  7  Jan  1991 


C  Modified  when  the  velocity  ws:  stored  in  files  for  the  iepedance 
C  calculations. 

CAU  IOSUB(NTMPVa,l,Va..N«DVEC) 

CALL  IQSUB(NDSPS,1,SP,NWDVEC] 

IF(ISCATR.NE.8.]G0  to  2M 
DO  lEi  1=1, NUMARS 

1 58  POf  ER=P0f  ER-REAL  (SP  (I)  »C0N  JQ  (Va  (I) ) )  •  AREAS  (I) 

288  CONTINUE 

CAa  IOSUB(NTMPVa,18,Va,0) 

CALL  IOSUB(NOVaS.18.Va,8) 

C  IRITE  MATICES  CONTAINED  IN  NTMPVa  TO  NOVaS 

DO  225  I=1,NBS0LN 

CALL  I0SUB(NTMPVa,2,Va,NIDVEC) 

CALL  lOSUB (Novas, l,Va,NIOVEC) 

225  CONTINUE 

••aaaaaaaaeeeeata  ^nd  of  7  Jan  1991  aods  ************************************ 

IF(IM0DAL.NE.8)THEN 
CALL  I0SUB(NDTQP,18,SP,8) 

DO  388  IS0LN=1,NBS0LN 

CALL  IOSUB(NDTaiP,2,SP,MWOVK) 

388  CALL  I0SUB(NDSPS,1,SP,N1DVEC) 

ENDIF 

IF(IRHSYM.EQ.  l)P0IER=NaKS«P01ER 
P0fER=PafER/1888. 

IEXTFG=8 

RETURN 

END 
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C3;  MTXSWZS 

The  subroutine  MTXSWZS  creates  the  modified  coefficient  matrices  of  Eq. 
(A1.4)  from  those  of  Eq.  (A1.2).  Note  that  this  formulation  allows  for  NUMIPS 
number  of  interior  points,  and  consequently  the  matrices  A  and  B  in  this 
subroutine  are  not  necessarily  square.  Further  note  that  the  B  matrix  of  Eq. 
(A1.2)  is  not  modified  to  create  the  exact  right  hand  coefficient  matrix  of 
Eq.  (A1.4).  This  is  legitimate  provided  one  does  not  input  zero  values  for 
the  corresponding  positions  in  the  velocity  vector  on  the  right-hand  side  of 
Eq.  (A1.4).  Also  note  that  the  impedance-coated  regions  are  always  associated 
with  the  latter  storage  positions  in  matrices  A  and  B  and  in  the  arrays  SP  and 
VEL  regardless  of  where  the  regions  actually  are  on  the  surface  of  the 
radiator . 


c -  - 

c 

C  CHIEF  Version  3  10  8-FEB-1991 
C 

C  MTXSWZS. FOR  «as  added  to  CHIEF  in  Septeiber  1990  by  W.  Thoapson 
C - - 

SUBROUTINE  MTXSWZS(A,  B,  N,  NZ,  NUUIPS,  ZC) 

(_ttt*tttttHtt*t%*tttt*t**tt*****f*^******t**********f  ****************** 

C  required  in  CHIEF  for  iapedance  coating  cases. 

C  NOTE:  aatrix  B  in  eqn  AP^^BV  is  not  aodified. 

C»«««««  «««««*««••«««««•«««<«*«*«««*•*«««•«***«**«««*••«««**•*•••«*••••*« 

COMPLEX  A(N. NUMIPS,  N) ,  B(N«NUUIPS,  N) ,  ZC(600) 

NN  =  N  •  NUMIPS 

00  20  I  =  I,  NN 

00  10  J  =  N  -  NZ  ♦  1,  N 

A(I,  J)  =  -A(I,  J)  »  ZC(J)  -  B(I.  J) 

10  CONTINUE 

20  CONTINUE 

RETURN 

END 


C4:  ZoAiZaJ. 

The  subroutine  ZRADMI  is  called  from  the  CHIEF  driver  and  is  used  to 
calculate  the  radiation  matrix  only  when  impedance  surfaces  are  present.  When 
impedance  surfaces  are  present,  both  ZRADMX  and  ZRADMI  are  called. 

c - - - - — 

C 

C  CHIEF  Version  3.10  8-FEB-1991 
C 

C - - - . - . - - - - - 

SUBROUTINE  ZRA0MI(ZMTX,  NUMARS) 

C  ZRADMIO  IS  a  subroutine  that  corrects  the  radiation  iapedance  aatria 
C  when  an  arbitrary  iapedance  layer  exists  on  the  surface. 

C 

C  Subroutines  used:  MATRIXINV() 

PARAMETER  (MXARS  =  500) 

C  iapedance  layer  coaaons.  . 

COMPLEX  ZCOAT(60I) 
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CQUUDN  /  COATING  /  ZCQAT,  NUMZ 

COMPLEX  ZRR(UXAftS,  MXARS) 

COMPLEX  ZIR (MXARS,  MXARS) 
COMPLEX  ZRI (MXARS,  MXARS) 
COMPLEX  ZIl (MXARS,  MXARS) 

COMPLEX  INVZII(MXARS,  MXARS) 
COMPLEX  TBP2 (MXARS,  MXARS) 
COMPLEX  TaiP3 (MXARS,  MXARS) 

REAL  TaM>4(MXARS  *  2,  MXARS  «  2) 
REAL  TBr6(MXARS  «  2,  MXARS  «  2) 
INTEGER  ITa«>(MXAfiS) 

INTEGER  IFIRST,  lUST 
INTEGER  RFIRST,  RLAST 

COMPLEX  ZMTX 

DIMBISION  ZMTX(NUMARS,  NUHARS) 


C  The  radiation  iepedance  eatrix  is  coeputed  assuaing  there  are  no 
C  iapedance  surfaces.  This  aatrix  ZMTX  has  diaensions  NUMARS  x  NUHARS. 
C  ZMTX  is  partitioned  into  4  aatrices,  ZRR,  ZRI,  ZIR,  and  ZII  ahose 
C  sizes  depend  upon  the  nuaber  of  iapedance  surfaces. 

RFIRST  =  1 

RLAST  NUMARS  -  NUMZ 

IFIRST  =  NUMARS  -  NUMZ  «  1 
ILAST  =  NUHARS 

DO  3f  J  >  RFIRST,  RLAST 
DO  U  K  =  RFIRST,  RLAST 
ZRR(K,  J)  =  ZMTX(K,  J) 

16  CONTINUE 

DO  26  K  =  IFIRST,  ILAST 

ZIR(K  -  RLAST,  J)  =  ZMTX(K,  J) 

26  CONTINUE 

36  CONTINUE 

DO  86  J  =  IFIRST,  ILAST 
DO  46  K  =  RFIRST,  RLAST 

ZRI(K,  J  -  RLAST)  =  ZMTX(K,  J) 

46  CONTINUE 

DO  66  K  =  IFIRST,  ILAST 

ZII(K  -  RLAST,  J  -  RLAST)  =  ZMTX(K,J) 

56  CONTINUE 

66  CONTINUE 

C  The  iapedance  coatings  are  added  to  the  diagonal  eleaents  of  the 
C  ZII  aatrix. 

DO  76  J  =  1,  NUMZ 

ZII(J,  J)  =  ZII(J,  J)  -  ZC0AT(J  ♦  RLAST) 

76  CONTINUE 


In  order  to  invert  the  coaplex  aatrix  ZII,  we  partitioned  it  into 
a  real  aatrix,  TBiP4,  that  is  4  tiaes  larger  as  follows: 


I  Real  ZII  -laaginary  ZII  | 

(  I 

j  laaginary  ZII  Real  ZII  j 


and  inverted  this  larger  real  aatrix. 

DO  96  J  -  1,  NUMZ 
DO  86  K  =  1,  NUMZ 

Ta6»4(  K,  J)  =  REAL(ZII(K,  J)) 

Ta6'4(NUMZ  ♦  K,  NUMZ  ♦  J)  =  REAL(ZII(K,  J)) 

TB6>4(NUMZ  ♦  K,  J)  *  AIMAG(ZII(K,  J)) 

TBi>4(  K,  NUMZ  «  J)  =  -AIMAG(ZII(K,  J)) 

86  CONTINUE 
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90  CONTINUE 

CALL  UATRIXINV(TEMP4,  TB0P5,  ITEUP,  NUMZ*2,  MXARS.2) 


C  After  inversion,  the  result  is  still  partitioned  in  the  sale  way 
C  The  real  and  iaaginary  parts  of  the  inverse  are  now  recoabined  into 
C  the  coaplex  array  INVZII. 

C 

C 

C  I  «*  Real  INVZII  •*  -Iaaginary  INVZII  | 

C  1  1 

C  I  «♦  Iaaginary  INVZII  »«  Real  INVZII  j 

C 

DO  110  J  =  1,  NUMZ 
DO  100  K  =  1,  NUMZ 

INVZII  (K,  J)  =  CUPU(Ta(P6(K,  J) ,  TEMP6(NUUZ  ♦  K.  J)) 

100  CONTINUE 
110  CONTINUE 

C  ZRI  is  aultiplied  by  the  inverse  of  ZII  and  put  into  the  TEMP2  array 

DO  140  J  1,  NUMZ 

DO  130  K  =  RFIRST,  RLAST 
TaiP2(K,  J)  =  0.0 
DO  120  L  =  1,  NUMZ 

TaiP2(K,  J)  =  TEMP2(K,  J)  -  ZRI(K,  L)  *  INVZII(L,  J) 

120  CONTINUE 

130  CONTINUE 

140  CONTINUE 

C  TeBp2  is  aultiplied  by  ZIR  and  put  into  the  array  TEMP3 

DO  170  J  =  RFIRST,  RLAST 
DO  180  K  =  RFIRST,  RLAST 
TEMP3(X,  J)  =  0  0 
DO  150  L  =  1, 

TEMP3(K,  J)  =  TaP3(K,  J)  ♦  Ta(P2(K,  L)  •  ZIR(L,  J) 

150  CONTINUE 

160  CONTINUE 

170  CONTINUE 

C  ZMTX  is  now  the  correct  radiation  iapedance  aatrix  and  is  passed 
C  back  to  the  CHIEF  driver.  (Note  the  size  of  ZMTX  is  larger  than 
C  the  size  actually  needed. 


00  190  J  =  RFIRST,  RLAST 
DO  180  K  =  RFIRST,  RLAST 

ZMTX(K,  J)  =  ZRR(K,  J)  -  TaiP3(K,  J) 
180  CONTINUE 
190  CONTINUE 

RCTURN 

END 


05:  MATRIUNV 


The  subroutine  MATRIXINV  is  called  from  ZRADMI.  This  subroutine  inverts  a 
matrix  and  calls  the  subroutine  DECOMP  and  SOLVER. 


c - - 

c 

c  CHIff  Version  3.10  8-FEB-1991 
C 

C - 


SUBROUTINE  MATRIXINV(MAT,  INV,  INDEX,  N,  S) 


C  Inverts  a  aatrix  (MAT) 
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C  Subroutines  used:  DECQ«P(),  S0LVER() 

INT6GQR  N  !  defines  size  of  trruy  to  be  inverted  (N  x  N) 

INTEGB^  S  'size  of  storage  area  (S  x  S)  (aay  be  larger  than  array) 

INTEGER  INDEX  !  used  in  bookeeping. 

REAL  MAT  !  array  to  be  inverted  (destroyed  ehen  done) 

REAL  INV  !  holds  inverted  array  when  done. 

INTEGER  D  !  «!  or  -1  depending  on  even  or  odd  nuaber  of  row  changes 

DIMENSION  MAT (S,  S),  INV(S,  S),  INDfX(S) 

DO  2«  I  =  1,  N 
DO  II  J  =  1,  N 
INV(I.  J)  =  0.0 
10  CONTINUE 

INV(I.  I)  =  1.1 
20  CONTINUE 

CALL  DECOMP(MAT.  N,  S,  INDEX.  0) 

DO  30  J  =  1  N 

CAU  S0LVER(MAT,  N,  S,  INDEX,  INV(1.  J)) 

30  CONTINUE 

RETURN 

END 


C6:  DBCQMP 


The  subroutine  DCOkP  performs  a  matrix  decomposition. 


c . . 

c 

c  CHIff  Version  3.10  0-FEB-1991 
C 

C . 


SUBROUTINE  0ECaHP(A,  N,  NP,  INDX,  D) 


C  Matrix  Oecoaposition  routine. 


PARAMETER  (NMAX  =  100,  TINY  =  1.0E  -  20) 

DIMENSION  A(NP,  NP),  INDX(N),  W(NMAX) 

D  =  1.0 

00  12  I  =  1,  N 
AAMAX  =  0.0 
DO  11  J  =  1,  N 

IF  (ABS(A(I,  J))  GT.  AAMAX)  AAMAX  =  A8S(A(I.  J)) 


11 

CONTINUE 

IF  (AAMAX  .ER.  0,0)  PAUSE  'Singul 

ar  aatrix 

W(I) 

=  1.0  /  AAMAX 

12 

CONTINUE 

00  19  J  : 

=  1,  N 

IF  (J 

GT.  1)  THEN 

DO 

14  I  =  1,  J  -  1 

SUM  =  A(I,  J) 

IF  (I  .GT.  1)  THEN 

DO  13  K  =  1,  I  -  1 

SUM  =  SUM  -  A(I,  K)  « 

A(K,  J) 

13 

CONTINUE 

A(I,  J)  =  SUM 

ENDIF 

14 

CONTINUE 

ENDIF 
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AAUAX  =  00 


DO  16  I  =  J,  N 
SUM  =  A(I,  J) 

IF  (J  .GT.  1)  THEN 
DO  15  K  =  1,  J  -  1 

SUM  =  SUM  -  A(I,  K)  «  A(K,  J) 
15  CONTINUE 

A(I,  J)  =  SUM 
ENOIF 

DUM  =  VV(I)  «  ABS(SUM) 

IF  (DUM  GE,  AAMAX)  THEN 
IMAX  ^  I 
AAMAX  =  DUM 
ENOIF 

18  CONTINUE 

IF  (J  NE.  IMAX)  THEN 
DO  17  K  =  1.  N 
DUM  =  A(IMAX,  K) 

A(IMAX,  K)  =  A(J,  K) 

A(J,  K)  =  DUM 
17  CONTINUE 

D  =  -D 

W(IMAX)  =  W(J) 

ENDIF 

INDX(J)  =  IMAX 
IF(J  .NE.  N)  THEN 

IF(A(J,  J)  EH.  0.0)  A(J,  J)  =  TINY 
DUM  =  1.0  /  A(J,  J) 

DO  18  I  =  J  •  1,  N 

A(I,  J)  =  A(I,  J)  *  DUM 


18 

CONTINUE 

ENDIF 

19 

CONTINUE 

IF(A(N,  N)  .ER.  0.0)  A(N,  N)  =  TINY 

RETURN 

END 


C7:  SOLVER 

The  subroutine  solves  a  matrix  equation . 


c - - - 

c 

C  CHIEF  Version  3  10  0-FEB-199I 
C 

C— - . - . - 


SUBROUTINE  S0LVER(A,  N,  NP,  INDX,  B) 


C  This  routine  solves  the  set  of  N  linsir  equations  A  .  X  =  B 
DIMENSION  A(NP,  NP) ,  INDX(N),  B(N) 

II  =  0 

00  12  I  =  1,  N 
LL  =  INDX(I) 

SUM  =  B(LL) 

B(LL)  =  B(I) 

IF  (11  NE.  0)  THEN 
DO  11  J  =  II,  I  -  1 

SUM  =  SUM  -  A(I,  J)  «  B(J) 

11  CONTINUE 


81 


r>  o  r>  r>  o 


Thompson,  Dubus 
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aSE  IF  (suy  .NE.  0  0)  THEN 
II  =  I 
ENOIF 

B(I)  =  SUM 

12  CONTINUE 

DO  14  I  =  N,  1,  -1 
SUM  =  B(I) 

IF(I  .LT.  N)  THEN 
DO  13  J  =  I  ♦  1,  N 

SUM  =  SUM  -  A(I,  J)  .  B(J) 

13  CONTINUE 
ENOIF 

B(I)  =  SUM  /  A(I.  I) 

14  CONTINUE 

RETURN 

END 


C8:  mTCM 


CHia  Va'aion  3.10  a-FEB-1991 


SUBROUTINE  INITCM 

PARAMETER  (MXSREG=500) 

COtAUN/CONST/RHO.C 
COMMON/PRTCOM/NUNPRT,  .s'JNERR 

C  3  Jan  1991  -  add«d  tsaporary  logical  unit  nuibors  NTBA’va,  NTBIPl  -  NTE^P3. 

COIMON/N  jASQ/NDdPTS ,  NOPMXS ,  NDVMXS ,  NDDECM ,  NDVaS ,  NKPS , 

*  NOPMXF , NDVMXF , NOPMXN , NOVMXN , NOPSSP , NOEXPR , NDCOMV , 

«  NOTBiP,  NOZROe,  N0PAT8,  NTMFVa .  NTBPl ,  NTaiP2 .  NTBIPS 

COMUON/SVALb/NSRB:,NSEqNS(MXSRH:)  ,SUL(MXSREQ)  ,SUU(MXSREG) , 

*  SVL (MXSREG) , SVU (MXSREG) , NSU (HXSREG) , NSV (UXSREQ) , 

*  CCS(10, MXSREG) , TRNSS (3, MXSREG) , IZAX(MXSREG) , 

«  IOROU(MXSREG),IORDV(MXSREG),NCCEIIS 

CCMMON/aWINP/AINC ,  THTINC ,  PHIINC ,  ISCATR 
COMMN/BAFFLE/INFRIG 
COIMON/MODAL/IMOOAL,  lUSaq 


RHO=1000. 

C=1600. 


C««««  <••«*««««  changed  NCCEI^S  to  13  to  add  cone  -  19  Sep  1990  «•••««•••« 

C  NUMBER  OF  EQUATION  SETS  IN  CCUNMD 
NCCE1)S=13 

Cfa**********************************************^*********************'' 

C  UNIT  DESIGNATORS  FOR  MASS  STORAGE 
NOQPr  -10 
NOPMXS=ll 
NDVMXS^IZ 
NOOECM=13 
NDVaS=14 
NDSPS=16 
N0PMXF=1B 
N0VMXF=17 
N0PMXN=18 
NDVMXN>ig 
NDPSSP=20 
NOEXPft=21 
NDC0MV=22 
NDTBff>^23 
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N0ZRDB=24 

C  3  Jan  1991  -  jdded  teaporary  logical  unit  numbers  TEUPVEL,  TEMPI 

NTMPVa  =  26 
NTEMPl  ^  V 
NTEMP2  =  28 
NTEMP3  =  29 

C  PRINT  UNIT  DESIGNATORS 

NUNPRT=:6 
NUNEJRR=6 

C  aANE  lAVt  va;_ues 

AiNC^0  8 

isca:r=0 

C  INPINITE  rigid  BAFPLE  VALUE 

:NFRiG=e 

C  MODAL  MATRIX  CQMP'JTATION  FUG 

IMOOA.^-  id 
iJSB-au 


re;  URN 
EiNU 


-  TEMP3. 
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